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Atmospheric  Transmittance/Radiance: 

Computer  Code  Lowtran  6 


1.  INTRODUCTION 


This  report  describes  a  Fortran  computer  code,  LOWTRAN  G,  designed  to 
calculate  atmospheric  transmittance  and  radiance  for  a  given  atmospheric  path 
at  moderate  spectral  resolution.  This  code  is  an  extension  of  the  current 
LOWTRAN  atmospheric  code,  LOWTRAN  5*  (and  its  predecessors  LOWTRAN  4,“ 
LOWTRAN  3B,  3  LOWTRAN  3, 4  and  LOWTRAN  25).  All  the  options  and  capabil¬ 
ities  of  the  LOWTRAN  5  code  have  been  retained.  Solar/lunar  scattered  radiation 
(Received  for  publication  2 0  July  1983) 

1.  Kneizys,  P.X.,  Shet.tle,  E.  P.  ,  Gallery.  W.O.  ,  Chetwynd,  Jr.,  J.  II.  , 

Abreu,  L.  W.  ,  Selby,  J.E.A.,  Penn,  R.W.,  and  McClatchey,  R.A.  (198O) 

A tmo sphe r i c  Trans m i ttance / Rad la nee:  Computer  Code  LOWTRAN  5 , 
AFGL-TR  -80-0067 ,  "AD  A0882  15. ' 

2.  Selby,  J.E.A.,  Kneizys,  F.X.,  Chetwynd,  Jr.  ,  J.  II.  ,  and  McClatchey,  R.A. 

(1978)  Atmospheric  Transmittance/Radiancc :  Computer  Code  LOWTRAN  4, 
AFGL-TR -  7 8-0033 ,  AD  A058643. 

3.  Selby,  J.E.A.,  Shettle,  E,  P.  ,  and  McClatchey,  R.A.  (1970)  Atmospheric 

Transmittance  from  0.2b  to  28.  5  am:  Supplement  LOWTRAN  3R , 

AFGL-TR -78-0258,  AD  A040701. 

4.  Selby,  J.E.A.,  and  McClatchey,  R.A.  (1978)  Atmospheric  Transmittance 

from  0.2  5  to  2  8.  5  urn:  Computer  C  ode  L0WTITAN_3 ,  AFCR  L-TR  -7 8-0 2 5 5 , 
"AD  Ao  17734. 

5.  Selby,  J.L.A.,  and  McClatchey,  R.A.  (1972)  Atmo:  phcric  Transmittance 

from  0.25  to  28.  5  nm:  Computer  Code  LOWTRAN  9,  AFCR  L-TR  -72 -074a, 
ADA763721. 
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has  been  added  to  the  code,  ns  well  ns  new  spherical  refractive  geometry  sub¬ 
routines  and  an  improved  water  vapor  continuum  model.  Other  new  features  of 
the  code  include  a  wind-dependent  maritime  aerosol  model,  a  vertical  structure 
aerosol  model,  a  cirrus  cloud  model,  and  a  rain  model. 

The  LOWTRAN  code  calculates  atmospheric  transmittance  and  radiance, 
averaged  over  20-cin  ^  intervals  in  steps  of  5  cm  ^  from  330  to  40,000  cm"'' 

(0.  2  5  to  28.  5  Mm).  The  code  uses  a  single-parameter  band  model  for  molecular 
absorption,  and  includes  the  effects  of  continuum  absorption,  molecular  scatter¬ 
ing,  and  aerosol  extinction.  Refraction  and  earth  curvature  are  Included  in  the 
calculation  of  an  atmospheric  slant  path.  The  code  contains  representative 
atmospheric  and  aerosol  models,  and  the  option  to  i  epiace  them  with  user-derived 
or  measured  values. 


In  this  report,  the  new  spherical  refractive  geometry  subroutine  is  described 
in  Section  2.  The  improved  water  vapor  continuum  model  is  presented  in  Section 
3.  Following  this  is  a  discussion  of  the  calculation  of  single -scattered  solar  or 
lunar  radiation  in  Section  4.  In  Sections  5  through  8,  the  wind-dependent  maritime 
model,  the  vertical  structure  algorithm  for  boundary  layer  aerosols,  a  cirrus 
cloud  model,  and  a  rain  model  are  presented.  The  structure  of  the  new  computer 
code  is  given  in  Section  9. 

User  instructions  for  the  LOWTRAN  code  arc  given  in  Section  10  and  exam¬ 
ples  of  the  output  cf  the  program  in  Section  11.  A  plotting  program  and  a  filter 
program  are  now  included  with  the  LOWTRAN  code,  and  are  explained  in  Appen¬ 
dices  A  and  13.  Appendix  C  defines  the  single  scattering  geometry  while  Appendix 
D  explains  the  development  ot'  the  standard  LOWTRAN  phase  functions.  The  pro¬ 
gram  listing  along  with  the  plot  program  and  filter  program  listing  are  available 
in  a  supplement  to  the  LOWTRAN  0  report. 

In  addition  to  the  modeling  work  of  the  Air  Force  Geophysics  Laboratory, 
special  acknowledgement  should  be  given  to  the  contributions  made  by  Sonicraft 
Incorporated,  the  Naval  Research  Laboratory,  the  Army  Atmospheric  Sciences 
Laboratory,  and  the  National  Oceanic  and  Atmospheric  Administration  to  the 
LOWTRAN  6  model  and  code  presented  in  this  report. 

The  LOWTRAN  6  code  will  be  made  available  from  the  National  Climatic 


Center,  Federal  Building,  Asheville,  NC  28801  (see  Section  9  for  complete 
details  on  availability).  It  is  requested  that  users  receiving  the  code,  send  their 
name,  affiliation,  and  address  to  us  to  update  the  AFCL  LOWTRAN  mailing  list 
and  for  notification  to  users  of  changes  in  the  code.  Correspondence  should  be 
mailed  to  F.X.  Kneizys,  AFGL/OPI,  Ilanseom  AFJ3,  MA  01731, 
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2.  AIR  MASS  COMPUTATION  (SPHERICAL  REFRACTIVE  GEOMETRY) 


2.1  Introduction 


The  transmittance  and  radiance  along  a  path  through  the  atmosphere  depend 
upon  the  total  amount  and  the  distribution  of  the  absorbing  or  scattering  species 
along  the  path.  The  integrated  amount  along  a  path  is  known  by  various  names, 
including  "column  density",  "equivalent  absorber  amount",  and  "air  mass". 

While  the  term  "air  mass"  applies  specifically  to  the  total  amount  of  gas  along 
tlie  path,  ic  will  be  used  here  to  refer  loosely  to  the  integrated  amounts  for  all 
the  different  species  relative  to  the  amount  for  a  vertical  path.  The  calculation 
of  air  mass  for  realistic  atmospheric  paths  requires  that  the  eai'th's  curvature 
and  refraction  be  taken  into  account. 

The  model  for  calculating  air  mass  has  been  greatly  improved  in  LOWTRAN 
G.  The  previous  model  assumed  that  the  index  of  refraction  was  constant  between 
layer  boundaries.  The  new  model  assumes  a  continuous  profile  for  the  refractive 
index,  with  an  exponential  profile  between  layer  boundaries.  It  is  more  accurate 
than  the  previous  model  and  works  for  all  paths.  All  the  options  from  the  previ¬ 
ous  model  for  specifying  the  slant  paths  have  been  retained. 

This  section  describes  the  new  model  for  calculating  air  mass  and  presents 
calculations  of  air  mass  for  several  representative  atmospheric  paths.  Tor  a 
more  complete  description  of  the  method  used  here,  see  Reference  G. 


2.2  Defining  Expiations 

The  atmosphere  is  modeled  as  a  set  of  spherically  symmetric  shells  with 
boundaries  at  the  altitudes  z^,  j  =  1,  N.  The  temperature,  pressure,  and  absorb¬ 
er  (gas  and  aerosol)  densities  are  specified  at  the  layer  boundaries.  Between 
boundaries,  the  temperature  profile  is  assumed  linear  while  the  pressure  and 
densities  are  assumed  to  follow  exponential  profiles.  For  example,  the  density 
p  at  an  altitude  z  between  z.  and  +  ^  is  given  by 

p(z)  =  pj  exp  [-(z-zA/IIp] 
where  the  density  scale  height  H  is 

"p  '  Vi "  V/ln  'TVi’  ■ 

The  scale  height  varies  with  the  layer  and  is  different  for  pressure  and  density. 

G.  Gallery,  W.O.,  Kneizys,  F.X.,  and  Clough,  S.A.  (1983)  Air  Mass  Computer 
Program  for  Atmospheric  Transmittance  Radiance/ Calculations:  FSCATM, 
AFG.L-TR  -83  -GOG  !5 . 
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Consider  an  optical  path  through  the  atmosphere  from  point  a  to  b  as  shown 
in  Figure  1.  The  path  is  defined  by  the  initial  and  final  altitudes  z  and  z^  and 
by  the  zenith  angle  at  a.  The  other  path  quantities  are:  k,  the  curved  path 
length  from  a;  p,  the  earth-centered  angle;  p,  the  zenith  angle  at  b;  and  </,,  the 
total  refractive  bending  along  the  path. 


Figure  1.  Slant  Path  Through  the  Atmosphere  From 
Point  a  to  Point  b 


The  integrated  amount  u  of  an  absorber  of  density  p(z)  is  given  by: 


a 


(1) 


(2) 
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At  any  point  along  the  path 


ds/dz  =  (cos  9)'1  ,  (3) 

where  6  is  the  local  zenith  angle. 

Due  to  the  curvature  of  the  earth  and  to  refraction,  Q  varies  along  the  path. 
However,  if  the  zenith  angle  is  less  than  about  80°,  the  variation  of  6  along  the 
path  is  negligible  and  Eq.  (2)  can  be  written  as 

b 

u  =  (cos  y0)_1  f  P  dz  .  (4) 

a 

Eq.  (4)  is  called  the  secant  approximation  and  is  equivalent  to  assuming  a  plane- 
parallel  atmosphere.  The  integral  in  Eq.  (4)  has  a  particularly  simple  form  for 
an  exponential  density  distribution: 

b 

f  P  dz  =  Hp  [p(za)  -p(zb)]  ,  (5) 

a 

where  Hp  is  the  density  scale  height.  If  the  path  extends  over  several  layers, 
each  with  a  different  scale  height,  then  the  integral  in  Eq.  (5)  must  be  broken 
into  separate  parts,  one  for  each  layer. 

For  the  general  case,  curvature  and  refraction  must  be  taken  into  account 
in  Eq,  (2).  This  is  accomplished  by  a  detailed  numerical  integration  of  Eq.  (2) 
as  follows.  The  interval  from  za  to  z b  is  divided  into  a  number  of  sub-intervals 
defined  by  z^,  z^,  .  .  .  z^.  The  integral  in  Eq.  (2)  is  approximated  by  the  sum: 


N- 1 

u  =  E  ^si  '  (6) 

i=  1 


where 


As. 

1 


(8) 


z. 


i+1 


/ 


z. 

1 


(ds/dz)  dz 


Since  the  density  is  assumed  to  follow  an  exponential  profile,  the  integral  in 
Eq.  (7)  can  be  written  analytically  as 


H 

p.  = 

Az. 

i 


[p(z.)  -  p(zi+i)] 


(9) 


where  the  scale  height  is  constant  over  the  layer  from  z^  to  z^.  The  integral 
in  Eq,  (8)  can  be  obtained  numerically  as  shown  in  the  next  section. 

The  number  and  spacing  of  the  intervals  z.  are  chosen  so  that  Eq.  (6)  is  a 
good  approximation  to  Eq.  (2)  as  will  also  be  shown.  Again,  if  the  path  extends 
over  several  layers,  with  different  scale  heights  in  each,  then  the  path  integral 
must  be  performed  separately  for  each  layer.  In  the  discussion  that  follows,  it 
will  be  assumed  that  the  path  is  confined  to  a  single  layer  in  which  the  scale 
heights  are  constant  with  altitude. 

2.3  Atmospheric  Refraction 

The  governing  equation  for  a  ray  passing  through  the  atmosphere  is  Snell's 
Law  for  a  spherically  symmetric  medium,  given  by: 

n(r)  r  sin  9  =  C  ,  (10) 

where  n  is  the  index  of  refraction,  r  is  the  radius  to  a  point  along  the  ray,  9  is 
the  zenith  angle  at  that  point,  and  C  is  a  constant  of  the  particular  path.  If  the 
ray  is  horizontal  at  a  point  ivp,  0  is  -  90°  at  ilia t  point,  and  C  -  n(iVp)r  j,;  the 
altitude  at  that  point  is  called  the  tangent  height. 

The  index  of  refraction  n  is  conveniently  written  as: 

n(r)  =  1  +  N  (r)  ,  (11) 

where  N(r)  is  called  the  refractivity.  *  N  is  wavenumber  dependent  and,  in  the 

visible  and  the  infrared,  N  is  also  very  nearly  proportional  to  the  total  air 

-4 

density.  At  sea  level,  in  the  infrared,  N  is  of  the  order  of  3  x  10  .  In  these 


*See  Section  2. 6  for  a  discussion  of  the  index  of  refraction. 
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calculations ,  we  assume  the  N  follows  an  exponential  profile  with  a  scale  height 
is  determined  separately  for  each  atmospheric  layer. 

The  effect  of  refraction  is  to  bend  the  path  in  the  direction  of  increasing  N. 
The  radius  of  curvature  K  of  the  refracted  ray  can  be  shown  to  be 

K  =  ~(n'/ n)  sin  0  ,  (12) 

where  n'  =  dn/dr.  It  is  useful  to  define  the  quantity  R(r)  as 


R(r) 


(13) 


JR  is  simply  the  ratio  of  r  to  the  radius  of  curvature  of  a  ray  tangent  at  r .  R  is  a 
property  of  the  atmospheric  profile  and  not  the  particular  path,  and  is  a  good 
measure  of  the  importance  of  refraction  at  a  particular  altitude.  For  example, 
for  the  U.S.  Standard  Atmosphere,  R  is  approximately  0.  16  at  sea  level  and 
decreases  exponentially  with  altitude  with  a  scale  height  of  about  10  km. 

To  trace  a  ray  through  the  atmosphere,  consider  the  path  shown  in  Figure  l: 
9  is  the  zenith  angle  at  z  ,  <t>  is  the  zenith  angle  at  z  ,  ft  is  the  earth -centered 
angle,  and  ip  is  the  bending  along  the  path.  Let  s  be  the  length  of  the  path  from 
point  a.  At  any  point  the  differential  path  quantities  are  given  by: 


ds 

=  n  dr 

cos  0 

(14) 

d  13 

=  tan  0  —  dr  , 
r 

(15) 

where  0  is  the  zenith  angle  at  the  point.  Substituting  for  cos  9  from  Fq.  (10)  into 
Eq.  (14)  gives 


ds  = 


1  - 


C 


.2 


-1/5 


2  2 
n  r 


dr 


(1C) 


Eq.  (16)  is  the  basic  atmospheric  ray  trace  equation.  If  the  function  n(r)  is 
known,  then  Eq.  (16)  can  be  integrated  numerically  along  the  path. 

However,  the  difficulty  with  integx  ating  Eq.  ( 1  (j)  is  that  it  has  a  singularity 
at  0  =  90°,  that  is,  at  the  tangent  height,  where  C  =  n(i  ,p)r,^,.  A  simple  change 
of  variables  will  remove  this  singularity  and  also  provide  some  insight  into  the 
importance  of  refraction.  Define  a  new  independent  variable  x  as: 
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x  =  r  cos  6 


(17) 


(x  can  be  interpreted  as  the  straight  line  distance  to  the  geometric  tangent  point). 
Differentiating  Eq.  (17)  gives: 

dx  =  [cos  6  -  r  sin  6  (d0/dr)]  dr  .  (18) 

Differentiating  Eq.  (10)  and  using  Eq.  (13)  gives 

dO/dr  =  -(l  -  R)  tan  0/r  .  (19) 

Substituting  Eq,  (19)  into  Eq.  (18)  gives 

dx  =  (l  -  R  sin2  6)  dr/ cos  6  ,  (20) 

Comparing  Eq.  (20)  with  Eq.  (14)  gives 

ds  =  (1  -  R  sin2  0)  1  dx  .  (21) 

In  this  form  of  the  equation  for  ds,  the  right-hand  side  is  a  well-behaved  function 
of  r  for  all  paths,  including  vertical  and  horizontal  paths  (except  in  the  unusual 
case  where  R  >  1  and  the  path  curves  back  toward  the  earth,  that  is,  looming). 
The  intermediate  variable  x  =  r  cos  0  is  also  well  defined  for  all  paths.  In 
practice,  the  numerical  integration  of  Eq.  (2l)  is  driven  in  steps  of  r,  from  r  to 
r  +  Ar.  The  corresponding  increment  in  x  is  calculated  from  Eq.  (17).  The 
integration  of  s  from  Eq.  (2l)  is  then  straightforward. 


2.4  Numerical  Algorithm 


The  numerical  algorithm  used  to  evaluate  Eq.  (6)  is  as  follows: 

1.  Find  the  minimum  and  maximum  altitude  HMIN  and  HMAX  along  the  path 
and  the  zenith  angle  0  at  HMIN.  If  the  path  goes  through  a  tangent  point,  then 
solve  Eq.  (10)  iteratively  for  the  tangent  height. 

2.  From  the  given  atmospheric  profile,  construct  a  new  profile  at  the  layer 
boundaries  from  HMIN  to  HMAX,  interpolating  the  pressure,  temperature,  and 
densities  where  necessary, 

3.  Starting  with  the  lowest  layer,  trace  the  path  through  each  layer: 

a.  Divide  the  layer  into  sublayers  defined  by  the  altitudes  z.  such  that 

■X  —X-  J 

Az^  =  As  cos  CF  _  ^ ,  where  As  is  a  nominal  path  length  t5  km)  and  F  _  ^  is  the 
zenith  angle  at  zj.j. 


1(1 


Az.  =  As  cos  0.  , 

3  J-l 


z.  =  z.  ,  +  Az. 
3  3"1  3 


r-  =  r  +  z.  (r  is  radius  of  the  earth) 
J  “  3  ® 


n.  =  1  +  N(z.) 

J  J 


sin  0.  = 


C 


1  n.r 
J  3  3 


2  l/2 

cos  0.  =  (1  -  sin  0.) 

3  J 


x.  -  r.  cos  0. 
3  3  3 


Ax.  =  x.  -  x.  , 
3  3  3  -1 


R,  =  -r,/[(dN/dr |  ,)/n.l 
v  J  J  J 


..•2  ,-l 


ds/dx|  ,  =  (1  -  R.  sin  0.) 
J  3  3 


ASj  =  l/2  (ds/dx|^_^  +  ds/dx|j)  Axj 


t.  For  each  species,  integrate  the  density  p 


pj  2  p(zj> 


pi  =  hp  <pj ' 


N-l 

u  ’  S  pj  asj  • 


3  =  1 
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2.5  Examples  of  Air  Mass  Calculation 


This  section  will  present  plots  of  air  mass  values  for  three  classes  of  slant 
paths.  The  term  "air  mass  values"  refers  to  the  integrated  amount  of  air  along 
a  path  compared  to  the  amount  for  a  vertical  path  from  ground  to  space.  For 
example,  the  air  mass  value  for  a  path  from  the  ground  to  space  with  the  zenith 
angle  at  the  ground  of  90°  is  38.  1  for  the  U,  S.  Standard  Atmosphere  {one  air 
mass  equals  2.  15  X  102°  molecules  cm  2  or  1.034  x  102  gm  cm”2).  These  three 
paths  are  described  by  the  initial  altitude  "Hi"  and  the  zenith  angle  "ANGLE"  at 
Hi.  The  other  end  of  the  path  is  the  top  of  the  atmosphere,  here  taken  to  be 
100  km.  The  three  classes  of  paths  are:  1)  Hi  =  0  km  for  ANGLE  varying  from 
0  to  90°;  2)  ANGLE  =  90°  for  Hi  varying  from  0  to  50  km;  3)  111  =  30  km  for 
ANGLE  varying  from  85°  to  95.  1°  at  which  point  the  path  intersects  the  earth. 
The  wavenumber  for  these  calculations  was  2000  cm  1  (5  mm).  The  dependence 
of  air  mass  on  wavenumber  in  the  infrared  is  small. 

In  addition  to  the  air  mass  value,  the  amounts  of  water  vapor  and  of  ozone 
relative  to  the  amounts  for  a  vertical  path  from  ground  to  space  are  also  shown. 
The  relative  amounts  of  these  gases  depend  upon  their  vertical  distribution;  in 
these  cases,  the  1962  U.S.  Standard  Atmosphere  density  profiles  aroused. 

Since  the  distributions  of  these  gases  in  the  atmosphere  are  so  variable,  the 
relative  amounts  for  other  profiles  could  be  significantly  different.  The  values 
for  water  vapor  and  ozone  shown  here  should  be  taken  to  be  illustrative  only. 

Figure  2  (a  and  b)  shews  the  air  mass  value,  relative  water  vapor,  and 
ozone  amounts  for  path  1.  Also  shown  in  Figure  2(b)  is  the  secant  of  the  zenith 
angle.  B’or  a  large  zenith  angle,  the  relative  amount  of  water  vapor  is  greater 
than  the  air  mass  value  while  the  amount  of  ozone  is  less.  This  effect  is  due  to 
the  fact  that  for  large  zenith  angles,  the  greater  part  of  the  path  is  near  the 
ground.  Water  vapor  is  concentrated  in  the  lower  layers  so  that  the  relative 
amount  of  water  vapor  is  large  compared  to  the  vertical  path.  Ozone,  however, 
is  concentrated  in  the  stratosphere,  which  contains  a  relatively  small  part  of  the 
path.  Note  that  the  secant  agrees  to  better  than  one  percent  with  the  air  mass 
value  up  to  72°,  up  to  80°  for  water  vapor,  but  only  up  to  60°  for  ozone.  The 
discrepancy  is  due  mainly  to  the  effect  of  the  earth's  curvature  and  not  refrac¬ 


tion;  by  including  curvature  but  neglecting  refraction,  the  relative  amounts  can 
be  calculated  to  better  than  one  percent  up  to  84°  for  air,  86°  for  water  vapor, 
and  82°  for  ozone. 

Figure  3  shows  the  air  mass  values  and  relative  amounts  for  path  2.  These 
curves  mimic  the  density  profiles  of  air,  water  vapor,  and  ozone  themselves 
since  the  bulk  of  the  gas  is  located  within  a  few  kilometers  (vertically)  of  the 
observer  altitude. 
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Figure  2.  Relative  Absorber  Amounts  vs  Zenith  Angle  for  Path  1. 
<a)  0  to  90°  and  (b)  74  to  90°  (also  shown  is  the  secant  of  the  zenith 
angle) 


Figure  3.  Relative  Absorber  Amounts  vs 
Observer  Altitude  (Ul)  for  Path  2 
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The  relative  amounts  shown  in  Figure  4  correspond  to  path  3,  which  is 
typical  of  a  stratospheric  balloon-borne  experiment  looking  at  the  sun  as  it  sets. 
Also  shown  on  the  right-hand  axis  is  the  tangent  height  vs  zenith  angle  and  the 
angular  diameter  of  the  sun.  If  the  sun  is  used  as  the  source  for  a  measurement, 
the  ail'  mass  value  io  different  points  on  the  face  of  the  sun  can  vary  by  a  factor 
of  2  for  large  zenith  angles.  The  variation  in  air  mass  due  to  this  effect  can  be 
a  major  source  of  uncertainty  in  the  measurement  and  must  be  considered 
carefully. 


ZENITH  ANGLE  (DEG) 


Figure  4.  ReLat.ive  Absorber  Amounts  vs  Zenith 
Angle  for  Path  3.  Also  shown  against  the  right- 
hand  axis  is  the  tangent  height  vs  zenith  angle  and 
the  angular  diameter  of  the  sun 


Two  other  quantities  of  interest  for  atmospheric  profiles  arc  the  tangent 
height  and  the  refractive  bending.  The  difference  in  tangent  height  between  an 
unrefracted  and  a  refracted  ray  coming  in  from  space  is  shown  as  a  function  of 
the  refracted  tangent  height  in  Figure  5  for  three  atmospheric  profiles  (the 
geometry  is  shown  schematically  in  the  inset).  The  total  refractive  bending  for 
paths  1  and  2  are  shown  in  Figures  (i  and  7  for  three  atmospheric  profiles.  Note 
that  the  total  bending  for  a  path  from  the  ground  to  space  at  90°  for  the 
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BENDING  (DEG 3 


^  % 


Figure  5.  Unrefracted  Tangent  Height  Minus 
Refracted  Tangent  Height  vs  Refracted  Tan¬ 
gent  Height  for  a  Ray  Coming  in  From  Space 
for  Three  Atmospheric  Profiles.  The  figure 
in  the  inset  illustrates  the  paths 


Figure  (j.  Refractive  Bending  vs  Zenith  Angle  for  Patli  1,  for  Three 
Atmospheric  Profiles,  (a)  0  to  90°  and  (b)  74  to  90° 
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OBSERVER  ALTITUDE  (km) 


Figure  7.  Refractive  Bending  vs  Observer 
Altitude  (Hi)  for  Path  2,  for  Three  Atmos¬ 
pheric  Profiles 


U.S.  Standard  Atmosphere  and  the  Tropical  Atmosphere  is  about  0.  5°,  which 
is  the  same  as  the  solar  diameter. 


r  .1  P 

,  /  1  W 

+  co  •  F"  ' 

L  ->  o 

where  v  is  the  wavenumber  in  cm  P  is  the  total  pressure  in  mb,  P  is  the 
partial  pressure  of  water  vapor,  PQ  is  1013.25  mb,  T  is  the  temperature  in 
Kelvin,  and  the  constants  a,  b,  and  c  are 

7.  Edlen,  K.  (1966)  The  refractive  index  of  air,  Metrologia  2:12, 
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ao  =  83.  43,  ax  =  185.  08,  a2  =4.11 

b1  =  1.  140  X  lO5,  b.,  =  G .  24  X  104 

c  =  43.  49,  c  =  1.  70  X  104  . 

O  1 


The  formula  used  in  previous  versions  of  LOW  THAN  was  a  simplified  version  of 
this  expression. 


3.  WATLK  VAPOR  CONTINUUM 


A  new  water  vapor  continuum  model  has  been  added  to  LOW  TRAN  G.  This 

model  for  the  continuum  contribution  from  water  vapor  absorption  was  originally 

8 

developed  by  Clough  et  al  l'ox-  use  with  the  line-by-line  transmittance  and  radi¬ 
ance  atmospheric  code,  FAS  CODE ,  ® 

For  atmospheric  applications  it  is  advantageous  to  express  the  density 
dependence  of  the  water  vapor  continuum  absorption  in  terms  of  a  self  and  foreign 
component.  The  continuum  contribution  to  the  absorption  coefficient  k^(r),  is 
given  by  the  expression 


kr(r)  =  p  v  tanh(hcv/2kT) 


C>,T)  + 

to 


(22) 


where  T  is  the  temperature  <°K),  v  the  wavenumber  (cm  4),  hc/k  =  1.4387  9 
°K/cra"  (p  /p  )  and  (pf/p  )  are  the  number  density  ratios  for  the  self  and  foreign 
continuum;  and  C  and  | (cm  1  mol/ cm")  ]  are  wavenumber  dependent  contin¬ 
uum  absorption  parameters  for  the  self  and  foreign  components.  The  density  Ps 
is  the  density  of  the  water  vapor  and  p^.  is  the  density  of  all  other  molecular 
species;  consequently,  p  +  p^.  represents  the  total  density.  The  quantity,  p^,  is 
the  reference  number  density  defined  at  1013  mb  and  2  9RK.  The  present  formula¬ 
tion  in  terms  of  density  has  the  advantage  that  the  continuum  contribution  to  the 
absorption  coefficient  decreases  with  increasing  temperature  through  the  number 


8.  Clough,  S.A.,  Kneizys,  F.X„,  Davies,  R. ,  Gamache,  R.  ,  and  Tipping,  R.II. 
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Atmospheric  Water  Vapor,  A.  Deepak,  T.D.  Wilkerson,  and  L.  II,  Rulmke, 
Eds. ,  Academic  Press,  New  York. 

9.  Clough,  S.A.,  Kneizys,  F.X.,  Rothman,  L.S.,  and  Gallery,  W.O.  ( 1  p 8 1 ) 

Atmospheric  spectral  transmittance  and  radiance:  FAS  CODE  lB,  Pro¬ 
ceedings  of  SFIE,  The  Inter.  Soc,  for  Opt.  Eng.  ,  277,  Atmospheric 
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density  ratio  term.  The  quantities  C  and  Cf  ioi'  water  vapor  are  stored  m  the 
program  for  the  spectral  range  0  to  20,000  cm 

The  values  for  C'  for  water  vapor  at  2  9GK  are  shown  in  l/igure  8  together 
with  the  experimental  values  obtained  by  Burch  et  al.  ^  The  strong  tempera 
lure  dependence  of  the  self  density  dependent  water  vapor  continuum  is  treated 
by  storing  values  of  ( ' ,  al  2G0K  and  2yt>K  and  linearly  interpolating  between  the 
2C0K  and  2PtiK  values.  The  2 0 0 K  j'esult  was  obtained  by  extrapolating  the  fits  to 

i  ■> 

the  338K  and  2  pfiK  data  of  Bur  eh  et  al.  "  The  results  for  2  GQK  and  2  9GK  are 
Shown  in  Figure  9. 


WAVENUMBER  (CM- 1 ) 

Figure  8.  The  Self  Density  Dependent  Continuum  Values, 
for  Water  Vapor  as  a  Function  oi  Wavenumber.  The  exper -  J 
imental  values  are  from  Burch,  et  al 
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Figure  9.  The  Self  Density  Dependent  Continuum  Values,  C  , 
for  Water  Vapor  as  a  Function  of  Wavenumber  at  260K  and 
296K,  The  values  from  296K  are  fits  to  experimental  results; 
the  2  60K  result  is  extrapolated 


Only  values  near  room  temperature  are  available  for  the  foreign  dependence 
of  the  water  vapor  continuum.  The  continuum  values  at  2  9GK  are  shown  in 
Figure  10  ^:,d  have  been  obtained  by  a  fit  to  the  data  of  Burch.  There  is  still 
considerable  uncertainty  in  the  foreign  values  for  the  spectral  wind  w  regions  at 
1000  and  2500  cm 

In  the  LOWTRAN  code,  the  total  optical  depth  due  to  water  vapor  continuum 
absorption  for  an  atmospheric  slant  path  of  N  layers  is  given  by 
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WAVENUMBER  (CM- 1) 


figure  10.  The  Foreign  Density  Dependent  Continuum  Values, 
C^,  for  Water  Vapor  as  a  Function  of  Wavenumber.  The  ex¬ 
perimental  values  are  from  Burch,  et  al 


where  ds  is  the  incremental  path  length,  T.  is  the  temperature  of  the  i'th  layer, 
and 


C  (*,296)  =  v  tanh  (hcv/2k(2  96))  £  (v, 2  96) 

s 

C  (*,260)  =  *  tanh  <hc*/2k(2G0))  £  (*,260)  (24) 

Cf  (*,296)  =  *  tanh  (hc*/2k(296))  Cf  (*,296) 

Calculations  of  atmospheric  slant  path  transmittance  using  the  new  water 
vapor  continuum  absorption  coefficients  will  result  in  approximately  the  same 
attenuation  as  in  the  LOWTRAN  5  model  for  the  atmospheric  window  regions  from 
8  to  12  Mm  and  3.  5  to  4.2  Mm.  However,  for  other  spectral  regions,  particularly 
from  4.  5  to  5.0  Mm,  significant  improvement  in  atmospheric  transmittance  calcu¬ 
lations  has  been  made  with  the  inclusion  of  the  contribution  of  continuum  absorp- 


4.  SOLAR/LUNAR  SINGLE  SCATTERING  MODEL 
W.L.  Ridgway,  R.A.  Moose,  A.C.  Cogley 
Sonicraft,  Iiirorporated,  Chicago,  Illinois 

4.1  Introduction 

The  radiation  propagating  through  the  atmosphere  originates  from  the  follow¬ 
ing  sources:  gaseous  emission  along  the  line-of-sight,  transmitted  extraterres¬ 
trial  (solar/lunar)  sources  or  background  emission  (earth  or  target),  and  radia¬ 
tion  scattered  into  the  line-of-sight  by  aerosols  or  molecules. 

Previous  versions  of  LOWTRAN  have  calculated  atmospheric  radiance  due 
only  to  gaseous  and  background  emissions.  While  in  many  cases  these  two 
sources  dominate  the  atmospheric  radiance,  there  are  other  cases  of  interest 
where  the  scattered  radiance  is  of  equal  or  greater  importance.  Until  now 
LOWTRAN  has  treated  scattering  only  as  a  loss  mechanism. 

A  number  of  techniques  exist  that  include  the  full  effects  of  scattering  on 
atmospheric  radiance:  these  include  for  example,  Monte  Carlo  and  "adding/ 
doubling"  techniques.  These  techniques  take  into  account  multiple  scattering 
and  can  include  both  external  and  internal  sources.  These  techniques  however, 
tend  to  be  computationally  expensive  and  some  of  them  are  incompatible  with  the 
spectrally-averaged  band  model  used  in  LOWTRAN. 

In  many  situations,  a  complete  multiple  scattering  calculation  is  not  neces¬ 
sary  and  the  scattered  radiation  is  dominated  by  solar  radiation  that  has  been 
scattered  only  once.  Calculation  of  the  single  solar  (or  lunar)  scattered  radiation 
is  relatively  simple  and  fits  well  within  the  context  and  structure  of  LOWTRAN. 

Calculation  of  single  scattering  has  been  incorporated  as  an  option  in 
LOWTRAN  G.  The  next  two  sections  of  this  chapter  develop  the  algorithm  for 
single  scattering  and  show  the  verification  of  the  LOWTRAN  6  calculations  against 
other  methods.  Then  the  phase  functions  for  scattering  by  atmospheric  aerosols 
and  molecules  will  be  explained.  Next,  sample  calculations  of  solar  scattering 
are  shown  that  illustrate  the  conditions  where  the  singly  scattered  radiation 
becomes  significant  compared  to  the  emitted  radiation.  Finally,  recommenda¬ 
tions  are  given  concerning  the  range  of  applicability  of  the  single  scattering 
model.  For  a  more  detailed  discussion  of  the  single  scattering  model,  see 
Reference  14. 


14.  Ridgway,  WL  L.  ,  Moose,  R.A.,  and  Cogley,  A.  C.  (1982)  Single  and  Multiple 
Scattered  Solar  Radiation,  AFGL-TR -82 -02 99,  AD  A12G323, 


4.2  Radiative  Transfer 


Before  proceeding  further,  it  will  be  helpful  to  introduce  the  following 
nomenclature : 

SUPERSCRIPTS: 

A  aerosol 

M  molecular 


SUBSCRIPTS: 

e 

a 

s 

ps 

op 


extinction 

absorption 

scattering 

primary  solar  path  (sun  to  scattering  point) 
line-of-sight  optical  path  (scattering  point  to  observer) 


OTHER  QUANTITIES: 


k 


f(y) 

XSUN 

tMOON 


monochromatic  volumetric  extinction,  absorption,  or 
scattering  coefficient 

monochromatic  transmittance  over  a  homogeneous  path  length 
s  due  to  extinction  absorption,  or  scattering 

scattering  phase  function  for  a  scattering  angle  7 

solar  extraterrestrial  intensity 

lunar  extraterrestrial  intensity 


Note  that  the  dependence  of  most  quantities  on  the  spectral  frequency  v  will  be 
shown  by  a  subscript  v,  although  it  will  sometimes  be  suppressed  for  simplicity 
of  notation  when  the  concept  is  clear  from  context. 

The  monochromatic  intensity  (radiance)  seen  by  an  observer  looking  along  a 
particular  directional  path  is  the  sum  of  contributions  from  all  sources  lying 
along  the  line-of-sight.  The  sources  are  either  primary  sources  (infrared  emis¬ 
sion)  or  scattering  sources.  The  scattering  source  function  J  for  scattering 
points  along  the  observer's  line-of-sight  can  be  expressed  in  terms  of  the  local 
incoming  intensity  at  each  point  I  (n1)  by 


«)  --  f  1  tt')fpAkA  +  P«k«') 

1  J  V  \  S  VS  j 


dn 


(2  5) 


where  n  is  the  unit  vector  directed  toward  the:  observer  and  n1  (SI1)  is  to  be  inte¬ 
grated  over  the  solid  angle  denoted  by  Si1.  With  only  single  solar/lunar  scattering 
included,  the  incident  intensity  I  (n1)  is  given  by 


I 

v 


in')  -  I 


SUN 

v 


ta+m 

e,  ps 


6(n\  ns’> 


(2  5) 


2  8 


•where  ng'  is  the  direction  of  the  incident  solar/lunar  radiation  at  the  scattering 
point.  A  schematic  of  the  scattering  geometry  for  a  particular  sun/observer 
orientation  is  displayed  in  Figure  11.  The  path  that  the  sunlight/moonlight  takes 
in  passing  through  the  atmosphere  prior  to  being  scattered  at  any  scattering  point 
P  will  be  called  the  'primary  solar'  path.  Other  sources  besides  direct  extra¬ 
terrestrial  illumination  could  of  course  contribute  to  the  prescattered  intensity 
I  <u'>.  One  might  include  other  direct  sources  such  as  gaseous  emission  and 
boundary  surface  radiation  plus  previously  scattered  radiation,  but  only  unscat¬ 
tered  sunlight/moonlight  is  included  in  the  present  scattering  source  function. 

The  resulting  source  function  is  found  by  using  Eq.  (2  6)  in  Eq.  (2  5)  to  obtain 
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Note  that  P^,  P^,  k^,  and  vary  with  altitude  (atmospheric  density  and  com- 

position)  and  are  generally  slowly  varying  functions  of  frequency.  Note  also  that 

the  scattering  angle  y  =  arccos  (n  •  n  ')  would  be  constant  (independent  of  the 

particular  scattering  point)  along  a  line-of-sight  in  the  absence  of  refractive 

bending.  Both  the  primary  solar  path  and  the  line-of-sight  optical  path  actually 

bend  somewhat,  so  that  y  can  be  expected  to  vary  by  as  much  as  a  few  degrees 

A+M 

along  the  line-of-sight.  The  primary  solar  transmittance  t  depends  strongly 

e ,  ps 

on  the  optical  path  length  of  the  primary  solar  path  (prior  to  scattering),  so  that 
this  factor  can  be  expected  to  vary  considerably  from  one  scattering  point  to  the 
next. 


Figure  11.  Schematic  Representation  of  the 
Single  Scattering  Geometry.  The  scattering 
point,  at  Hi'  is  shown  for  an  observer  looking 
up  from  an  altitude  Hi 
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The  monochromatic  intensity  at  the  observer  due  to  all  of  the  single  scatter¬ 
ing  sources  within  the  line-of -sight  is  obtained  by  summing  over  the  optical  path 
the  product  of  the  source  function,  and  the  transmission  function  that  gives 
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The  scattering  optical  depth  increment  kg  ds^  can  be  expressed  in  terms  of  the 
incremental  transmittance  for  both  aerosol  and  molecular  scattering  as 
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with  X  being  either  A  or  M.  The  intensity  can  therefore  be  written  as 
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which  includes  two  separate  integrals  covering  aerosol  and  molecular  scattering 
effects.  The  above  equation,  which  provides  for  a  monochromatic  calculation  at 
any  frequency  v ,  is  now  adapted  for  use  with  the  molecular  band  transmission 
model  used  in  LOWTRAN.  The  spectrally  averaged  intensity  I  is  formally  defined 
in  terms  of  a  convolution  of  the  spectral  intensity  with  a  triangular  instrument 
shape  function  g(i/)  taken  over  a  spectral  width  (ha]f  width  at  half  maximum)  of 
approximately  bv  =  10  cm  ,  that  is. 
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The  spectrally  averaged,  scattered  intensity  can  be  expressed  in  terms  of  known 
LOWTRAN  quantities  provided  that  only  the  molecular  absorption  transmittance 
is  a  rapidly  varying  function  of  frequency.  All  other  quantities  are  assumed  to 
be  constant  over  the  spectral  interval  bv.  The  result  is 
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The  quantity  pS+0p  represents  the  spectrally  averaged  transmittance  that  is 
calculated  in  LOWTRAN,  Therefore,  the  molecular  band  models  and  aerosol 
models  of  LOWTRAN  provide  a  direct  means  of  calculating  the  path  transmittance 
required  for  each  of  the  scattering  points.  In  order  to  maintain  compatibility 
with  the  spherical  shell  atmosphere  of  LOWTRAN,  the  integral  over  the  path  of 
scattering  sources  is  replaced  by  a  layer-by-layer  sum  along  the  line-of-sight. 
For  an  optical  path  traversing  N  layers  in  an  upward  or  downward  direction  this 
process  gives 
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The  quantity  A  t  is  the  change  in  molecular  or  aerosol  scattering  transmittance 
in  passing  through  layer  j,  while  <  >j  denotes  an  average  value  for  that  layer. 

Evaluating  Eq.  (33)  requires  the  equivalent  absorber  amounts  for  both  the 
line-of-sight  and  the  primary  solar  paths  associated  with  each  scattering  point 
plus  the  scattering  angle  at  each  scattering  point.  The  calculation  of  these 
amounts  and  angles  is  described  in  Appendix  C.  The  layer-by-layer  sum  of  the 
singly  scattered  intensity  is  computed  simultaneously  with  the  existing  direct 
thermal  radiance  using  the  following  expression 
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where  n  is  the  number  of  scattering  points  (layer  boundaries).  The  layer  average 
<  >j  in  Eq.  (33)  has  been  evaluated  in  Eq.  (34)  using  the  properties  of  only  the 
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two  scattering  points  that  bound  each  layer  path  segment.  In  this  scheme,  the 
observer  position  coincides  with  i  =  1  and  the  end  of  the  line -of-sight  with  i  -  n. 

If  the  optical  path  intersects  the  earth,  then  Eq.  (34)  has  an  additional  term 
representing  the  sunlight  reflected  from  the  ground.  The  ground  is  assumed  to 
be  a  diffuse  Lambertian  reflector.  The  irradiance  at  the  ground  is  proportional 
to  cos  (6),  where  0  is  the  solar  zenith  angle  at  the  ground.  The  ground  reflected 
sunlight  is  given  by  the  term 


jFUN 
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e,  ps+op, n 


A  cos  (fj) / 2 zr 


(3b) 


where  A  is  the  ground  albedo. 

The  ground  albedo  is  assumed  to  be  independent  of  frequency  and  is  read  in 
as  an  input  to  LOVVTRAN  with  a  default  of  0  (no  reflection). 

The  extra  terrestrial  solar  intensity  is  obtained  from  the  data  compiled 

by  Thekaekara.  The  intensity  is  corrected  for  variation  in  the  earth-to-sun 
distance  due  to  the  earth's  elliptical  orbit.  The  lunar  extraterrestrial  intensity 
is  obtained  by  reflecting  the  solar  intensity  off  of  the  moon's  surface  as  in 
Reference  16, 
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Here  av  is  the  wavenumber  dependent  geometric  albedo  of  the  moon  '  while 
pMOON  the  moon!s  phase  function  giving  the  relative  intensity  as  a  function  of 
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the  phase  angle  y'  of  the  moon.  Note  that  P(y’  =  0)  =  1  for  a  full  moon. 


4.3  Verification  of  the  Single  Scattering  Algorithm 

To  ensure  that  the  single  scattering  algorithm  was  correctly  implemented, 
LOWTRAN  calculations  were  compared  to  calculations  of  single  scattered  radiance 
by  an  independent,  well-developed  model.  This  model  is  a  modification  of  the 
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32 


plane -parallel,  monochromatic  multiple  scattering  code  described  in  Reference 
20.  This  model,  based  on  the  adding/ doubling  technique,  was  modified  to  compute 
only  single  scattered  radiance.  Since  LOW  TRAN  gaseous  transmission  functions 
do  not  obey  Beer's  Law,  gaseous  absorption  was  deleted  in  the  calculations  by 
setting  the  gaseous  transmittances  in  LOWTRAN  equal  to  1.0.  Statements  were 
added  to  LOWTRAN  to  calculate  and  write  total  optical  depths  and  albedos  based 
on  the  remaining  attenuation  mechanisms.  This  data  was  then  used  in  the  modified 
adding/ doubling  programs.  Comparisons  were  limited  to  cases  where  both  the 
solar  and  line-of-sight  zenith  angles  were  within  the  range  where  the  plane  paral¬ 
lel  approximation  is  valid. 

An  example  of  the  comparisons  is  show;.  ’  ,  Table  ).  The  optical  path  in  this 
case  is  from  ground  to  space  with  a  zenith  angle  of  12.  95°,  the  solar  zenith  angle 
is  also  12.  95°.  A  Henycy-Greenstein  phase  function  was  used  [see  Eqs.  (37)  and 
(38)].  The  total  optical  depth  for  the  optical  path  from  ground  to  space  is  0.  183, 
The  table  presents  the  ratio  of  the  LOWTRAN  scattered  radiance  to  the  single 
scattered  plane -parallel  calculation  for  asymmetry  factors  g  of  0  and  0.  8  and  for 
relative  azimuthal  angles  of  0°,  90°,  and  180°.  The  radiances  shown  are  the 
upward  radiances  at  100  km,  both  the  upward  and  downward  radiance  at  2  km,  and 
the  downward  radiance  at  the  ground.  In  all  cases,  the  LOW ^'R AN  scattered 
radiance  is  within  one  percent  of  the  radiance  calculated  by  the  adding/doubling 
program.  This  agreement  demonstrates  that  the  LOWTRAN  single  scattering 
algorithm  has  been  properly  implemented. 

4.4  Phase  Functions  for  Scattering  by  Atmospheric  Aerosols  and  Molecules 

The  angular  scattering  of  light  by  the  atmosphere  is  specified  by  the  phase 
function  that  gives  the  differential  probability  of  the  scattered  radiation  going  in 
a  given  direction.  The  scattering  by  the  aerosols  and  air  molecules  are  treated 
separately  using  the  appropriate  phase  function  for  each.  The  angular  distribution 
from  the  two  types  of  scattering  are  combined,  weighted  by  the  corresponding 
scattering  coefficients. 

The  phase  functions  as  used  in  the  program  are  normalized  so  the  integral 
over  all  possible  scattering  directions  (that  is,  a  sphere)  is  unity: 

If  P(y)  d  S2  =  1  (38) 

4  r 


Sharma,  S.  (1930)  An  Accurate  and  Computationally  Fast  formulation  for 
Radiative  Fields  and  Ileat  Transfer  in  General,  Plane -Parallel,  Non-Grey 
Media  With  Anisotropic  Scattering,  PhD  Thesis,  University  of  Illinois  at 
Chicago. 
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Table  1.  Ratio  of  Single  Scattered  Radiance,  LOWTRAN  6  to 
Adding/Doubling  (See  Text  for  a  Detailed  Description).  The 
figure  at  the  bottom  illustrates  the  geometry 
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with  this  normalization,  P(y)AS2  is  the  fraction  of  the  scattered  radiation  that  is 
scattered  into  a  solid  angle  AS2  about  an  angle  y  relative  to  the  direction  of  the 
incident  light. 

4.4.  1  AEROSOL  ANGULAR  SCATTERING  FUNCTIONS 

The  LOW  TRAN  program  offers  the  user  three  choices  on  handling  the  aerosol 
phase  functions: 

(1)  They  can  use  the  standard  phase  functions  stored  in  the  program  for  the 
various  aerosol  models; 

(2)  They  can  use  a  Henyey-Greenstein  type  phase  function,  with  a  specified 
value  for  the  asymmetry  parameter; 

(3)  They  can  input  their  own  phase  functions  for  the  different  altitude  regions. 

4.4.2  STANDARD  LOWTRAN  PHASE  FUNCTIONS 

The  standard  aerosol  phase  functions  stored  in  the  LOWTRAN  program  cor¬ 
respond  to  the  different  aerosol  models  available  within  the  LOWTRAN  program. 

It  is  therefore  recommended  that  this  option  be  chosen  whenever  the  LOWTRAN 
aerosol  models  are  used  for  solar  scattering.  These  standard  phase  functions 
were  developed  to  approximate  the  exact  phase  functions,  within  about  2  0  percent, 
for  any  of  the  various  aerosol  models  available  in  LOWTRAN  as  a  function  of 
wavelength,  between  0.2  and  40  Hm.  The  development  of  this  standard  set  of 
approximate  phase  functions  is  discussed  in  Appendix  D,  along  with  details  of 
their  implementation  in  the  LOWTRAN  program. 

The  number  of  phase  functions  in  this  set  represents  a  compromise  between 
accuracy  and  memory  requirements.  The  nominal  accuracy  of  20  percent  is 
compatible  with  the  other  uncertainties  in  using  the  aerosol  models  (such  as 
determining  the  concentration  of  the  aerosols).  If  greater  accuracy  is  desired 
in  specifying  the  phase  functions,  the  phase  functions  for  all  aerosol  models,  for 
a  number  of  wavelengths,  are  provided  as  a  supplemental  datu  file  on  the 

LOWTRAN  6  tape.  They  are  also  tabulated  and  discussed  in  a  separate  report 
1  1 

by  Shettlc  et  al.  “  1 

4.4.3  HENYEY-GREENSTEIN  PHASE  FUNCTION 

In  addition  to  the  standard  LOWTRAN  phase  functions  corresponding  to  the 
different  aerosol  models  built  into  LOWTRAN,  the  user  has  the  option  of  specify¬ 
ing  a  Henyey-Greenstein  scattering  function  be  used.  The  Henyey-Greenstein 
function  is  given  by: 

2  1.  Shettle,  E.  P.  ,  Abreu,  I,,  W.  ,  and  Moose,  R.  (1983)  Angular  Scattering 

Properties  of  the  Atmospheric  Aerosols,  AFGL-TR-83-  (to  be  published). 

22.  Henyey,  L.  G.  ,  and  Greenstein,  J.  L.  (1941)  Diffuse  radiation  in  the  galaxy, 
Astrophys.  J.  93:7  0-83. 
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■where  7  is  the  scattering  angle  and  g  is  the  asymmetry  parameter, 

g  =  J~  f  COS  7  P(7)dll  ,  (38) 
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with  P(y)  normalized  as  in  Eq.  (36).  The  asymmetry  parameter  gives  a  measure 
of  the  asymmetry  of  the  angular  scattering.  It  has  a  value  of  +1  for  complete 
forward  scattering,  0  for  isotropic  or  symmetric  scattering,  and  -1  for  complete 
backs  ea  tiering. 

4.4.4  USER-DEFINED  PHASE  FUNCTIONS 

The  LOWTRAN  code  allows  the  user  to  input  their  own  phase  functions  for 
the  different  altitude  regions.  These  scattering  functions  can  be  defined  at  up  to 
50  different  angles,  as  specified  by  the  user.  The  same  angles  must  be  used  for 
all  four  altitude  regions. 

When  inputting  their  own  phase  functions  the  user  should  make  sure  they  are 
normalized  as  in  Eq.  (36).  The  literature  is  not  standard  on  this  convention,  and 
other  conventions  are  used,  the  most  common  alternate  form  has  the  integral  (36) 
equal  to  4  7r- 

4.4.  5  MOLECULAR  SCATTERING  PHASE  FUNCTION 

The  angular  distribution  of  light  scattered  by  the  air  molecules  is  described 
by  the  Rayleigh  scattering  phase  function: 

P(y)  =  yyy  •  yfy  [(l  +  5)  +  (l  -  6)  cos2  7  ]  .  (39) 


where  6  is  the  depolarization  factor  that  gives  the  correction  for  the  depolarization 
effect  of  scattering  from  anisotropic  molecules.  When  6  goes  to  zero,  that  is, 
no  depolarization,  or  symmetric  molecules,  Eq,  (39)  reduces  to: 

P(l)  r  77—  [l  +  cos"  7] 

16;r 


which  is  a  commonly  used  approximation  i'or  the  Rayleigh  phase  function. 


3  6 


(40) 


A  valut'  ol  6  -  0.  UJ  i*,>  (KlisUmi  )  List'd  m  Imj,  (St?)*  ^oung“  has  givt'ii  a 
value  of  6  -  U.  02  /  !>  fur  dry  air  based  on  more  remit  measurements.  With  all 
constants  evaluated,  lap  (Hi1)  becomes 

P(->)  -  0.  0t> 05 a  ■+  0.  ib) 7 OH  co«~  .  (41) 

4,3  Sample  Calculations 

1’  igures  12  and  Id  show  some  representative  calculations  from  LOW THAN 
using  the  single  seattei  mg  option.  In  these  figures,  tin  solid  line  represents  the 
radiance  emitted  by  tin*  atmosphere,  tin  Hashed  line  the  path  scattered  radiance 
and  tilt*  dotted  lint  i'r  ground  reflected  radiance  (If  am).  The  atmospheric  pro¬ 
file  in  all  cases  is  the  1*M'2  U.S.  Standard  Atmosphere,  and  tlie  Rural  aerosol 
model  with  10-km  meteorological  range  (VIS  =  ]0  km).  The  surface  albedo  is 
0.  05,  independent  of  wavenumber  corresponding  to  an  c.missivity  of  0,  ha.  The 
solar  zenith  angle  is  4 a'  ;  results  are  shown  for  relative  path  azimuth  angles  of 
both  0  and  180  ,  (See  Appendix  t  for  a  discussion  of  the  scattering  geometry.  ) 
In  Figure  12,  the  observer  is  on  toe  ground  looking  out  to  space  with  a  zenith 
angle  of  30  .  Ine  upper  dashed  line  'vr responds  to  a  relative  azimuth  of  0°  with 
a  scattering  angle  ot  la  ,  while  the  lever  dashed  Line  is  180°  relative  azimuth 
w ith  a  scattering  angle  ot  7 5  .  1'he  ngurt  show's  licit  for  the  relative  azimuth  of 
0  ,  the  scattered  radiance  dominaies  everywhere  above  2400  cm  except  around 
the  strong  CO.,  absorption  band  centered  at  2700  cm’  \  For  the  relative  azimuth 
of  180°,  the  scattered  radiance  is  more  than  an  order  of  magnitude  less.  This 
difference  is  due  entirely  to  the  difference  m  the  .scattering  angles:  the  luie-of- 
sight  path  and  the  scattering  point -to-sun  paths  are  otherwise  identical.  The 
reason  lor  this  diifereneo  is  as  lollows:  tin.'  bulk  of  1 1  u. •  scattering  occurs  in  the 
boundary  layer  when1  the  rural  aerosol  model  aoplios.  The  phase  function  for 
the  rural  aerosol  model  inis  n  strong  torvard  peak  lor  These  wavelengths.  The 
large  difierenee  m  the  phase  luiu'tion  between  a  scattering  angle  of  la1'  and  7a11 
accounts  lor  the  large  dillerenee  in  the  scattered  radiance. 

In  1*  igure  13,  the  observer  is  at  100  km  looking  down  at  the  earth  with  a 
zenith  angle  of  laU°  (this  is  the-  reverse  of  tin-  path  in  Figure  12).  Again  the 
upper  dashed  line  corresponds  to  the  relative  azimuth  of  0°  with  a  103°  scattering 
angle  and  the  lower  one  corresponds  to  lfiii1  relative  azimuth  with  a  scattering 

23.  Kasten,  1*’.  (lhbo)  Knyieigh-cabnnnes  slreuung  in  troekenei'  luft  unter 

bei’ueksieht  iguug  neuerer  dcpolarisu  Lons-nu.  ssungen,  f hit  i k,  27  : 173- pa;. 

24.  Young,  A.T.  (if'fiu)  Revised  depolarization  corrections  for  atmospheric 

extinction,  Appl .  Opt..  lh:3427  -342ii. 
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Figure  13.  Calculated  Radiances  for  the  Following 
Conditions:  Observer  at  100  km  Looking  at  the 
Ground  With  a  Zenith  Angle  at  100  km  of  150°,  So¬ 
lar  Zenith  Angle  is  45° ,  Relative  Azimuths  of  0° 
and  180°,  Ground  Albedo  of  0.  05,  U.  S.  Standard 
Atmosphere  19G2,  Rural  Aerosol  Model  VIS  = 

10  km.  Solid  line  is  atmospheric  emission,  dashed 
lines  are  path  scattered  radiances,  dotted  line  is 
ground  reflected  radiance,  (a)  2  000  cm'*  to 
5000  cm"*  and  (b)  5000  cm"*  to  20000  cm’* 
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angle  of  165°.  In  this  case  the  ground  reflected  solar  radiance  is  greater  than 
the  path  scattered  radiance  below  5000  cm  ^  and  greater  than  the  atmospheric 
emission  above  4000  cm  ^ .  The  path  scattered  radiances  for  the  two  relative 
azimuths  are  now  quite  similar  since  the  difference  in  the  phase  function  between 
105°  and  165°  is  small.  In  the  visible  (~  17  500  cm  *)  the  ground  reflected  radi¬ 
ance  is  more  than  an  order  of  magnitude  less  than  the  path  scattered  radiance  so 
that  the  ground  is  effectively  obscured  by  the  haze  above  it.  Note,  however,  that 
the  assumed  albedo  0.  05  is  low  for  the  visible  region  of  the  spectrum. 

4.6  Recommendations  on  Usage 

The  inclusion  of  single  scattered  solar  radiance  in  LOWTRAN  6  is  a  signifi¬ 
cant  improvement  over  previous  versions  that  calculated  the  atmospheric  emis¬ 
sion  only.  The  single  scattering  approximation  is  valid  over  a  broad  range  of 
conditions  found  in  the  atmosphere.  However,  there  are  also  conditions  of  interest 
in  the  atmosphere  where  multiple  scattering  and/or  internal  sources  must  be 
included  to  accurately  calculate  the  atmospheric  radiance.  There  is  no  simple 
indicator  that  predicts  the  conditions  for  which  the  single  solar  scattering  approx¬ 
imation  is  acceptable:  rather  the  range  of  applicability  depends  upon  a  large  set 
of  parameters  including  the  atmospheric  profile,  the  optical  path,  the  solar 

geometry,  the  aerosol  phase  function,  and  the  wavenumber  region.  The  user 

14 

will  find  some  guidance  in  Section  4  of  Ridgway  et  al.  The  following  general 
comments,  in  part  drawn  from  this  source,  may  be  useful.  However,  the  user 
must  be  aware  that  they  may  not  apply  in  all  cases  and  are  indicative  only.  Also, 
these  comments  apply  only  to  the  validity  of  the  single  scattering  approximation 
and  not  the  uncertainties  ir.  the  atmospheric  data. 

1.  The  single  scattering  approximation  always  underestimates  the  scattered 
radiance  compared  to  multiple  scattering. 

2.  The  single  scattering  approximation  becomes  less  valid  with  increasing 
scattering  optical  depth  and  with  increasing  single  scattering  albedo.  For  a 
scattering  optical  depth  of  less  than  about  0.7,  the  ratio  of  multiply  scattered  to 
singly  scattered  radiances  should  be  less  than  1.  5.  For  scattering  optical  depths 
greater  than  about  2,  the  ratio  may  be  much  larger1. 

3.  For  an  observer  in  space  looking  down  at  the  ground  in  a  window  region 
where  the  total  optical  depth  is  less  than  2,  the  l'atio  of  multiply  scattered  to 
singly  scattered  radiance  is  in  most  cases  less  than  2,0  and  in  many  cases  less 
than  1.5.  And  contrary  to  intuition,  the  ratio  decreases  as  the  solar  zenith  angle 
increases  and/or  the  path  zenith  angle  decreases  (note:  a  path  zenith  angle  of 
180°  is  straight  down). 
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4.  For  an  observer  at  the  ground,  the  ratio  of  multiple  to  single  scattered 
radiance  increases  with  both  the  path  and  the  solar  zenith  angle,  and  in  general, 
single  scattering  is  a  poor  approximation  (error  greater  than  a  factor  of  2)  for 
cases  where  both  the  zenith  angles  are  greater  than  70. 

5.  Multiple  scattering  effects  are  dominant  in  clouds  and  thick  f  gs. 

6.  Single  scattering  is  a  good  approximation  for  early  twilight  c  bps,  that 
is,  where  the  sun  is  just  below  the  horizon.  For  late  twilight  cases,  multiple 
scattering  becomes  significant. 

7.  Single  scattering  is  a  good  approximation  when  looking  near  the  sun,  since 
the  scattering  is  dominated  by  the  large  forward  peak. 

8.  In  general,  the  aerosol  scattering  optical  depth  increases  with  wavenumber 
so  that  scattering  in  the  infrared  is  less  than  that  in  the  visible. 

From  a  purely  mechanical  point-of-view,  the  single  scattering  option  should 
not  be  used  along  with  either  the  cirrus  cloud  model,  or  the  rain  model  since  the 
program  does  not  compute  the  scattered  radiance  from  these  aerosols.  For  an 
optical  path  looking  directly  at  the  sun,  the  single  scattering  model  includes  only 
the  scattered  radiance,  and  not  the  transmitted  solar  radiance. 


4.7  Directly-Transmitted  Solar  Irradiance 

An  additional  option  has  been  provided  to  allow  the  user  to  compute  the 
directly  transmitted  solar  irradiance  (flux),  that  is,  the  irradiance  measured  by 
an  observer  looking  directly  at  the  sun.  This  irradiance  is  given  by: 


-  jSUN  -A+M 
v  ~  v  e 


(42) 


Note  that  all  scattered  light  is  lost  and  that  no  scattering  into  the  path  is  included. 
Instructions  for  using  this  option  are  in  Section  10.2,3.2. 

An  example  of  the  directly  transmitted  solar  irradiance  is  given  in  Figure  14, 
The  dashed  line  is  the  solar  irradiance  at  the  top  of  the  atmosphere.  The  solid 
line  is  the  transmitted  irradiance  for  a  vertical  path  from  the  ground,  for  the 
U.S.  Standard  Atmosphere  19G2  and  no  aerosol  extinction. 
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Figure  14.  Solar  Radiance  (Dashed  Dine)  and 
Directly  Transmitted  Solar  Irradiance  (Solid 
Dine)  for  a  Vertical  Path,  From  the  Ground, 
U.  S.  Standard  Atmosphere  1962,  No  Aerosol 
Extinction 


5.  NAVY  MARITIME  AEROSOL  MODEL 

S.G.  Gathman,  Naval  Research  Laboratory,  Washington,  DC 


This  chapter  provides  a  brief  description  of  the  Navy  maritime  aerosol 
model  and  its  implementation  in  LOWTRAN.  A  complete  discussion  of  the  model 
is  given  by  Gathman.  Since  this  model  includes  an  explicit  dependence  on  wind 
speed  it  is  recommended  that  it  be  used  instead  of  the  L.OWTRAN  5  maritime 

2  0 

model,  which  assumed  moderate  wind  speeds  (Shettle  &  Fenn).  The  latter 
model  is  retained  in  LOWTRAN  6  for  comparison  purposes. 


5.1  Description  of  the  Model 

The  aerosol  population  found  over  the  world's  ocean  is  significantly  different 
in  composition  and  distribution  from  that  of  a  continental  origin.  These  aerosol 
are  largely  derived  from  the  sea.  They  are  produced  by  the  evaporation  of  sea 
spray  and  from  jet  and  film  droplets.  Jet  droplets  are  ejected  into  the  air  by  the 
bursting  of  small  air  bubbles  at  the  sea  surface.  The  bursting  of  the  bubble  film 
leaves  behind  many  smaller  film  droplets  that  may  also  be  diffused  into  the  air. 

25.  Gathman,  S.G.  (1983)  Optical  properties  of  the  marine  aerosol  as  predicted 

by  the  Navy  aerosol  model.  Opt.  Eng.  22 :57 -62 . 

26,  Shettle,  E,  P.  ,  ana  Fenn,  R.W.  (197  9)  Models  of  the  Aerosols  of  the  Lower 

Atmosphere  and  the  Effects  of  Humidity  Variations  on  Their  Optical 
Properties,  AFGL-TR  -79-U2  14,  ADA085951. 


42 


These  mechanisms  are  wind  dependent  and  require  whitewater  phenomenon  in 
order  to  produce  aerosol. 

Once  the  aerosol  droplets  are  airborne,  they  undergo  additional  sorting  and 
mixing  processes.  The  marine  boundary  layer  is  usually  capped  by  a  temperature 
inversion  and,  within  this  boundary  layer,  the  smaller  marine  aerosol  together 
with  any  background  aerosol  form  a  fairly  uniform  aerosol  spatial  distribution. 

Once  introduced  into  the  atmosphere  the  lifetime  of  an  aerosol  particle  is  dependent 
on  the  size  of  the  particular  aerosol  particle.  Those  with  very  small  sizes  have 
a  very  long  residence  time  in  the  boundary  layer  if  there  are  no  washout  processes 
taking  place.  On  the  other  hand,  those  with  very  large  sizes  have  a  short  resi¬ 
dence  time  and  do  not  contribute  to  the  stationary  long-term  aerosol  population. 

The  Navy  maritime  aerosol  model  differentiates  between  these  various  types 
of  aerosol  by  postulating  that  the  marine  atmosphere  is  composed  of  three  distinct 
populations  each  of  which  is  described  by  a  lognormal  size  distribution.  The 
parameters  that  describe  the  analytical  form  of  the  size  distribution  are  then 
related  to  both  recent  meteorological  history  and  current  meteorological  observa¬ 
tions. 

The  smallest  component  of  the  model  is  a  "continental"  component.  This  is 
the  background  aerosol  and  although  it  apparently  has  little  to  do  with  current 
wind  parameters,  it  is  dependent  on  the  elapsed  time  required  for  the  air  mass 
to  traverse  the  sea  from  the  continent  to  the  point  of  observation.  Quantification 
of  this  component  in  terms  of  routine  meteorological  measurements  is  difficult, 
but  for  convenience  an  integer  from  1  to  10  is  used  to  specify  the  ICSTL  param¬ 
eter,  which  gives  a  qualitative  indication  of  the  continental  contribution:  a 
value  of  1  representing  relatively  pure  maritime  aerosol,  and  a  value  of 
ICSTL  =  10  meaning  a  significant  continental  component. 

The  second  component,  the  "stationary"  component  of  the  maritime  aerosol, 
is  the  part  of  the  maritime  aerosol  that  depends  on  the  current  and  past  history 
of  the  wind  and  represents  that  portion  of  the  spectra  that  are  produced  by  the 
high  wind  and  whitewater  phenomenon  but  do  not  fall  out  rapidly.  The  amplitude 
of  this  component  is  related  to  the  average  wind  speed  over  the  past  24  h,  and  is 
specified  by  the  WHH  parameter. 

The  third  or  "fresh"  component  of  the  Navy  aerosol  model  is  a  lognormal 
population  of  aerosol  that  is  related  to  the  current  wind  speed  (specified  in  the 
program  by  the  WSS  parameter).  The  amplitude  of  this  component  is  a  function 
of  the  current  wind  speed  and  reflects  the  current  action  of  the  production  of  drops 
produced  by  whitewater  as  a  result  of  wind  wave  actions. 

The  amplitudes  of  both  the  second  and  third  components  of  the  aerosol  popula¬ 
tion  reflect  the  necessity  of  wind  speed  being  above  a  certain  minimum  value 
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before  whitewater  phenomena  are  observed  and  thus,  marine  aerosol  produced. 
This  minimum  value  is  2.2  m/s. 

The  model  is  also  responsive  to  the  current  relative  humidity.  It  is  well 
known  that  particles  composed  of  sea  salt  are  hygroscopic  and  change  their  com¬ 
position  and  size  as  a  function  of  the  relative  humidity.  The  model  uses  the 

ii  ii  27 

'swelling  factor"  proposed  by  Fitzgerald,  '  which  adjusts  the  mode  radii  of  the 
three  components  of  the  model,  but  does  not  alter  the  total  number  of  particles 
that  are  airborne.  The  model  also  adjusts  the  complex  index  of  refraction  of  the 

2  Q 

aerosol  based  on  the  volume  weighted  method  of  Hanel,  using  the  refractive 

2  9 

index  of  soluble  aerosol  (V olz)  for  the  dry  component  and  that  of  pure  water 

30 

from  Hale  and  Querry. 


5.2  Use  of  the  Navy  Maritime  Model 

As  discussed  in  the  preceding  section,  this  model  requires  three  parameters 
to  be  specified  in  addition  to  those  used  by  the  other  aerosol  models  (that  is 
visibility  and  relative  humidity).  These  additional  parameters  are:  (l)  ICSTL, 
which  indicates  the  degree  of  continental  influence,  (2)  WHH,  the  average  wind 
speed  over  the  past  day,  and  (3)  WSS,  the  current  wind  speed.  The  LOWTRAN 
program  will  use  default  values  for  any  unspecified  parameter. 

Three  methods  can  be  used  to  estimate  ICSTL,  The  first  is  by  plotting  the 
air  mass  trajectory  and  determining  the  elapsed  time,  t,  since  the  air  parcel 
left  land.  This  time  is  related  to  ICSTL  by  the  following  empirical  equation: 

ICSTL  =  INT(9  exp  l-t(days)/4]  +  l)  ,  (43) 


where  INT(x)  truncates  to  the  nearest  integer  less  than  x. 

Secondly,  if  measurements  of  the  current  radon  222  concentrations  in  the 
atmosphere  are  available  (Larson  and  Bressan),31  then  the  air  mass  parameter 
can  be  estimated  by  the  formula: 


27.  Fitzgerald,  J.W.  (1978)  On  the  Growth  of  Aerosol  Particles  With  Relative 

Humidity,  NRL  Memo  Rpt.  3847. 

28.  Hanel,  G.  ( 1 9 7  1 )  New  results  concerning  the  dependence  of  visibility  on 

relative  humidity  and  their  significance  in  a  model  for  visibility  forecasts, 
Contrib.  Atmos.  Phys.  44:137-167. 

29.  Volz,  F.E.  (1972)  Infrared  refractive  index  of  atmospheric  aerosol  sub¬ 

stances,  Appl,  Opt.  11:755-759. 

30.  Hale,  G.M.,  and  Querry,  M.R.  (1973)  Optical  constants  of  water  in  the 

200  nm  to  200  micrometer  wavelength  region,  Appl.  Opt.  12:555-563. 

31.  Larson,  R.E.,  and  Bressan,  D.J.  (1980)  Air  mass  characteristics  over 

coastal  areas  as  determined  by  radon  measurements.  Preprint  of  Second 
Conference  on  Coastal  Meteorology,  30  January  -  1  February  1980, 

Los  Angeles,  Calif .  ;  published  by  A.  M.  S.  ,  Boston,  Mass. 


ICSTL  =  INT(Rn/4)+l  , 


(44) 
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where  Rn  is  the  concentration  of  radon  222  in  pCi/m  .  This  relationship  can  be 
used  because  radon  222  is  introduced  into  the  atmosphere  only  by  processes 
occurring  over  land.  Therefore,  since  this  radioactive  substance  has  a  half-life 
of  3.86  days,  the  concentration  of  radon  222  is  then  related  to  the  time  since  the 
air  parcel  left  the  land. 

A  third  method  for  determining  ICSTL  consists  of  subjectively  choosing  an 
integer  between  1  and  10  to  determine  the  "quality"  of  the  air  mass,  with  a  value 
of  "l"  being  for  pure  oceanic  air  and  a  value  of  "lO"  if  the  air  has  recently  been 
ashore  over  a  polluted  industrial  area.  Values  in  between  can  be  used  to  specify 
the  various  grey  areas  between  these  two  extremes. 

If  the  user  does  not  input  a  value  for  the  parameter  ICSTL,  the  LOWTRAN 
code  will  use  a  default  value  of  ICSTL  =  3. 

The  current  wind  speed,  WSS,  and  the  average  wind  speed  over  the  past  24  h, 
WHH,  should  be  input  in  units  of  m/s.  If  the  average  wind  speed  is  input  as  0  or 
is  given  a  negative  value,  a  default  value  will  be  chosen  that  depends  on  the  model 
atmosphere  being  used,  (specified  by  the  parameter  MODEL).  These  default 
wind  speeds  are  shown  in  Table  2.  These  default  wind  speeds  are  based  on  aver¬ 
age  values  for  observations  made  in  the  indicated  region,  except  for  the  user- 
defined  ca.ses  (MODEL  =  0  or  7),  which  use  a  global  mean  value.  If  the  current 
wind  speed,  WSS,  is  not  specified,  it  is  set  equal  to  the  average  wind  speed,  WHH. 

Visibility  observations  at  sea  are  usually  only  estimates  because  of  the  lack 
of  targets  at  fixed  distances  from  the  observer.  Therefore,  it  is  suggested  in 
the  use  of  this  model  that,  unless  visibility  is  measured  accurately,  the  default 


Table  2.  Default  Wind  Speeds  for  Different  Model  Atmospheres 


MODEL 

Model  Atmosphere 

Default  Wind  Speed  (m/s) 

0 

User -defined  (Horizontal  Path) 

6.  9 

1 

Tropical 

4.  1 

2 

Midlatitude  summer 

4.  1 

3 

Midlatitude  winter 

10.2  9 

4 

Subarctic  summer 

6.  G9 

5 

Subarctic  winter 

12. 35 

6 

U.  S.  Standard 

7.2 

7 

User-defined 

6.  9 

visibility  condition  be  used  because  user-specified  visibility  inputs  adjust  all  of 
the  extinction  and  absorption  coefficients  in  the  calculations  in  order  to  force  the 
calculated  extinction  at  0.  55  m  to  agree  with  the  visibility  and,  if  inaccurate,  may 
introduce  excessive  error  into  the  calculations. 

The  Navy  model  is  designed  to  operate  accurately  within  certain  limits  of 
input  parameters.  While  parameter  values  outside  of  these  limits  are  permitted 
by  the  overall  LOWTRAN  6  program,  the  accuracy  of  the  predictions  outside  of 
these  limits  is  reduced.  The  design  limits  of  the  model  parameters  are: 

50  percent  ^Relative  Humidity  £  98  percent 

0  m/s  s  WSS  s  20  m/ s 

0  m/s  <  WHH  <  20  m/s 

and 

0.  8  km  :£  VIS  <  80  km 

For  relative  humidities  or  wind  speeds  outside  these  design  limits  the  pro¬ 
gram  will  internally  reset  the  value  to  the  nearest  limit,  rather  than  try  to 
extrapolate  the  aerosol  properties. 

5.3  Sample  Results 

This  section  briefly  presents  some  sample  results  of  transmittance  calcula¬ 
tions  using  the  maritime  aerosol  model.  Figure  15  shows  the  transmittance  for 
a  10 -km  horizontal  path  at  the  surface,  for  each  of  the  standard  model  atmos¬ 
pheres  in  LOWTRAN,  These  calculations  were  all  done  using  the  default  values 
for  all  the  parameters.  Thus,  ICSTL  was  set  to  3,  the  wind  speed  depended  on 
the  model  atmosphere  (see  Table  2),  and  the  maritime  aerosol  model  calculated 
the  visibility  based  on  the  aerosol  properties  with  these  parameter  values  and  the 
relative  humidity  for  the  model  atmosphere.  The  values  of  the  parameters  for 
these  different  cases  are  summarized  in  Table  3. 

The  differences  in  the  transmittances  shown  in  Figure  15  are  due  both  to 
variations  in  the  aerosol  properties  for  different  atmospheric  conditions  and  to 
the  different  amounts  of  water  vapor  in  the  various  model  atmospheres,  The 
transmittance  in  the  3  through  5  Mm  (2000  through  3300  cm  “)  window  is  more 
sensitive  to  changes  in  the  aerosol  properties  and  in  the  8  through  1.3  Mm  (7  50 
through  12  50  cm  *)  window,  the  transmittance  is  more  sensitive  to  the  variations 
in  the  water  vapor. 


Figure  15.  Atmospheric  Transmittance  for  a  10-km  Hori¬ 
zontal  Path  at  the  surface  With  the  Navy  Maritime  Aerosol 
Model,  (a)  for  the  tropical,  midlatitude  summer,  and  mid¬ 
latitude  winter  model  atmospheres  and  (b)  for  the  subarctic 
summer,  subarctic  w inter,  andU.S.  Standard  atmospheres 
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Table  3.  Conditions  for  Sample  Runs  of  the  Navy  Maritime  Aerosol  Model 


MODEL 

Atmosphere 

Wind 
(m/  s) 

Rel  Hum 
(percent) 

Vis 

(km) 

. PH?°, 

(gm/m°) 

1 

Tropical 

4.  10 

73 

4  9.  4 

19.  0 

2 

Midlatitude  summer 

4.  10 

7  0 

52.2 

14.  0 

3 

Midlatitude  winter 

10.29 

71 

17.  8 

3.  5 

4 

Subarctic  summer 

6.  69 

72 

28.  3 

9.  1 

5 

Subarctic  winter 

12. 35 

73 

14. 2 

1.2 

6 

U.S.  Standard 

7.2 
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39.  5 

5.  9 

6.  ARMY  VERTICAL  STRUCTURE  ALGORITHM 
M.G.  Heaps,  Army  Atmospheric  Sciences  Laboratory 
Wliite  Sands  Missile  Range,  New  Mexico 


6.1  Introduction 

An  algorithm  for  modeling  the  vertical  structure  of  aerosols  has  been  added 

to  the  LOWTRAN  fi  code.  It  was  developed  initially  to  describe  the  vertical 

distribution  of  the  atmospheric  aerosols  for  conditions  of  limited  visibility  and 
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beneath  low -lying  stratus  cloud  decks.  The  formalism  has  been  extended  so 
it  can  also  represent  cases  with  no  cloud  ceiling  and  moderate  to  high  visibility.33 
The  algorithm  will  generate  the  vertical  aerosol  profile  within  the  boundary  layer 
from  input  parameters,  such  as  surface  visibility  and  the  cloud  ceiling'  height. 
This  model  is  designed  for  use  within  the  lowest  2  km  of  the  atmosphere. 


6.2  The  Vertical  Profile  Model 

In  low  visibility  situations,  due  either  to  haze  or  fog,  increasing  numbers  of 
observations  show  that  the  measured  visibility  at  the  surface  is  not  representative 
of  conditions  a  few  hundreds  of  meters,  or  even  tens  of  meters,  above  the  surface. 
Thus,  the  "slant  path  visibility"  can  be  significantly  different  from  the  'horizontal 
visibility".  In  a  sigiificant  fraction  of  the  cases  the  visibility  becomes  worse  as 
the  height  above  the  surface  increases.  These  cases  are  of  special  concern  here. 


32.  Heaps,  M.G.  (1982)  A  Vertical  Structure  Algorithm  for  Low  Visibility/Low 

Stratus  Conditions,  ASL-TR-0111,  US  Army  Atmospheric  Sciences 
Laboratory,  White  Sands  Missile  Range,  N,  Mex, 

33.  Heaps,  M.  G.  ,  and  Johnson,  R.D.  (1983)  An  Empirical  Algorithm  for  the 

Vertical  Structure  of  Atmospheric  Extinction,  ASL-TK-0112,  US  Army 
Atmospheric  Sciences  Laboratory,  White  Sands  Missile  Range,  N.  Mex. 
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Detailed  data  on  the  vertical  structure  of  fogs  and  hazes  have  been  gathered 
in  the  Federal  Republic  of  Germany  on  several  different  occasions.  Droplet 

size  distributions  in  the  0.  5-  to  47 -Mm  range  have  been  measured  from  a  balloon- 
borne  instrument,  thus  yielding  vertical  profiles.1  Extinction  coefficients  at 
desired  wavelengths  or  the  liquid  water  content  can  be  calculated  from  these 
measured  droplet  size  distributions. 

The  vertical  structure  of  these  profiles  has  been  examined  previously  by 
37 

Duncan  et  al,  who  characterized  the  vertical  structure  in  the  form 

y  =  a'  x  +  b'  ,  (45) 

where  x  =  log^g  k(z),  y  -  log^Q  k(z  +  20),  a'  and  b'  are  coefficients  that  were 
chosen  to  fit  the  data,  and  k(z)  is  the  value  of  the  extinction  coefficient  at  the 
altitude  z;  k(z  +  20)  is  then  the  value  of  this  variable  at  an  altitude  of  z  +  20  m. 
Thus,  one  can  work  stepwise  from  the  surface  up  through  the  cloud  boundary 
layer.  Figure  16  shows  the  fit  of  Eq.  (45)  to  the  data.  It  should  be  noticed  that 
there  is  a  sharp  change  in  slope  at  a  value  of  about  7.  1  km  1  for  the  extinction. 
The  physical  significance  of  this  inflection  point  is  discussed  below. 

The  point  of  intersection  of  the  two  line  segments  physically  represents  the 
changes  in  extinction  due  to  changes  in  the  state  of  particle  growth  as  one  moves 
from  a  subsaturated  environment  (lower  line  segment),  where  relative  humidities 
are  less  than  100  percent,  to  a  supersaturated  environment,  (upper  line  segment). 
Thus,  this  point  of  intersection  will  be  taken  to  represent  the  cloudbase  or  lower 
cloud  boundary. 
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Figure  16.  Relationship  Between  the  Extinction  Coef¬ 
ficient  (at  0.  55  m m)  at.  Altitudes  z  and  z  +  20  m.  The 
vertical  linesare  the  error  bars  for  the  data  (after 
Duncan  et  al)37 


Since  x  =  log10  k(z )  in  Eq.  (45)  and  y  is  really  just  x  +  Ax  over  an  altitude 

interval  Az,  the  relation  for  the  extinction  as  a  function  of  altitude  can  be 
32 

expressed  as 

kg(0.  55  pm)  =  A  exp  [B  exp  (Cz)]  ,  (46) 


where  A,  B,  and  C  are  functions  of  preselected  boundary  values,  the  initial  or 
starting  value  of  extinction,  and  the  cloud  ceiling  height.  Note  that  since  there 
are  two  straight  line  segments,  the  coefficients  A,  B,  and  C  have  different 
values,  depending  on  which  part  of  the  data  curve  in  Figure  16  is  being  followed. 

The  rate  at  which  the  extinction  changes  with  altitude  below  the  cloudbase 
actually  depends  on  the  cloud  ceiling  height  z^.  The  lower  line  segment  in  Fig¬ 
ure  16  represents  an  average  of  several  sets  of  data  and  therefore  gives  a  single, 


average  value  for  the  coefficient  C  in  Kq.  (46).  The  explicit  dependence  of  C  on 
the  cloud  ceiling  height  can  be  incorporated  by  defining  the  coefficient  C  as 


C 


ln(E/A) 

ln(D/A) 


(47) 


where  F  is  the  value  of  the  extinction  at  the  eloudbase,  D  is  the  observed  value  of 
extinction  at  the  surface,  and  A  is  the  same  coefficient  used  in  £q.  (46).  In  this 
case  A  is  the  lower  limit  to  the  extinction  in  the  hazy/foggy  region  below  the 
cloud.  Figux'e  17  shows  the  visible  extinction  coefficient  plotted  as  a  function  of 
altitude  for  the  same  initial  surface  value,  but  several  different  cloud  ceiling- 
heights.  The  solid  line  represents  the  (average)  values  from  the  line  segments 
in  Figure  16.  The  dashed  vertical  line  represents  the  value  of  the  extinction  for 
the  eloudbase  given  by  the  intersection  of  the  line  segments  in  Figure  10.  The 
solid  line  to  the  right  of  the  dashed  vertical  line  represents  the  extinction  profile 
inside  the  cloud  and  can  be  appended  to  any  one  of  the  vertical  profiles  to  the  left. 


VISIBILITY  (km) 

10.  10  0.1 


Figure  17.  The  Vertical  Px-ofilc  of  the  0.  5 5 -pin  Extinc¬ 
tion  Coefficient  foi  Various  Cloud  Ceiling  Heights.  The 
solid  line  shows  the  average  profiles  from  Figure  1(1, 
the  dashed  vertical  line  represents  the  value  at  the  cloud 
boundary 
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Initially  the  algorithm  for  the  vertical  structure  of  hazes,  fogs,  and  clouds 
represented  by  Eq.  {4(0  was  developed  for  low  visibility/low  stratus  conditions 
and  is  based  on  inputs  of  the  surface  meteorological  range  (extinction  coefficient) 
and  the  cloud  ceiling  height.  This  algorithm  has  now  been  extended  to  eases 
where  there  may  be  no  cloud  ceiling  and  where  the  extinction  coefficient  decreases 
with  increasing  altitude. 

6.3  Applicability  of  I  lie  Vertical  Structure  Algorithm 

Three  initial  visibility  conditions  are  considered.  The  first  condition  is  for 
strstus  clouds  and  thick  fogs  (which  in  this  instance  may  be  treated  as  a  cloud  at 
the  ground);  the  second  condition  is  for  hazes  and  fog;  and  the  third  condition  is 
for  the  clear  to  hazy  atmosphere.  The  vertical  structure  of  visibility  can  be 
represented  by  foui  different  types  of  curves  as  illustrated  in  Figure  115.  Curves 
1  and  2  represent  the  cases  where  the  extinction  coefficient  increases  (that  is, 
visibility  degrades)  with  increasing  altitude;  these  cases  are  representative  of 
the  vertical  structure  of  extinction  for  thick  fogs  or  for  low  visibility/ low  stratus 
conditions.  Curves  3  and  4  represent  eases  where  the  extinction  coefficient 
decreases  (that  is,  visibility  improves)  with  increasing  altitude.  Each  of  thesi 
eases  will  now  be  briefly  outlined. 

Case  1:  This  curve  is  to  be  used  for  dense  fogs  at  ground  level  or  when  one 
is  at  the  cloudbase  or  in  the  cloud.  Physically  this  curve  represents  the  increase 
in  liquid  water  content  (LWC),  auci  consequently  the  increase  in  extinction  coef¬ 
ficient  and  decrease  in  visibility,  of  a  saturated  parcel  of  air  rising  at  the  wet 
adiabatic  lapse  rate.  This  curve  should  be  used  only  when  the  initial  extinction 
coefficient  (or  meteorological  range)  is  in  the  thick  fog/ cloud  region  shown 
between  the  two  dashed  lines  representing  boundary  values  on  the  right-hand 
side  of  Figure  18. 

Case  2  :  This  curve  is  to  be  used  for  low  visibility  conditions  beneath  the 
clouds  due  to  haze  or  fog  when  there  is  a  low  cloud  ceiling  present. 

Case  3:  This  curve  is  to  be  used  when  there  is  a  shallow  radiation  fog 
present  or  when  a  haze  layer  is  capped  by  a  distinct  (low-lying)  temperature 
inversion.  A  cloud  ceiling  is  not  present. 

Case  4:  This  profile  is  to  be  used  for  eases  where  there  is  reasonable 
vertical  homogeneity  for  visibility  in  a  clear  to  slightly  hazy  atmosphere  that 
may  have  a  shallow  haze  layer  near  the  surface.  A  cloud  ceiling  is  not  present. 


VISIBILITY  (km) 
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Figure  18.  Four  Different  Cases  Represented  by  the  Vertical 
Structure  Algorithm 


Profiles  of  the  0.  55-^m  extinction  coefficient  are  shown  in  Figure  18  for 
these  different  cases.  Two  examples  of  representative  profiles  are  shown  for 
case  1.  The  first,  example  is  for  a  thick  fog  at  the  surface,  which  is  represented 
by  an  extinction  coefficient  profile  that  increases  with  height.  When  the  depth  of 
the  fog  is  not  known  (which  is  usually  the  case  because  the  sky  is  obscured),  a 
default  depth  of  200  m  is  recommended.  The  second  example  is  for  a  low-lying 
stratus  cloud;  for  illustration  the  cloud  ceiling  height  is  taken  to  be  2  00  m.  This 
profile  should  only  be  used  from  the  eloudbase  to  the  cloud  top.  Again,  cloud 
thickness  is  usually  not  a  measured  quantity,  and  a  default  value  of  2  00  m  is 
recommended.  The  two  examples  shown  here  ai'e  actually  the  same  profile 
one  starting  at  the  surface  for  the  thick  fog  and  the  other  starting  at  the  cloud 
base  of  a  low-lying  stratus  cloud.  Within  region  1,  thick  fog/ cloud,  only  pro¬ 
files  of  the  case  1  type  should  be  used,  For  a  dense,  shallow  radiation  fog,  use 
a  profile  for  case  3  as  described  below. 


A  representative  profile  for  the  structure  beneath  a  stratus  cloud  is  shown 
for  case  2.  In  this  instance  the  visibility  conditions  at  the  surface  are  represent¬ 
ative  of  region  2,  haze/fog,  and  the  cloud  ceiling  height  is  200  m.  The  slope  and 
shape  of  the  vertical  structure  profile  beneath  the  cloud  deck  are  a  function  of  the 
initial  value  of  the  (surface)  visibility  and  the  cloud  ceiling  height.  For  haze/fog 
conditions  when  a  cloud  ceiling  height  less  than  2  km  is  present,  a  profile  of  the 
case  2  type  should  always  be  used. 

Often  a  low-lying  cloud  cover  is  present  when  the  surface  visibility  is  clear 
to  only  slightly  hazy.  In  this  instance  a  vertical  structure  profile  similar  to 
case  2  is  appropriate.  This  profile  is  denoted  as  case  2'  and  is  shown  in  Figure 
18  as  an  alternate  profile  for  the  instance  of  a  2  00-m  cloud  ceiling  height.  The 
only  difference  between  case  2  and  case  2'  is  the  manner  of  choosing  the  value  of 
the  coefficient  A,  which  in  turn  influences  the  shape  of  the  vertical  profile. 

A  shallow  radiation  fog  or  a  haze  layer  bounded  by  a  temperature  inversion 
can  be  represented  by  a  vertical  structure  profile  as  shown  in  Figure  18  for 
case  3.  The  boundary  layer  heights  for  such  occurrences  are  often  difficult  to 
estimate.  Temperature  inversion  heights  can  be  obtained  from  acoustic  sounders 
or  radiosonde  observations;  often  visual  sightings  can  be  used  to  estimate  depths 
of  shallow  fogs  or  haze  layers.  A  nominal  boundary  layer  height  of  2  00  m  has 
been  selected  for  illustrative  purposes.  For  radiation  fogs  where  the  depth  is 
not  known,  a  default  value  of  200  m  is  selected1';  for  inversion  layers  where  the 
height  of  the  inversion  or  boundary  layer  is  not  known,  a  default  value  of  2  km  is 
selected. 

Case  4  is  represented  by  a  profile  for  the  condition  where  the  vertical 
structure  is  essentially  constant  with  altitude,  with  the  exception  of  the  lowest 
hundred  meters  of  the  boundary  layer.  An  appropriate  default  value  is  the  nomi¬ 
nal  background  value  for  the  0.  55-qm  extinction  coefficient  for  the  fair  weather 
case.  Numerous  observations  have  shown  that  the  extinction  coefficient  is 
essentially  constant  within  the  planetary  layer  for  well  mixed  conditions.  Setting 
the  coefficient  C  equal  to  zero  in  Fq.  (40)  will  cause  the  algorithm  to  default  to 
the  observed  surface  value  while  providing  a  constant  vertical  profile. 

Table  4  gives  the  tabular  values  of  the  0.55-gm  extinction  coefficients  that 
are  to  be  used  as  boundary  values  for  the  different  cases  in  their  respective 
regions  of  applicability. 


Note  added  in  Proof.  A  more  realistic  value  fur  radiation  fogs  is  about  50  rn . 
To  override  the  default  value  in  the  program,  read  in  the  depth  of  the  radiation 
fog  for  this  case  with  ZINVSA  =  0.  05  km  (see  Section  10.2.2.  1,  page  00). 
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Table  4.  Summary  of  the  Conditions  and  Parameter  Values  for  Different  VSA  Cases 

k  (0.  55  ^m!  =  A  exp  [B  exp  (Cz)] 
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boundary  layer 
height  (m) 

A  nominal  background  value  for  the  lower  troposphere. 

When  the  fog  or  cloud  depth  is  not  known  (usually  the  case),  the  algorithm  should  be  used  only  up  to  200  m  above 
the  surface  or  cloud  boundary  (a  default  value  for  the  stratus  cloud  thickness). 


6.4  Implementation  of  the  Vertical  Structure  Algorithm  in  L0WTRAN6 

The  operation  of  the  Vertical  Structure  Algorithm  (VSA)  is  controlled  by 
three  parameters,  in  addition  to  the  Meteorological  Range  at  the  surface  (VIS) 
and  type  of  aerosol  (IHAZE)  for  the  boundary  layer.  These  three  additional 
parameters  are:  the  cloud  ceiling  height  (altitude  of  the  cloud  base),  the  thickness 
of  the  cloud  or  fog,  and  the  height  of  the  inversion  or  boundary  layer,  ZCVSA, 
ZTVSA,  and  ZINVSA  respectively.  The  type  of  aerosol  vertical  profile  generated 
depends  on  the  values  input  for  these  parameters.  The  different  cases  or  profile 
types  selected  are  summarized  in  Table  5.  Note  that  the  value  of  ZINVSA  will  be 
ignored  unless  ZCVSA  <  0. 

The  VSA  defines  the  aerosol  extinction  at  nine  heights,  from  the  ground  to 
the  top  of  the  cloud  (ZCVSA  +  ZTVSA)  or  the  top  of  the  boundary  layer  (ZINVSA). 
Ten  meters  above  this  cloud-top  or  boundary-layer  level,  the  aerosol  profile 
reverts  to  the  standard  LOWTRAN  aerosol  vertical  distribution  (or  the  user- 
supplied  profile  for  a  MODEL  =  7  case).  For  these  nine  heights  the  air  pressure, 
temperature,  and  the  ozone  concentration  are  found  by  interpolation  from  the 
model  atmosphere  indicated  by  the  parameter  MODEL  (see  the  user  instructions. 
Section  10). 

The  relative  humidity  for  the  LOWTRAN  model  atmospheres  does  not  con¬ 
sider  the  presence  of  clouds  (that  is,  all  the  model  atmospheres  have  RH  £80  per¬ 
cent  at  all  altitudes).  To  correct  for  this,  the  VSA  estimates  the  relative  humidity 
as  a  function  of  the  visible  (X  -  0.  5  5  pm)  extinction  for  the  nine  levels: 

•  lntke(z)]  ,  kp  <  7.064  km  1 

,  k  >7.064  km"1  (48) 

e 

If  the  user  inputs  their  own  relative  humidity  profile  (MODEL  =  7),  that  will  be 
used  instead  of  Eq.  (48). 


86.407  +  6. 953 


RH  (z) 


100% 


50 


Table  5.  Data  Inputs  and  Default  Values  for  the  Different  VSA  Cases 


7.  CIRRUS  CLOUD  MODEL 

F.F.  Hall,  Jr.,  M.J.  Post,  R.A.  Richter,  G.M.  Lerfald,  and  V.E.  Derr 
National  Oceanic  and  Atmospheric  Administration 
Boulder,  Colorado 

7.1  Introduction 

This  section  of  the  report  describes  the  development  of  a  cirrus  cloud  trans¬ 
mittance  model  and  computer  subroutine  for  use  with  LOWTRAN.  Previous 

38-44 

measurements  and  models  of  cirrus  optical  transmittance  have  indicated 

large  variations  in  the  attenuation  coefficient  for  cirrus  in  the  near  ultraviolet, 
visible,  and  infrared  spectral  regions.  In  the  model  developed  here,  the  cirrus 
attenuation  coefficient  is  shown  to  be  proportional  to  the  cloud  thickness  and 
independent  of  wavelength  from  the  ultraviolet,  0.317 -pm,  to  the  infrared,  10-pm 
atmospheric  window.  The  model  is  discussed,  worldwide  cirrus  statistics  are 
presented,  and  thickness  statistics  are  tabulated. 


7.2  Worldwide  Cirrus  Climatologies 

A  comprehensive  study  of  worldwide  cloud  statistics  has  been  done  by  Chang 
45 

and  Willand.  Their  study  used  reports  of  the  surface -based  observer  in  record- 

4  g 

ing  the  presence  or  absence  of  high  clouds.  However,  as  pointed  out  by  Stone, 
"where  low  and  middle  clouds  are  frequent,  the  ground  observations  of  cirrus 
probably  record  but  50%  of  the  true  frequency".  Therefore,  the  ground-based 

38.  Fritz,  S.  ,  and  Rao,  P.  K.  (1967)  On  the  infrared  transmission  through 

cirrus  clouds  and  the  estimation  of  relative  humidity  from  satellites, 

J.  Appl.  Meteorol.  6:1088-1096. 
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observer  climatologies  of  Reference  45  were  multiplied  by  1.25  to  give  a  more 
representative  probability  of  cirrus  present.  In  Table  6  the  recommended  cirrus 
climatologies  for  the  tropics  (inter -tropical  convergence  zone),  temperate,  and 
subarctic  regions  are  presented. 


Table  6.  A  Recommended  Cirrus  Climatology  Based  on  Surface-based 
Observations,  Compensated  for  Obscuration  by  Lower  Clouds 


Latitude  zone 

Tropics  (ITCZ) 
10°S  -  2  0°N 

Temperate  or 
Midlatitude 

Subarctic  or 
Polar  latitudes 

Summer 

Winter 

Summer 

W  inter 

Fraction  of  time  cirrus 
clouds  are  present 

0.  80 

0.  40 

0.  50 

0.  45 

0.40 

Even  though  ground  observations  may  show  an  apparent  diurnal  variation  in 
the  amount  of  cirrus  coverage,  this  is  most  probably  caused  by  the  difficulty  of 
observing  cirrus  during  the  nighttime  hours.  Therefore,  we  recommend  following 
the  conclusion  of  Stone48  that,  "there  is  probably  no  real  diurnal  variation  in 
cirrus  except  where  air  mass  thunderstorms  or  orographic  effects  are  important". 
Also,  the  long  lifetimes  of  cirrus,  attributable  to  the  slow  evaporation  of  the 
typically  large  ice  crystals,  probably  result  in  little  diurnal  variation. 

7.3  Cirrus  Height  and  Thickness  Statistics 

The  distinguishing  characteristic  of  cirrus  and  cirrostratus  clouds,  according 

47 

to  the  International  Cloud  Atlas,  is  that  these  clouds  are  "composed  almost 
exclusively  of  ice  crystals".  Because  of  the  usual  scarcity  of  freezing  nuclei 
active  above  -2  0°C  in  the  atmosphere,  cirrus  clouds  are  not  usually  found  at 
temperatures  higher  than  this.  Except  in  polar  regions  and  in  wintertime  tem¬ 
perate  zones  near  the  Arctic  or  Antarctic,  these  low  temperatures  generally 
occur  only  in  the  middle  or  upper  troposphere.  Thus,  cirrus  are  usually  classed 
as  high  clouds,  although  by  definition  ice  fogs  (in  German,  "Diamantstaub")  can 
also  be  called  cirrus. 

In  the  Compendium  on  Cirrus  Forecasting  by  Stone,  48  cirrus  thickness 
statistics  measured  by  the  British  Meteorological  Research  flights  and  by  the 
U.  S.  Air  Force  Project  WIBACK  are  given.  In  addition,  Weickmann48  and  Alt49 
present  cirrus  thickness  measurements.  When  these  various  statistics  are 

References  47  to  49  will  not  be  listed  here.  See  References,  page  151. 
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plotted  on  a  graph  of  thickness  vs  the  percentage  of  time  that  a  given  thickness 
occurs,  a  histogram  as  shown  in  Figure  19  results.  This  highly-skewed  histogram 
shows  a  median  cirrus  thickness  at  1.  0  km.  When  the  logarithms  of  the  thick¬ 
nesses  shown  in  Figure  19  are  plotted  vs  probability  of  occurrence,  a  straight 
line  results,  except  for  the  most  thin  and  thick  extremum  points.  Cirrus  thick¬ 
nesses  thus  seem  to  occur  in  a  truncated  lognormal  distribution.  This  is  not 
unexpected  since  Lopez"  '  has  shown  that  other  (convective)  cloud  properties 
are  so  distributed  Most  of  the  moisture  at  cirrus  levels  must  reach  these 
heights  through  convective  clouds.  Cirrus  with  the  thickness  statistics  shown 
in  Figure  19  can  occur  throughout  the  middle  and  upper  troposphere  and  occasion¬ 
ally  into  the  lower  stratosphere  where  the  moisture  may  have  reached  such  levels 
through  breaks  in  the  tropopause.  In  Table  7  the  height  ranges  where  cirrus  may 
be  expected  are  presented  based  on  the  average  temperature  profiles  for  the 
latitudes  shown. 


0  I  2  3  4  5 

CIRRUS  THICKNESS  L  (km) 


Figure  19.  Histogram  of  Cirrus  Thickness,  Showing  the  Percentages 
of  Time  that  Different  Thicknesses  (in  500-m  Increments)  are  Found 
When  Cirrus  are  Present.  Median  thickness  is  1.0  km 


50.  Lopez,  R.L.  (1977)  The  lognormal  distribution  and  cumulus  cloud  popula¬ 

tions,  Mon.  Wea,  Rev.  10J5 :8G5-872 . 

51.  Lopez,  R.E.  (1978)  The  determination  of  convective  shower  populations 

from  radar  data,  Preprint,  18th  Conf,  Radar  Mctcorol.  ,  Amer,  Meteorol. 
Soc.  ,  Boston,  Mass.  ,  155-158. 


Table  7.  Cirrus  Height  Ranges 


,1 


Tropics 

7.  50 

-  16.  50  km 

Temperate,  winter 

4.  50 

-  14.00  km 

Temperate,  summer 

7.30 

-  13.  50  km 

Polar,  winter 

2.  50 

-  10.  00  km 

Polar,  summer 

4.  50 

-  9.  50  km 

7.4  Cirrus  Transmittance  Models 

A  number  of  investigators  have  constructed  transmittance  models  for  water 
droplet  clouds^2  and  most  of  these  calculations  are  based  on  extinction  by 
Mie  theory  scattering  and  absorption  in  the  spherical  cloud  droplets  for  known 
cloud  droplet  size  distribution.  Only  Plass  and  Kattawar  considered  ice  clouds 
in  their  models.  Their  ice  crystals  were  assumed  to  be  spherical  and  the  num¬ 
bers  and  sizes  of  ice  spheres  used  are  not  too  realistic  based  on  measurements 
that  have  since  become  available.  Knowledge  of  actual  cirrus  ice-crystal  size 
distributions  has  been  greatly  improved  through  a  series  of  measurements  made 
by  the  Air  Force  Geophysics  Laboratory  in  the  late  1970's.  Derr43  has 

used  these  recent  cirrus  ice  crystal  sizes  and  number  density  information  to 
model  the  transmission  through  cirrus,  assuming  ice  spheres  of  equivalent 


52.  Deirmendjian,  D.  (1964)  Scattering  and  polarization  properties  of  water 

clouds  and  hazes  in  the  visible  and  infrared,  Appl.  Opt.  3:187-196. 

53.  Plass,  G.  N.  ,  and  Kattawar,  G.  W.  (197 1)  Radiative  transfer  in  water  and  ice 

clouds  in  the  visible  and  infrared  region,  Appl.  Opt.  10:738-748. 

54.  McClatchey,  R.A.,  Fenn,  R.W.,  Selby,  J.E.A.,  Volz,  F.E.,  and 
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56.  Manton,  M.  J.  (1980)  Computations  of  the  effect  of  cloud  properties  on  solar 

radiation,  J.  Recherches  Atmospheriques,  14:1-16. 
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volume  to  the  nonspherical  particles.  If  the  attenuation  coefficient  calculated  by 
Derr  is  plotted  vs  the  cirrus  ice  crystal  number  density,  (particles  per  cubic 
meter),  the  plot  of  Figure  20  results.  It  is  clear  that  for  the  realistic  particle 
number  densities  and  sizes  for  the  clouds  so  modeled  there  appears  to  be  no 
relationship  between  number  density  and  attenuation  coefficient.  If,  on  the  other 
hand,  the  attenuation  coefficient  is  plotted  vs  the  equivalent  LWC  for  the  cirrus, 
we  obtain  the  rather  linear  plot  of  Figure  21.  This  same  linear  relationship 
between  extinction  and  LWC  has  been  found  to  hold  in  the  infrared  for  liquid 
droplet  clouds.  ’  A  number  of  the  AFGL  cirrus  particle  measurements 
give  both  cirrus  LWC  and  cloud  thickness.  When  a  plot  is  made  of  these 

two  variables,  the  graph  shown  in  Figure  22  results.  It  appears  that  LWC  is 
directly  proportional  to  cloud  thickness  and  since  the  attenuation  coefficient  is 
also  directly  proportional  to  LWC,  then  the  attenuation  coefficient  must  be  pro¬ 
portional  to  cloud  thickness  also. 


Figure  2  0.  Calculated  Cirrus  Attenuation  Coefficients  for 
Solar  Radiation  Plotted  vs  Ice  Crystal  Number  Density. 
The  numbers  on  the  data  points  refer  to  model  numbers 
used  by  Derr 


References  61  to  63  will  not  be  listed  here.  Sec  References,  page  151, 
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CIRRUS  LIQUID  WATER  CONTENT  OR  DENSITY  (g  m"s) 

Figure  21.  Calculated  Cirrus  Attenuation  Coefficients  for 
Solar  Radiation  Plotted  vs  Cloud  Liquid  Water  Content 
(LWC).  Again,  the  numbers  are  the  cirrus  model  num¬ 
bers  used  by  Derr 


M  l0° 


O  O  o 


AFGL  LWC  vs 
CLOUD  THICKNESS 


CIRRUS  LIQUID  WATER  CONTENT  (g  nrf5) 


Figure  22.  Cirrus  Cloud  Thickness  vs  Liquid  Water  Content 
(LWC)  From  AFGL  Measurements 


7.5  Grrus  Transmittance  and  Thickness  Data 

As  part  of  a  Department  of  Energy  study,  NOAA  investigated  the  transmit¬ 
tance  of  all  cloud  types  using  a  nine -channel,  filtered  photometer/radiometer  and 

64 

a  pulsed  ruby  lidar  to  measure  cloud  height  and  thickness.  Data  were  obtained 
at  Erie  and  Grover,  Colorado;  White  Sands  Missile  Range,  New  Mexico;  Colstrip, 
Montana,  and  later  at  Point  Mugu,  California.  Time  lapse,  all-sky  camera 
coverage  allowed  choosing  cases  when  cirrus  clouds  alone  contributed  to  solar 
radiation  attenuation.  It  was  necessary  to  use  the  camera  because  the  lidar 
pointed  vertically  to  monitor  cloud  height  and  thickness,  but  the  photometer/ 
radiometer  tracked  the  sun.  Forty  independent  cases  were  selected  for  analysis 
and  these  cases  were  chosen  for  clouds  with  trajectories  that  passed  through  the 
zenith,  then  eclipsed  the  sun,  or  for  clouds  that  eclipsed  the  sun  and  then  passed 
through  the  zenith.  Of  course  it  was  impossible  to  determine  that  the  exact  same 
portion  of  the  adverting  cirrus  clouds  affected  the  photometer  and  lidar  with 
perfect  time-lagged  spatial  correlation.  Cirrus  are  notorious  for  large  varia¬ 
tions  of  transparency  even  when  extensive  cirrostratus  sheets  are  encountered. 

This  variability  led  Valovcin  ’  to  despair  of  ever  obtaining  a  cirrus  transmittance 
model.  Nonetheless,  when  the  attenuation  coefficient  calculated  for  the  measured 
cirrus  thickness  (using  the  photometer  channel  at  0.  97  5  Mm)  was  plotted  vs  the 
cloud  thickness  for  the  forty  cases,  we  obtained  the  plot  of  Figure  23.  For  cases 
where  cloud  thickness  was  not  rapidly  varying  with  time,  single  circles  are 
plotted.  For  cases  where  scattered  cirrus,  that  seemed  to  be  fibrous  and  cover¬ 
ing  the  entire  sky  but  could  not  be  sampled  through  the  same  portion  of  the  cloud 
with  the  lidar  and  photometers,  X's  are  plotted.  For  cases  when  multiple  layers 
were  measured  with  the  lidar,  the  number  of  significant  layers  so  measured  is 
indicated  by  the  number  above  the  line  connecting  the  two  extremes  in  thickness 
of  the  clouds  for  the  several  minutes  when  that  same  portion  of  the  cloud  eclipsed 
the  photometer.  Although  there  is  significant  scatter  to  the  data,  the  stx'aight  line 
shown  on  Figure  2  3  represents  the  best  experimental  fit  to  the  forty  case  data  set. 

For  most  of  the  thinner  clouds  with  total  optical  thickness  (attenuation  coef¬ 
ficient  multiplied  times  the  cloud  thickness),  less  than  0.  5,  very  little  spectral 
dependence  could  be  found  over  the  entire  range  examined,  0.  317  Mm  to  the 
8  through  12  Mm  window  region.  For  the  thicker  clouds,  the  optical  thickness 
seemed  to  increase  at  the  longer  wavelengths,  but  this  could  well  have  been  an 

64.  Lerfald,  G,  M.  ,  Derr,  V.E.,  Abshire,  N.  L.  ,  Cupp,  R.E.,  and  Ericson,  H.L. 

(1980)  Optical  properties  of  clouds  and  aerosols  derived  from  ground-based 
remote  sensing  methods,  Role  of  Electro -Optics  in  Photovoltaic  Energy 
Conversion,  SP1E  vol.  248,  166-17  1. 

65.  Valovcin,  F.R.  (1968)  Infrared  measurements  of  jet-stream  cirrus,  J.  Appl. 

Meteorul.  7:817-826. 
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Figure  23.  NOAA  Measurements  of  Cirrus  Transmittance. 

Solar  irradiance  was  measured  using  a  nine-channel  photom¬ 
eter/radiometer  and  cloud  thickness  with  a  pulsed  ruby  lidar 

artifact  of  the  larger  field-of-view  of  the  short  wavelength  photometer.  The 
smaller  field-of-view  of  the  infrared  radiometer  did  not  sample  as  much  forward 
scattered  radiation  and,  of  course,  the  smaller  effective  size  parameter  at  t.he 
infrared  wavelengths  would  lead  to  a  wider  scattering  diagram.  Our  best  estimate 
is  that  the  attenuation  coefficient  for  cirrus  is  independent  of  the  wavelength  over 
the  range  measured  during  our  study,  or  0.  317  to  10  Mm. 

7.6  An  Analytical  Model  of  Cirrus  Transmittance 

The  best  fit  to  the  experimental  data  shown  on  Figure  2  3  is  expressed  by  the 
relationship,  k  =  0.  14  L,  where  the  attenuation  coefficient,  k,  is  in  the  units  per 
kilometer  and  the  cloud  thickness,  L,  is  measured  in  kilometers.  When  the 
experimental  data  is  plotted  on  the  same  graph  as  the  AFGL  LWC  model  based 
on  Derr's  calculations,  the  graph  of  Figure  24  results.  It  is  seen  that  the  LWC 
model  leads  to  a  slightly  lower  attenuation  coefficient  for  most  cirrus  clouds. 
Because  of  the  uncertainties  in  the  model,  especially  the  assumption  of  the 
equivalent  size  spheres,  it  is  felt  that  the  experimental  data  should  be  used.  It 
is  encouraging,  however,  to  see  that  the  model  and  experimental  data  agree 
rather  closely  in  spite  of  the  many  assumptions  in  the  model  and  the  difficulty  in 
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ATTENUATION  COEFFICIENT,  (km‘l) 


Figure  24.  A  Comparison  of  the  NOAA  Experimental 
and  AFGE  Eiquid  Water  Content  (LWC)  Calculated 
Attenuation  Coefficient  for  Cirrus.  The  experimental 
data  curve  is  chosen  for  the  LOW  TRAN  model,  as 
explained  in  the  text 


making  the  measurements  that  lead  to  the  experimental  curve.  We  thus  have  an 
equation  for  cirrus  normal  transmittance,  r  ,  of  the  form 

r  =  exp  -  (0.  14  L2)  .  <49) 

40 

This  expression  closely  duplicates  the  double  exponential  model  of  Davis  for 
most  realistic  cirrus  thicknesses.  Even  though  a  constant  attenuation  coefficient 
k  =  0.  2  km  1  is  nearly  as  good  a  fit  to  the  experimental  data,  the  LWC  dependence 
on  cloud  thickness  argues  in  favor  of  the  adopted  expression. 


7.7  Cirrus  Computer  Subroutine 


A  cirrus  subroutine  (CIRRUS)  was  added  to  the  DOWTRAN  6  code  to  include 
the  extinction  due  to  cirrus  clouds.  As  indicated  in  the  user  instructions  in 
Section  10,  the  cirrus  subroutine  is  called  by  setting  ICIR  =  l  on  CARD  2.  The 
default  value  for  the  cirrus  thickness  is  the  median  value  of  1  km  and  for  the 
cirrus  base  altitude  is  11,  10,  8,  7,  and  5  km  for  the  tropical,  midlatitude  sum¬ 
mer,  midlatitude  winter,  subarctic  summer,  and  subarctic  winter  models  respec¬ 
tively.  The  user  may  override  the  default  values  by  specifying  CTIIIK  (cirrus 
thickness)  and  CALT  (cirrus  base  altitude). 

If  the  user  sets  the  random  number  initialization  flag  ISEED  f  0.  the  program 
will  generate  random  cirrus  thicknesses  according  to  the  probability  distribution 
shown  in  Figure  19.  (The  appropriate  computer  system  random  number  generator 
must  be  included.  )  This  feature  may  be  of  use  in  electro-optical  systems  studies 
to  establish  transmittance  probability  statistics.  To  use  this  option,  the  program 
must  be  run  repeatedly  using  the  control  parameter  IRPT  >  0  on  CARD  4. 


8.  RAIN  MODEL 

8.1  Intro<Iuct:on 

In  developing  a  rain  model  to  be  used  in  the  LOW  TRAN  code  the  first  con¬ 
sideration  was  relating  the  transmission  over  a  given  path  to  the  most  directly 
obtainable  meteorological  parameter.  This  parameter  is  the  rain  rate  in  mm/h, 
reported  by  worldwide  weather  stations  on  a  six  hourly  basis. 

The  Marshall-Palmer  (M-P)f>()  raindrop  size  distribution  was  chosen  because 
the  two  main  components  are  rain  rate  and  drop  diameter,  and  the  M-P  raindrop 

size  distribution  is  widely  accepted  in  the  research  community.  The  M-P  distribu- 

O'  7 

tion  is  the  same  one  being  used  in  the  millimeter  region  (Falcone  et  al;  ‘  by 
the  FASCODE  (Clough  et  al)'  high  -resolution  atmospheric  transmittance/ radiance 
modeling  code. 


88.  Marshall,  J.S.,  and  Palmer,  W.  M.  K.  (1948)  The  distribution  of  raindrops 
with  size,  J,  Meteorol.  5:165-188, 

67,  Falcone,  Jr.,  V.J.,  Abreu,  L.  W.  ,  and  Shuttle,  E.F.  ( l  ;)7  D)  Atmospheric 

Attenuation  of  Millimeter  and  Submillimeter  Waves:  Models  and  Computer 
Code,  AFGL-TR -7  9-0258,  AD  A084485. 
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8.2  Formulation  of  the  Mode! 


The  M-P  dropsize  distribution  is  given  in  Eq.  (50). 


dN 

dD 


n  (D)  =  exp  (-AD) 


(50) 


where 


N  ~  8,000.  mm  ^  m  ^  (5 la) 

A  =  4.  1  R'0-21  (51b) 

R  =  rain  rate  (mm  hr  *) 

D  =  drop  diameter  (mm) 

From  Mie  theory  we  can  write  the  extinction  coefficient,  kext: 

00 

kext  =  /  f  D*  QextUD/X'  mM]  dD  dD  *  <52) 

o 

where  Qgxt  is  the  Mie  Extinction  Efficiency,  A  is  the  wavelength,  and  m(A)  is  the 
complex  refractive  index  of  water.  Since  for  rain  0.  1  <  D  <  5  mm,  in  the  visible 
and  infrared  we  have  (D  »  A).  Therefore,  Q  ~  2,  independent  of  the  wave¬ 
length.  Using  this  assumption  and  Eq.  (50)  in  Eq.  (52)  we  have: 


00 

kext  '  f  No  /  U“  eJipCAL>)  dD  . 
o 

Carrying  out  the  integration,  this  simplifies  to: 


k  .  7T  N  A 
ext  o 


Substituting  Eq.  (51)  in  Eq.  (54)  yields: 


kext  *  °-365  U°'  G3  (k'n"1) 


(53) 


(54) 


(55) 


08 


We  should  note  here  that  this  derivation  shows  that  the  extinction  due  to  rain 
is  independent  of  wavelength,  assuming 

\  «  D  ~  0.  1  to  10  mm 

In  practice  this  assumption  applies  throughout  the  visible  and  the  IR  windows. 

The  transmittance  over  path  length,  s,  in  km,  can  be  written  as 

r  =  exp  (-ks)  ,  (56) 

or  using  Eq.  (55) 

r  =  exp  (-0.  365  R0'  63  s)  .  (57) 

It  should  be  recognized  that  the  extinction  [Eq.  (55)]  or  transmittance 

[Eq.  (57)]  measured  by  a  transmissometer  will,  in  general,  have  to  be  corrected 

for  forward  scattering  effects.  However,  since  this  correction  is  a  function  of 

the  receiver  and  source  geometries  it  is  beyond  the  scope  of  the  current  LOWTRAN 
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model.  A  discussion  of  this  correction  is  given  by  Shettle  et  al. 

8.3  Other  Raindrop-Size  Distribution  Formulations 

Several  researchers  have  attempted  to  relate  the  parameters  in  the  exponen¬ 
tial  rain  dropsize  distribution  to  the  type  of  rainfall  (for  example,  Joss  and 
69  7  0 

Waldvogel,  a  and  Sektion  and  Srivastava  ).  These  different  parameterizations 
lead  to  an  expression  similar  to  Eq.  (55)  for  the  extinction: 

k  -  A  •  RB  .  (58) 

ext 

For  the  convenience  of  the  LOWTRAN  users  who  may  wish  to  modify 
SUBROUTINE  TNRAIN,  to  implement  one  of  these  other  models,  the  parameters 
for  size  distributions  and  extinction  coefficients  are  summarized  in  Table  8. 


68.  Shettle,  E.  P.  ,  Fenn,  R.W.,  and  Mill,  J.  D.  (1983)  The  Optical  and  Infrared 

Properties  of  Atmospheric  Particulates,  AFGL-TR-83-  (to  be 

published). 

69.  Joss,  J.  ,  and  Waldvogel,  A.  (1969)  Raindrop  size  distributions  and  sampling 

size  errors,  J.  Atmos.  Sci.  26:566-560- 

70.  Selchon,  R.S,  ,  and  Srivastava,  R.  C.  (107  l)  Doppler  radar  observations  of 

drop-size  distributions  in  a  thunderstorm,  J,  Atmos.  Sei,  215:983-984. 
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Table  8.  Parameters  Relating  Size  Distribution  [Eq.  (50)]  Extinction 
Coefficient  [Eq.  (58)]  to  Rain  Rate  for  Different  Types  of  Rain 


Type  of  Rain 

N  (mm  ^  m  ^ ) 
o 

A  (mm-1) 

A 

B 

Marshall-Palmer6  6 

8,  000 

4.  1 

r-0.2! 

0.  365 

0.  63 

Drizzle  (Joss  and 
Waldvogel)69 

30, 000 

5.7 

r-0.21 

0.  509 

0.  63 

Widespread  (Joss  and 
Waldvogel)”® 

7,000 

4.  1 

r-0.21 

0.  319 

0.  63 

Thunderstorm  (Joss  and 
Waldvogel)"" 

1,400 

3.  0 

r-0.21 

0.  163 

0.  63 

Thunderstorm  (Sekhon  and 
Srivastava)^0 

7, 000  R0,  37 

3.  8 

r-0.14 

0.  401 

0.  70 

The  divergence  between  the  two  different  thunderstorm  models  indicates  the 
difficulty  in  making  such  parameterizations  and  the  uncertainty  in  the  parameter 
values  given. 

8.4  Sample  Output  of  Typical  Rain  Cases 

The  atmospheric  transmittance ,  using  the  M-P  model  with  the  LOWTRAN 
code  for  rain  rates  varying  from  1  to  100  mm/hr  is  shown  in  Figure  2  5  for 
400  through  4000  cm  1  and  in  Figure  26  for  4000  through  40000  cm  . 


WAVENUMBER  (CM-1) 


Figure  2  5.  Atmospheric  Transmittance  for  Different  Rainrates  (RR)  and 
for  Frequencies  From  400  to  4000  cm-1.  The  measurement  path  is  300  m 
at  the  surface  with  T  =  Td  =  10°C,  with  a  meteorological  range  of 
23  km  in  the  absence  of  rain 


WAVENUMBER  ( Cfl-U 


Figure  26.  Atmospheric  Transmittance  for  Different  Rainrates  (RR)  and 
for  Frequencies  I'rom  4000  to  40000  cm  .  The  measurement  path  is 
3  00  m  at  the  surface  with  T  =  Tdew  =  io°C,  with  a  meteorological  range 
of  23  km  in  the  absence  of  rain 


9.  DISCUSSION  OF  LOW  TRAN  6  CODE 


This  section  describes  the  LOWTRAN  6  program  structure,  including  the 
structure  of  several  modules  of  subroutines.  Five  tables  are  presented  that 
include  names  and  descriptions  of  each  subroutine  as  they  appear  in  the  program. 

Section  9.2  discusses  program  portability,  B’ortran  compatibility,  precision, 
and  some  specific  comments  relating  to  implementing  the  program.  Specific 
instructions  are  given  in  Section  9.  3  enabling  the  user  to  reduce  the  execution 
field  length.  Lastly,  the  availability  of  the  program  package  is  detailed  in 
Section  9.4. 

9.1  LOWTRAN  6  Code  Structure 

A  graphical  representation  of  the  LOWTRAN  6  main  program  structure  is 
depicted  in  Figure  27.  The  structure  for  non-standard  model,  air  mass,  single 
scattering  geometry,  and  transmittance  modules  are  shown  in  Figures  2  8  through 
31  respectively.  Descriptions  of  the  major  executable  subroutines  shown  in 
Figure  27  are  given  in  Table  9.  The  subroutines  in  the  air  mass,  single  scatter¬ 
ing  geometry,  and  transmittance  modules  are  explained  in  Tables  10,  11,  and  12 
respectively.  Table  13  contains  a  brief  description  of  the  Block  Data  subroutines. 


I  SUB-  I 

I  ROUTINES  t 
I _ i 

Figure  27.  LOWTRAN  6  Main  Program  Structure.  The  boxes  enclosed 
by  dashes  arc  modules  of  subroutines  for  the  calculation  of  non-standard 
models,  air  mass  geometry,  single  scattering  geometry,  and  trans¬ 
mittance 
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Figure  28.  Program  Structure  for  the  Non-standard 
Model  Subroutines 


Figure  2  9.  Program  Structure  for  the  Air  Mass  Geometry  Subroutines 
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Figure  30.  Program  Structure  for  the  Single  Scattering  Geometry  Subroutines 


Figure  31.  Program  Structure  for  the  Transmittance  Subroutines.  The  dashed 
box  labelled  Molecular  Attenuation  includes  the  following  subroutines:  ClDTA, 
C2DTA,  C3DTA,  C4DTA,  C6DTA,  SLF296,  SLF260,  FRN2  96.  and  HNQ3 
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Table  8.  Description  of  LWTRN6  Subroutines 


LWTRNG 

-  Main  driver  program.  Reads  control  cards. 

MARINE 

-  Determines  aerosol  extinction  and  absorption 
coefficients  for  the  Navy  maritime  model. 

CIRRUS 

-  Generates  altitude  profiles  of  cirrus  cloud  density. 

Irandom 

-  Calls  machine -dependent  function  RANF,  that  is  a 
uniform  random  number  generator. 

VSA 

-  Army  vertical  structure  algorithm  of  aerosol 
extinction  and  relative  humidity  for  low  visibility/ 
low  ceiling  conditions. 

NSMDL 

-  For  user-defined  atmospheric  or  aerosol  profiles. 

LAYVSA 

-  Restructures  the  atmospheric  profile  for  finer 
layering  near  the  ground  for  use  with  the  VSA  option. 

RDNSM 

-  Reads  model  7  data  for  Army  vertical  structure- 
algorithm. 

SUBSOE 

-  Calculates  the  subsolar  point  angles  based  upon 
time  and  day. 

STDIvIDL 

-  Sets  up  atmospheric  profiles  of  attenuator  densities. 

1 AERPRF 

-  Computes  scaling  factor  profiles  for  aerosols. 

GEO 

-  Driver  for  air  mass  subroutines.  Calculates 
attenuator  amounts  for  the  slant  path. 

SSGEO 

-  Obtains  attenuator  amounts  from  scattering  points 
along  optical  path  to  the  extraterrestrial  source. 

EXABIN 

-  Loads  aei’osol  extinction  and  absorption  coefficients 
for  the  appropriate  model  and  relative  humidity. 

TRANS 

-  Calculates  transmittance,  atmospheric  radiance, 
and  solar/lunar  scattered  radiance  for  slant  path. 
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Table  10.  Description  of  Air  Mass  Subroutines 


GEO 

-  Driver  for  air  mass  subroutines.  Calculates 
attenuator  amounts  for  the  slant  path. 

GEOINP 

-  Interprets  geometry  input  parameters  into  the 
standard  form  Hi,  H2,  ANGLE,  and  LEN. 

FNDHMN 

-  Calculates  HMIN,  the  minimum  altitude  along 
the  path  and  PHI,  the  zenith  angle  at  H2. 

REDUCE 

-  Eliminates  slant  path  segments  that  extend 
beyond  the  highest  profile  altitude. 

FDBETA 

-  Calculates  angle,  given  Hi,  H2,  and  BETA  by 
iteration. 

RFPATH 

1 

I 

-  Determines  the  refracted  path  and  the  absorber 
amounts  through  all  the  layers. 

FILL 

-  Defines  the  boundaries  of  the  slant  path  and 
interpolates  densities  at  these  boundaries. 

LAYER 

-  Calculates  the  path  and  amounts  through  one 
layer. 

RADREF 

-  Computes  radius  of  curvature  of  the  refracted 
ray  for  a  horizontal  path. 

FINDSH 

-  Finds  layer  boundaries  and  scale  height  at  ground 
for  index  of  refraction. 

SCALHT 

-  Calculates  scale  height  of  index  of  refraction. 

ANDEX 

-  Computes  index  of  refraction  at  a  specific  height. 

EX  PINT 

-  Performs  exponential  interpolations  for  the 
geometry  routines. 

Table  11.  Description  of  SSGEO  Subroutines 


SSGEO 

-  Obtains  attenuator  amounts  from  scattering  points  along 

optical  path  to  the  extraterrestrial  source. 

FSIDEL 

-  Calculates  the  relative  azimuth  between  the  line  of  sight 
and  the  direct  solar /lunar  path. 

PSI 

-  Returns  solar  azimuth  relative  to  line-of-sight  at 
current  scattering  location. 

DEL 

-  Returns  solar  zenith  angle  at  any  point  along  optical  path. 

GEO 

-  Driver  for  air  mass  subroutines.  Calculates  attenuator 
amounts  for  the  slant  path. 

SCTANG 

-  Returns  the  scattering  angle  at  any  point  along  the 
optical  path. 

9.2  Portability 

LOWTRAN  6  was  developed  on  a  CDC  6600  (a  60  bit-per-word  machine)  with 
the  FORTRAN  77  (FORTRAN  5)  compiler.  A  major  effort  has  been  made  to 
make  LOWTRAN  6  as  compatible  as  possible  with  both  ANSI  standard 
FORTRAN  66  (FORTRAN  4)  and  77,  and  to  run  in  single  precision  on  a  32  bit- 
per-word  computer.  For  the  most  part,  this  goal  has  been  achieved.  The  only 
non-ANSI  FORTRAN  77  feature  is  that  type  CHARACTER  is  not  used,  rather 
character  data  is  stored  as  Hollerith  data  in  real  variables,  four  characters 
per  word.  This  feature  is  compatible  with  FORTRAN  66  and  most  FORTRAN  77 
compilers  allow  it  as  well.  For  FORTRAN  66  compilers,  three  non-ANSI 
features  used  are  the  OPEN  statements,  named  BLOCK  DATA,  and  single  quotes 
in  FORMAT  statements.  It  should  be  a  very  simple  matter  to  modify  the  ap¬ 
propriate  statements. 

LOWTRAN  6  uses  three  files: 

1.  INPUT,  read  on  UNIT  =  5,  containing  LOWTRAN  input  directions. 
Maximum  record  length  is  80  characters. 

2.  OU'IPUT,  written  on  UNIT  =  6,  containing  the  standard  LOWTRAN  output. 
Maximum  record  length  is  12  0  characters. 

3.  TAPE? ,  written  on  UNIT  =  7,  containing  copies  of  the  input  cards  and  the 
spectral  results  (transmittance  and/or  radiance).  Maximum  record  length  is 
120  characters. 

The  unit  numbers  for  these  files  are  stored  in  the  variables  IRD(  =  5),  IPK(=6), 
and  IPU(=7),  which  are  carried  in  the  common  block  IFIL.  These  files  are  ac¬ 
cessed  using  OPEN  statements  in  the  main  program. 


Table  12.  Description  of  TRANS  Subroutines 


TRANS 

-  Calculates  transmittance,  atmospheric  radiance, 
and  solar/lunar  scattered  radiance  for  slant  path. 

AEREXT 

-  Interpolates  aerosol  attenuation  coefficients  to 
required  wavenumber. 

HN03 

-  Determines  nitric  acid  absorption  coefficient  at 
required  wavenumber. 

TRANFN 

-  Calculates  transmittance  for  ozone,  uniformly- 
mixed  gases,  and  water  vapor. 

SOURCE 

-  Contains  solar  intensity  data  and  calculates 
lunar  intensity. 

TNRAIN 

-  Calculates  transmittance  of  rain  as  a  function  of 
rain  rate  and  slant  range. 

SSRAD 

! 

-  Performs  the  layer  by  layer  single  scattering 
radiance  sum. 

PHASEF 

-  Chooses  correct  phase  function  based  on  relative 
humidity,  frequency,  scattering  angle,  and  model. 

INTERP 

-  Performs  linear  or  logarithmic  interpolation. 

PF 

-  Returns  the  appropriate  phase  function  from  the 
stored  database. 

ClDTA 

-  Returns  water  vapor  band  absorption  coefficient 
at  required  wavenumber. 

C2DTA 

-  Returns  uniformly-mixed  gases  absorption 
coefficient  at  required  wavenumber. 

C3DTA 

-  Returns  ozone  band  absorption  coefficient  at 
required  wavenumber. 

C4DTA 

-  Returns  Ng  continuum  absorption  coefficient  at 
required  wavenumber. 

C6DTA 

-  Returns  molecular  scattering  attenuation 
coefficient  at  required  wavenumber. 

C8DTA 

-  Returns  ozone  UV  and  visible  absorption 
coefficient  at  required  wavenumber. 

SLF2  96 

-  Loads  self -broadened  water  vapor  continuum  at  29G°K, 

SLF2  60 

-  Loads  self -broadened  water  vapor  continuum  at  260°K. 

FRN2  96 

i 

-  Loads  foreign -broadened  water  vapor  continuum  at  2  9(1  K. 

1 

1  SINT 

-  Performs  interpolation  for  water  vapor  continuum. 

7  8 


Table  13.  Description  of  Block  Data  Subi’outines 


MDTA 

-  Modei  atmospheric  data. 

TITLE 

-  Titles  for  outputs. 

PRFDTA 

-  Aerosol  profile  data. 

EXTDTA 

-  Aerosol  extinction  and  absorption  data. 

SF2  9  b 

-  Self-broadened  absorption  coefficients  for  water 
continuum  at  2  96°K. 

vapor 

SF260 

-  Self-broadened  absorption  coefficients  for  water 
continuum  at  260°K, 

vapor 

BFH20 

-  Foreign-broadened  absorption  coefficients  for  water 
vapor  continuum  at  2  96°K. 

TRFN 

-  LOW  TRAN  transmittance  functions. 

CiD 

-  Water  vapor  band  model  absorption  coefficients, 

C2D 

-  Uniformly-mixed  gases  band  model  absorption 
coefficients. 

C3D 

-  Ozone  band  model  absorption  coefficients. 

04  D 

-  Nitrogen  continuum  absorption  coefficients  and 
UV  ozone  absorption  coefficients. 

MARDTA 

-  Navy  marine  aerosol  extinction  and  absorption  data. 

PHSDTA 

-  70  Averaged  phase  functions  and  truth  table  identifying 
correct  phase  function. 

For  computers  with  virtual  memory,  the  program  must  be  compiled  with 
"GLOBAL  SAVE"  or  the  local  equivalent.  Otherwise,  information  relating  to 
the  phase  function  is  initialized  in  a  first  call  to  a  subroutine,  but  lost,  in  sub¬ 
sequent  calls. 

The  cirrus  cloud  model  includes  an  option  to  generate  cirrus  clouds  at 
random  altitudes.  This  option  calls  the  random  number  generator  function  sub¬ 
routine  RANF,  which  is  called  from  RANDOM.  The  user  will  have  to  supply  the 
local  equivalent  of  RANF. 

The  user  may  find  it  useful  to  have  the  date  and  time  printed  at  the  beginning 
of  each  case.  The  statements  required  to  do  this  are  commented  out.  in  the  main 
program  with  'CO'  in  columns  1  and  2.  The  subroutines  FDATE  and  FCLOt’K, 
which  return  the  date  and  time  in  ad  format,  are  commented  and  located  at  the 
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end  of  the  main  program  LWTRN  6.  The  user  will  have  to  modify  these  sub¬ 
routines  as  necessary. 

The  sample  output  included  on  the  LOW  TRAN  6  tape  was  generated  on  a 
CDC  6600  that  has  14  decimal  digits  of  precision.  The  output  for  these  same 
cases  generated  in  a  32  bit-per-word  machine  with  about  seven  decimal  digits  of 
precision  (for  example,  an  IBM  370  or  4341  or  a  VAX)  should  agree  with  the 
sample  to  within  about  four  decimal  digits.  This  precision  is  better  than  the 
accuracy  of  the  models,  so  that  going  to  DOUBLE  PRECISION  to  increase  the 
precision  would  produce  no  real  benefit.  When  calculating  the  radiance  where 
the  whole  path  is  at  high  altitude  (40  km  or  more),  the  truncation  error  of  a 
32  bit-per-word  machine  might  cause  the  emission  to  be  zero. 

When  run  on  an  IBM  370  or  434  1  using  the  FORTRAN  G  compiler,  sample 
case  5  (Section  ll)  generates  two  boundaries,  both  at  1.  8  km,  and  then  produces 
divide  checks  in  subroutine  LAYER.  However,  the  final  answers  agree  with  the 
sample  output  found  in  this  report.  We  have  not  found  a  solution  to  this  problem. 

9.3  Execution  Field  Length 

On  a  CDC  6600,  LOWTR AN  6  requires  136000g  of  core  memory  to  run 
( 144000g  to  load).  There  are  several  techniques  available  to  the  user  to  reduce 
the  executable  field  length  if  desired.  On  CDC  CYBER  systems,  the  program 
may  be  run  using  the  SEGMENT  loader,  which  effectively  creates  an  overlay 
structure.  File  2  of  the  LOWTRAN  f>  tape  (Section  9.4)  gives  the  SEGMENT 
loader  directives  that  reduce  the  field  length  to  less  than  7  5000g.  Users  of 
other  systems  who  wish  to  create  overlays  should  consult  the  structure  chart  in 
Figure  27  and  the  SEGMENT  input  directives  for  guidance.  If  certain  program 
options  are  not  required,  the  field  length  can  also  be  reduced  by  not  including 
the  corresponding  subroutines.  For  example,  if  the  Solar  Scattering  option  is 
not  required  the  subroutines  SSGEO,  PSIDEL,  PSI,  DEL,  SCTANG,  SUBSOL, 
SSRAD,  PHASEF,  INTERP,  PF,  SOURCE,  and  PHSDTA  and  the  common  blocks 
SOLS  and  MNMPHS  can  be  deleted  saving  about  17u00g  of  field  length. 

9.4  Availability 

The  LOWTRAN  6  package  is  available  from: 

National  Climatic  Center,  NOAA 
Environmental  Data  Services 
Federal  Building 
Asheville,  NC  28801 

(704)  258-2850,  ext,  682  (Ms.  Yolanda  Goodge) 


The  package  is  normally  distributed  on  magnetic  tape  with  the  following 
characteristics : 

1.  9  track,  ltiOO  BPI 

2 .  unlabeled 

3.  ASCII 

4.  fixed-length  record"  140  characters-per-record,  one  record  per  block. 

The  tape  has  12  files.  ne  contents  of  these  files  arc  as  follows: 

1.  LOWTRAN  6  source  code 

2.  SEGMENT  input  direction  (relevant  to  CDC/NOS/BE  systems  only) 

3.  INPUT  file  for  test  cases 

4.  OUTPUT  file  for  test  cases 

5.  TAPE7  file  for  test  cases 

6.  Source  code  for  plot  program  (see  Appendix  A) 

7.  INPUT  file  for  test  cases  of  plot  program 

8.  OUTPUT  file  for  test  cases  of  plot  program 

9.  Source  code  for  filter  program  (see  Appendix  B) 

10.  INPUT  file  for  test  cases  of  filter  program 

11.  OUTPUT  file  for  test  cases  of  filter  program 

12.  List,  of  all  LOWTRAN  phase  functions. 

The  record  length  for  all  the  files  is  140  characters;  however,  for  files 
containing  source  code,  INPUT  files,  and  the  phase  functions,  only  the  first 
80  characters  contain  information.  Users  who  require  a  different  format  should 
contact  the  National  Climatic  Center.  The  tape  presently  costs  $99.  00. 

10.  INSTRUCTIONS  FOR  USING  LOWTRAN  6 

The  instructions  for  using  LOWTRAN  6  are  similar  to  those  for  previous 
versions.  However,  mati)  new  parameters  have  been  added  necessitating  the 
addition  of  a  fifth  control  card.  The  new  parameters  are  principally  required 
for  the  Navy  maritime  model.  Army  Vertical  Structure  Algorithm  (VSA), 

NOAA  Cirrus  Cloud  model,  and  inclusion  of  the  Solar/Lunar  Single  Scattering 
model. 

In  general,  foi  standard  atmospheric  models,  five  input  cards  are  required 
to  run  the  program  for  a  given  problem.  For  a  specific  problem  utilizing  any  of 
the  aforementioned  algorithms,  a  combination  of  several  of  the  ten  additional 
optional  control  cards  arc  possible.  The  formats  for  the  five  main  cards,  ten 
optional  cards,  and  definitions  of  the  input  p.arameters  are  given  below. 


10.1  Input  Data  and  Formats 

The  data  necessary  to  specify  a  given  problem  are  given  on  the  input  cards 
as  follows: 


CARD  1  MODEL,  ITYPE,  IEMSCT,  Ml,  M2,  M3,  IM,  NOPRT,  TBOUND,  SALB 
FORMAT (81 5,  2F10.  3) 


CARD  2  1HAZE,  ISEA3N,  IVULCN,  ICSTL.  ICIR,  IVSA,  VIS,  WSS,  WHH,  RAINRT 
FORMAT (615, 4F 10.  3) 


CARD  2 A  CTHIK,  CALT,  ISEED  (ICIR  =  l) 

FORMAT(2FlO.  3,110) 

CAP.D2B  ZCVSA ,  ZTVSA,  ZINVSA  (IVSA  =  l) 

FOR  MAT  (3F 10.  3) 


CARD  2 C],  ML,  TITLE 

FOR  MAT  (1 5,  18A4) 


(MODEL,  =  7) 
(IM  =  1) 


CARD  2C2(ML  Cards)  Z,  P,  T,  DP,  RH,  WH,  WO,  AKAZE,  VISl, 

I  HA  1 ,  ISEA 1,  IVUL1 

FOR MAT(3F  10.  3,  2F5.  1,  3E 10.  3,  F7 .3,  311) 

CARD  2D(10  Cards)  (DUMMY,  EXTC(l,I),  ABSC(l,I),  1=1,40) 

(I  HAZE  =  7) 

FORMAT  (4(F6.2,2F7.  5)) 


CARD  3  Hi,  H2,  ANGLE,  RANGE,  BETA,  RO,  LEN 
FORMAT  (6F  10.3, 15) 


ALTERNATE  CARD  3  (MODEL  =  0,  ITYPE  =  l) 


CARD  3  III,  P,  T,  DP,  RH,  WH,  WO,  RANGE 

FORMAT(3F  10.  3.2F5.  l,2ElO.  3,FlO.  3) 

CARD  3Al  IPARM,  IPII,  IDAY,  ISOURC  (IEMSCT  =  2) 
FOR  MAT  (415) 


CARD  3A2  PARMl,  PAR  M2,  PARM3,  PARM.4,  TIME,  PSIPO,  ANGLEM,  G 
(IEMSCT  =  2) 

FOR MA T (8 F 1 0 .  3) 

CARD  3B1  NANGLS  (IPII  =  l) 

FORMATU5) 

CARD  3B2  (1  TO  NANGLS)  (IPII  =  1) 

(ANGF(I),  F(l,  I),  F (2 ,  I),  F(3, 1),  F(4, 1),  1=1,  NANGLS) 
FORMAT(F  10.  3,  4E10.  3) 


CARD  4  Vl,  V2,  DV 

FOR  MA  T  (3F 10 .  3) 

CA  R  D  5  1R  PT 

FOR  MAT  (1 5) 


Definitions  of  these  quantities  will  be  discussed  in  Section  10.2. 


82 


10.2  Uasic  Instructions 


The  various  quantities  to  be  specified  on  each  of  the  five  control  cards  along 
with  the  ten  possible  optional  cards  (summarized  in  Section  10.  l)  will  be  dis¬ 
cussed  in  this  section. 

10.2.1  CARD  l  MODEL,  ITYPE,  IEMSCT,  Ml,  M2,  M3,  IM,  NOPR.T, 
TBOUND,  SALB 

The  parameter  MODEL  selects  one  of  six  geographical  model  atmospheres 
or  specifies  that  user-defined  meteorological  data  are  to  be  used  in  place  of  the 
standard  models.  ITYPE  defines  one  of  three  types  of  atmospheric  paths  for  a 
given  problem.  IEMSCT  selects  the  mode  of  program  execution  (transmittance, 
radiance,  or  radiance  with  solar/lunar  scattering).  Ml,  M2,  M3,  IM,  SALB, 
and  TBOUND  are  additional  input  parameters  for  non-standard  cases.  NOPRT 
is  a  user  option  to  suppress  printing  of  profiles  and  tables  in  the  output. 

MODEL  =  0  if  meteorological  data  are  specified  (for  horizontal  paths  only). 

=  1  selects  TROPICAL  MODEL  ATMOSPHERE. 

-  2  selects  MIDLATITUDE  SUMMER. 

=  3  selects  MIDLATITUDE  WINTER. 

=  4  selects  SUBARCTIC  SUMMER. 

=  5  selects  SUBARCTIC  WINTER. 

=  G  selects  1962  U.S.  STANDARD. 

=  7  if  a  new  model  atmosphere  (or  radiosonde  data)  is  to  be  inserted. 

ITYPE  indicates  the  type  of  atmospheric  path. 

ITYPE  =  1  for  a  horizontal  (constant-pressure)  path, 

=  2  for  a  vertical  or  slant  path  between  two  altitudes. 
s  3  for  a  vertical  or  slant  path  to  space. 

IEMSCT  determines  the  mode  of  execution  of  the  program. 

IEMSCT  =  0  program  calculates  transmittance. 

=  1  program  calculates  atmospheric  radiance. 

=  2  program  calculates  atmospheric  and  single  scattered  solar/Tunar 
radiance. 

=  3  program  calculates  directly  transmitted  solar  irradiance. 

Ml,  M2,  and  M3  are  used  to  modify  or  supplement  the  altitude  profiles  of  tem¬ 
perature  and  pressure,  water  vapor,  and  ozone  respectively. 

Ml  =  M2  -  M3  =  0  for  normal  operation  of  program. 


Ml  =  1  selects  the  TROPICAL  temperature  and  pressure  altitude  profiles. 

=  2  selects  the  MIDLATITUDE  SUMMER  temperature  and  pressure 
altitude  profiles. 


=  6  selects  the  1962  U.S.  STANDARD  temperature  and  pressure  altitude 
prof.' -es. 

M2  =  1  selects  the  TROPICAL  water  vapor  altitude  profile. 

=  2  selects  the  MIDLATITUDE  SUMMER  water  vapor  altitude  profile. 


=  6  selects  the  1962  U.S.  STANDARD  water  vapor  altitude  profile. 

M3  =  1  selects  the  TROPICAL  ozone  altitude  profile. 

=  2  selects  the  MIDLATITUDE  SUMMER  ozone  altitude  profile. 


=  6  selects  the  1962  U.S.  STANDARD  ozone  altitude  profile. 

IM  =  0  for  normal  operation  of  program  or  when  subsequent  calculations 
are  to  be  run  with  MODEL  =  7, 

=  1  when  radiosonde  data  are  to  be  read  in  initially. 

NOPRT  =  0  for  normal  operation  of  program. 

-  1  to  suppress  printing  of  the  atmospheric  profiles  and  the  ray 
trace  output. 

TBOUND  =  boundary  temperature  (°K),  used  in  the  radiance  mode  (IEMSCT 
=  1  or  2)  for  slant  paths  that  intersect  the  earth  or  terminate  at 
a  grey  boundary  (for  example  cloud,  target).  If  TBOUND  is 
left  blank  and  the  path  intersects  the  earth,  the  program  will 
use  the  temperature  of  the  first  atmospheric  level  as  the 
boundary  temperature. 

SA.LB  =  surface  albedo  of  the  earth  at  the  location  and  average  frequency 
of  the  calculation  (0.  0  to  1.  0).  If  SALB  is  left  blank,  the  pro¬ 
gram  assumes  the  surface  is  a  blackbody. 

Table  14  summarizes  the  use  of  the  four  control  parameters  MODEL, 

ITYPE,  IEMSCT,  and  NOPRT  on  CARD  l. 

10.2.2  CARD  2  IIIAZE,  ISEASN,  IVULCN,  ICSTL,  I  OR,  I.VSA,  VIS,  WSS, 
W11H,  RAINRT 

IIIAZE,  ISEASN,  IVULCN,  and  VIS  select  the  altitude-  and  seasonal- 
dependent  aerosol  profiles  and  aerosol  extinction  coefficients.  IIIAZE  specifies 
the  aerosol  model  used  for  the  boundary -layer  (0  to  2  km)  and  a  default-surface 
meteorological  range.  The  relative  humidify  dependence  of  the  boundary-layer 
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Table  14.  LOWTRAN  CARD  1  Input  Parameters:  MODEL,  ITYPE,  IEMSCT 
and  NOPRT 


MODEL,  ITYPE,  IEMSCT,  Ml,  M2,  M3,  IM,  NOPRT, 
CARD  1  TBOUND,  SALB 

FORMAT(8I5,2F10.  3) 


COL 

5 

MODEL 

0 

User- 

defined* 

1 

Tropical 

2 

Midlatitude 

summer 

3 

Midlatitude 

winter 

4 

Subarctic 

summer 

5 

Subarctic 

winter 

6 

1962  U.S. 
standard 

7 

User- 

defined* 

ITYPE 


Horizontal 

path 


Slant  path 
Hi  to  H2 


Slant  path 
to  space 


COL 

15 

IEMSCT 

0 

Transmittance 

1 

Radiance 

2 

Radiance  with 

solar/lunar 

scattering 

3 

Transmitted 
solar  irra- 
diance 

COL  NOPRT 
40 


Normal 

output 


Ml,  M2,  M3,  IM,  TBOUND,  and  SALB  are  leit  blank  for  standard  cases, 
"  Options  for  non-standard  models. 

Refer  to  Section  10.3  for  non-standard  cases. 


aerosol  extinction  coefficients  is  based  on  the  water  vapor  content  of  the  model 
atmosphere  selected  by  MODEL.  ISEASN  selects  the  seasonal  dependence  of  the 
profiles  for  both  the  tropospheric  (2  to  10  km)  and  stratospheric  (10  to  30  km) 
aerosols.  IVULCN  is  used  to  select,  both  the  profile  and  extinction  type  for  the 
stratospheric  aerosols  and  to  determine  transition  profiles  above  the  stratosphere 
to  100  km.  VIS,  the  meteorological  range,  when  specified,  will  supersede  the 
default  meteorological  range  in  the  boundary-layer  aerosol  profile  set  by  1HAZE. 


IHAZE  -  0  no  aerosol  attenuation  included  in  the  calculation. 

=  1  RURAL  extinction,  default  VIS  =  23  km. 

=  2  RURAL  extinction,  default  VIS  =  5  km. 

=  3  NAVY  MARITIME  extinction,  sets  own  VIS. 

=  4  MARITIME  extinction,  default  VIS  =  23  km  (LOWTRAN  5  model). 

=  5  URBAN  extinction,  default  VIS  =  5  km. 

=  6  TROPOSPHERIC  extinction,  default  VIS  =  50  km. 

=  7  USER -DEFINED  extinction,  default  VIS  =23  km  (read  into  the 
program  as  CARD2D), 

=  8  FOGl  (advection  fog)  extinction,  default  VIS  =  0.2  km. 

=  9  FOG2  (radiation  fog)  extinction,  default  VIS  =  0.5  km. 

As  noted,  IHAZE  selects  the  type  of  extinction  and  a  default  meteorological 
range  for  the  boundary-layer  aerosol  models  only.  If  VIS  is  also  specified  it  will 
override  the  default  IHAZE  value.  Interpolation  of  the  extinction  coefficients 
based  on  relative  humidity  is  performed  only  for  the  RURAL,  MARITIME, 
URBAN,  and  TROPOSPHERIC  coefficients  used  in  the  boundary  layer  (0-  to  2 -km 
altitude). 

ISEASN  =  0  season  determined  by  the  value  of  MODEL; 

SPRING-SUMMER  for  MODEL  =  0,  1,  2,  4,  6,  7 
FALL-WINTER  for  MODEL  =  3,  5 

=  1  SPRING -SUMMER 

=  2  FALL-WINTER 

ISEASN  selects  the  appropriate  seasonal  aerosol  profile  for  both  the  tropo¬ 
spheric  anu  stratospheric  aerosols.  Only  the  tropospheric  aerosol  extinction 
coefficients  are  used  with  the  2-  to  10-km  profiles. 

IVULCN  =  0,  1  BACKGROUND  STRATOSPHERIC  profile  and  extinction 

=  2  MODERATE  VOLCANIC  profile  and  AGED  VOLCANIC 
extinction 

=  3  HIGH  VOLCANIC  profile  and  FRESH  VOLCANIC  extmctiori 

=  4  HIGH  VOLCANIC  profile  and  AGED  VOLCANIC  extinction 

=  5  MODERATE  VOLCANIC  profile  and  FRESH  VOLCANIC 
extinction 

The  parameter  IVULCN  controls  both  the  selection  of  the  aerosol  profile  as 
well  as  the  type  of  extinction  for  the  stratospheric  aerosols.  It  also  selects 
appropriate  transition  profiles  above  the  stratosphere  to  100  km.  Meteoric  dust 
extinction  coefficients  are  always  used  for  altitudes  from  3  0  to  100  km. 
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VIS  =  surface  meteorological  range  (km)1'  (when  specified,  supersedes 
default  value  set  by  IHAZE). 

ICSTL  is  the  air  mass  character  (l  to  10),  only  used  with  the  Navy  maritime 
model  (IHAZE  =  3).  Default  value  is  3. 

ICSTL  =  l  open  ocean 


=  10  strong  continental  influence 

See  Section  5  for  a  detailed  description  of  determining  the  value  of  ICSTL, 

ICIR  selects  the  inclusion  of  cirrus  cloud  attenuation  in  the  calculation. 

ICIR  =  0  no  cirrus 

=  1  use  cirrus  profile 

See  Section  7  for  a  description  of  the  cirrus  cloud  model. 

IVSA  selects  the  use  of  the  Army  (VSA)  for  aerosols  in  the  boundary  layer. 

IV SA  =  0  not  used 

=  1  vertical  structure  algorithm 
See  Section  6  for  a  description  of  utilizing  the  Army  VSA. 

WSS  =  current  wind  speed  (m/s).  Only  used  with  the  Navy  maritime 
model  (IHAZE  =  3) 

WHH  =  24-h  average  wind  speed  (m/s).  Only  used  with  the  Navy  maritime 
model  (IHAZE  =  3) 


!,!The  terms  "meteorological  range"  and  "visibility"  are  not  always  used  correctly 
in  the  literature.  Correctly, r<  1»?2  visibility  is  the  greatest  distance  at  which  it 
is  just  possible  to  see  and  identify  with  the  unaided  eye:  (a)  in  the  daytime,  a 
dark  object  against  the  horizon  sky;  and  (b)  at  night,  a  known  moderately  intense 
light  source.  Meteorological  range  is  defined  quantitatively,  eliminating  the 
subjective  nature  of  the  observer  and  the  distinction  between  day  and  night. 
Meteorological  range  V  is  defined  by  the  Koschmieder  formula 


V 


3.  912 
(3  ' 


where  0  is  the  extinction  coefficient,  and  e  is  the  threshold  contrast,  set  equal 
to  0.  02.  As  used  in  the  LOWTRAN  computer  code,  the  inputs  are  in  terms  of 
meteorological  range,  with  0,  the  extinction  coefficient,  evaluated  at  0.55  urn. 
If  only  an  observer  visibility  V(^)S  is  available,  the  meteorological  range  can 
be  estimated  as  V  ~  (1.3  +  0.  3;  ’°V 

7  1.  Huschke,  R.E.,  Ed.  (1959)  Glossary  of  Meteorology,  American  Meteoro¬ 
logical  Society,  Boston,  Mass. 

72.  Middleton,  W.E.K.  (1952)  Vision  Through  the  Atmosphere,  University  of 
Toronto  Press,  2  50  pp. 
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If  WSS  =  WHH  -  0,  default  wind  speeds  set  by  MODEL. 


MODEL  WSS,  WHH  (m/s) 

0  6.  90 

1  4.  10 

2  4.  10 

3  10.29 

4  6. 69 

5  12.35 

6  -  7.20 

7  6.  90 

The  use  of  WSS  and  WHH  are  further  defined  in  Section  5. 

RAINRT  =  rain  rate  (mm/h).  Default  value  is  zero. 

The  use  of  the  parameter  RAINRT  is  explained  in  Section  8.  Table  15 
summarizes  the  use  of  aerosol  control  parameters  on  CARD  2. 

10.2.2.  1  Optional  Cards  Following  Card  2 

Optional  input  cards  after  CARD  2  selected  by  the  parameters  1CIR,  IVSA, 
MODEL,  and  I  HAZE  on  CARDS  1  and  2. 

CARD  2A  CTHIK,  CALT,  ISEED  (ICIR  =  l) 

Input  card  for  cirrus  altitude  profile  subroutine  when  ICIR  =  1. 

CTHIK  =  cirrus  thickness  (km) 

0  -  use  thickness  statistics 
,  NE.  0  =  user  defines  thickness 
CALT  =  cirrus  base  altitude  (km) 

0  =  use  calculated  value 
.NE.  0  =  user  defines  base  altitude 

ISEED  =  random  number  initialization  flag.  0  =  use  default  mean  values 
for  cirrus,  .NE,  0  =  initial  value  of  seed  for  RANF  function. 
Change  seed  value  each  run  for  different  random  number  se¬ 
quences.  This  provides  for  statistical  determination  of  cirrus 
base  altitude  and  thickness. 

CARD  2B  ZCVSA,  ZTVSA,  ZINVSA  (IVSA  =  l) 

Input  card  for  Army  VSA  subroutine  when  IVSA  =  l.  The  case  is  determined 

by  the  parameters  VIS,  ZCVSA,  ZTVSA,  and  ZINVSA. 

CASE  1:  cloud/fog  at  the  surface;  increasing  extinction  with  height  from 
cloud/fog  base  to  cloud/fog  top.  Selected  by  VIS  £  0.  5  km  and 
ZCVSA  5:  0. 

Use  case  2  or  2 '  below  cloud  and  case  1  inside  it. 


CASE  2:  hazy/lightly  foggy;  increasing  extinction  with  height  up  to  the 
cloudbase.  Selected  by  0.  5  <  VIS  ^  10  km,  ZCVSA  >  0. 

CASE  2':  clear/hazy;  increasing  extinction  with  height,  but  less  so  than 

case  2,  up  to  the  cloudbase.  Selected  by  VIS  >  10  km,  ZCVSA  >0. 

CASE  3:  no  cloud  ceiling  but  a  radiation  fog  or  an  inversion  or  boundary 
layer  present;  decreasing  extinction  with  height  up  to  the  height 
of  the  fog  or  layer.  Selected  by  ZCVSA  <  0,  ZINVSA  >  0. 

CASE  4:  no  cloud  ceiling  or  inversion  layer;  constant  extinction  with  height. 
Selected  by  ZCVSA  <  0  and  ZINVSA  <  0. 

ZCVSA  =  cloud  ceiling  height  (km): 

>0.0  the  known  cloud  ceiling  height; 

=  0.0  height  unknown:  the  program  will  calculate  one  for  case  2, 
and  default  is  1.  8  km  for  case  2';  or 

<  0.  0  no  cloud  ceiling  (cases  3  and  4). 

ZTVSA  =  the  thickness  of  the  cloud  (case  2)  or  the  thickness  of  the  fog  at 
the  surface  (case  l)  (km): 

>0.0  the  known  value  of  the  cloud  thickness; 

=  0.  0  thickness  unknown:  default  is  0.2  km. 


ZINVSA  =  the  height  of  the  inversion  or  boundary  layer  (km); 

>0.0  the  known  height  of  the  inversion  layer; 

=  0,0  height  unknown:  default  is  2  km,  0.  2  km  for  fog; 

<  0.  0  no  inversion  layer  (case  4,  if  ZCVSA  <0.0  also). 

CARD  2C 1  ML,  TITLE  (MODEL  =  7) 

(IM  =  l) 

A  new  atmospheric  profile  can  be  read  in  provided  the  parameters  MODEL 
and  IM  are  set  equal  to  7  and  1  respectively  on  CARD  i. 

ML  =  number  of  atmospheric  boundaries  to  be  read  in.  (maximum  of  34) 


TITLE  =  identification  of  new  atmospheric  profiles,  up  to  72  characters. 


CARD  2C2  (ML  cards) 

Z,  P,  T,  DP,  RII,  WII,  WO,  AHAZE,  VISl,  IHAl,  ISEAl,  IVULl 


Z 

P 

T 

DP 

RH 

WH 

WO 

AIIAZE 

VLSI 


altitude  (km) 

pressure  (mb) 

ambient  temperature  (°C) 

dewpoint  temperature  (°C) 

relative  humidity  (%) 

water  vapor  density  (gm  m  3) 

ozone  density  (gm  m  ^ ) 

aerosol  scaling  'actor  (normalized  by  the  user  to  the  required 

meteorology.  M  range  using  the  LOW  TRAN  extinction  coefficients) 

meteorological  range  (km)  for  the  altitude,  Z 


00 


IHAl  =  aerosol  extinction  and  meteorological  range  control  for  the 
altitude,  Z 

ISEA1  =  aerosol  season  control  for  the  altitude,  Z 

IVULl  -  aerosol  profile  and  extinction  control  for  the  altitude,  Z 

See  Section  10.  3  for  a  more  detailed  explanation  of  using  MODEL  =  7. 

CARD  2D  (10  cards)  (DUMMY,  EXTC(i,I),  ABSC(l,I),  1  =  1,40)  (I HAZE  =  7) 
User-defined  aerosol  extinction  and  absorption  coefficients  when  IIIAZE  =  7 
on  CARD  2. 

DUMMY  =  wavelength  of  aerosol  coefficient  (pm)  (not  used  by  program 
but  should  correspond  to  wavelengths  defined  in  array  VX2 
in  SUBROUTINE  EXTDTA) 

EXTC(l,I)  =  user-defined  aerosol  extinction  coefficient 
ABSC(l.I)  =  user-defined  aerosol  absorption  coefficient 

10.2.3  CARD  3  -  Hi,  H2,  ANGLE,  RANGE,  BETA,  RO,  LEN 

CARD  3  is  used  to  define  the  geometrical  path  parameters  for  a  given 
problem. 

Hi  =  initial  altitude  (km) 

H2  =  final  altitude  (km)  (for  1TYPE  =  2) 

H2  =  tangent  height  (km)  (for  ITYPE  =  3) 

It  is  important  to  emphasize  here  that  in  the  radiance  mode  of  program  execution 
(IEMSCT  =  1  or  2),  Hi,  the  initial  altitude,  always  defines  the  position  of  the 
observer  (or  sensor).  Hi  and  H2  cannot  be  used  interchangeably  as  in  the  trans¬ 
mittance  mode. 

ANGLE  =  initial  zenith  angle  (degrees)  as  measured  from  III 
RANGE  =  path  length  (km) 

BETA  =  earth  center  angle  subtended  by  III  and  H2  (degrees) 

RO  =  radius  of  the  earth  (1cm)  at.  the  particular  latitude  at  which  the 
calculation  is  to  be  performed. 

If  RO  is  left  blank,  the  program  will  use  the  midlatitude  value  of 
(1371.23  km  if  MODEL  is  set  equal  to  7.  Otherwise,  the  earth 
radius  for  the  appropriate  standard  model  atmosphere  (specified 
by  MODEL)  will  be  used. 

For  an  ITYPE  =  2  path  for  which  III  >  112  (and  by  necessity, 
ANGLE  >  90°),  two  paths  are  possible:  the  long  path  from  Hi 
through  the  tangent  height  to  H2  and  the  short  path  from  Hi  to  112. 
LEN  selects  the  type  of  path  in  these  cases. 

LEN  =  0  short  path  (default). 

=  1  long  path,  through  the  tangent  height. 

It  is  not  necessary  to  specify  every  quantity  given  above;  only 
those  that,  adequately  describe  the  problem  according  to  the 
parameter  ITYPE  (as  described  below'). 
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(1)  Horizontal  Paths  (ITYPE  =  ]) 

(a)  specify  Hi,  RANGE 

(b)  if  non-standard  meteorological  data  are  to  be  used,  that 
is,  if  MODEL  =  0  on  CARD  1,  then  refer  to  Section  10.  3 
for  parameters  and  format  of  CARD  3. 

(2)  Slant  Paths  Between  Two  Altitudes  (ITYPE  =  2) 

(a)  specify  Hi,  H2,  and  ANGLE 

(b)  specify  Hi,  ANGLE,  and  RANGE 

(c)  specify  Hi,  H2,  and  RANGE 

(d)  specify  Hi,  H2,  and  BETA 

(3)  Slant  Paths  to  Space  (ITYPE  =  3) 

(a)  specify  Hi  and  ANGLE 

(b)  specify  Hi  and  H2  (for  limb -viewing  problem  where  H2  is 
the  tangent  height  or  minimum  altitude  of  the  path 
trajectory). 

For  cases  2(b)  and  2(c.),  the  program  will  calculate  H2  and  ANGLE  respec¬ 
tively,  assuming  no  refraction;  then  proceed  as  for  case  2(a).  The  actual  slant 
path  range  will  differ  from  the  input  value.  This  method  of  defining  the  problem 
should  be  used  when  refraction  effects  are  not  important;  for  example,  for 
ranges  of  a  few  tens  of  km  at  zenith  angles  less  than  80°.  For  case  2(d),  the 
program  will  determine  the  proper  value  of  ANGLE  (including  the  effects  of 
refraction)  through  an  iterative  procedure.  This  method  can  be  used  when  the 
geometrical  configuration  of  the  source  and  receiver  is  known  accurately,  but 
the  initial  zenith  angle  is  not  known  precisely  due  to  atmospheric  refraction 
effects.  Beta  is  most  frequently  determined  by  the  user  from  ground  range 
information. 

Table  16  lists  the  options  on  CARD  3  provided  to  the  user  for  the  different 
types  of  atmospheric  paths. 

10.2.  3.  1  Alternate  CARD  3  for  Horizontal  Paths  (MODEL  =  0,  ITYPE  =  l) 

If  meteorological  data  are  to  be  used  for  horizontal  path  atmospheric  trans¬ 
mittance  calculations,  then  set  MODEL  =  0  on  CARD  1.  The  following  param¬ 
eters  can  then  be  specified  on  CARD  3, 

CARD  3  Hi,  P,  T,  DP,  EH,  WH,  WO,  RANGE 

Where  these  parameters  refer  to  altitude  (km),  pressure  (mb),  ambient 
temperature  (°C),  dewpoint  temperature  (°C),  relative  humidity  (%),  water 
vapor  density  (gm  m  ),  ozone  density  (gm  m  ),  and  path  length  (km) 
respectively. 

See  Section  10.  3  for  a  more  detailed  explanation  of  MODEL  =  0. 


Table  16.  Allowable  Combinations  of  Slant  Path  Parameters 


Case 

ITYPE 

HI  H2  ANGLE  RANGE  BETA  LEN 

(Optional) 

2a(l) 

2 

(2) 

2b' 

2 

=1=  * 

2e<3) 

2 

*  *  m 

2d<4) 

2 

5|<  i|<  >;< 

3a 

3 

>;<  )! t 

3b<51 

3 

(1)  LEN  option  is  available  only  when  Hi  >  H2  and  ANGLE  >  90°.  Otherwise, 
LEN  is  set  in  the  program. 

(2)  H2  calculated  assuming  no  refraction.  Calculated  RANGE  will  differ  from 
the  input  value. 

(3)  ANGLE  calculated  assuming  no  refraction.  Calculated  RANGE  will  differ 
from  the  input  value. 

(4)  Exact  ANGLE  is  calculated  by  iteration  of  the  path  calculation. 

(5)  H2  is  interpreted  as  the  tangent  height.  If  H2  and  ANGLE  are  both  zero. 
Case  3a  is  assumed  with  ANGLE  -  0  (that  is,  vertical  path).  For  a  path 
tangent  at  the  earth's  surface,  read  in  a  small  number  for  112,  for  example, 
0.  00 1  km. 


10.2.3,2  Alternate  CARD  3  for  Transmitted  Solar  Irradianee  (IEMSCT  =  3) 

For  calculating  directly  transmitted  solar  irradianee,  an  ITYPE  =  3  path  is 
assumed  and  CARD  3  has  the  following  form: 

Hi,  H2,  ANGLE,  IDAY,  HO 

FORMAT  (3F10.3,  15,  5X,  FlO.3) 

Hi  =  altitude  of  the  observer 

H2  -  tangent  height  of  path  to  sun 

ANGLE  =  apparent  solar  zenith  angle  at  Hi 

IDAY  -  day  of  the  year,  used  to  correct  for  variation  in  the  earth -to -sun 
distance 

RO  =  radius  of  earth  (default  according  to  MODEL) 

Either  H2  or  ANGLE  should  be  specified.  If  both  are  given  as  zero,  then  a 
vertical  path  (ANGLE  =  0°)  is  assumed.  If  I  DA  1  is  not  specified,  then  the  mean 
earth  to  son  distance  is  assumed. 
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If  the  apparent  solar  zenith  angle  is  not  known  for  a  particular  case,  then 
the  solar  scattering  option  (IEMSCT  =  2)  may  be  used  along  with  for  instance, 
the  observers  location,  day  of  the  year,  and  time  of  day  to  determine  the  solar 
zenith  angle  (see  Section  10.2.  3.  3  of  the  user  instructions).  Note  that  the 
apparent  solar  zenith  angle  is  zenith  angle  at  H 1  of  the  refracted  path  to  the  sun 
and  is  less  than  the  astronomical  solar  zenith  angle.  The  difference  between 
the  two  angles  is  negligible  for  angles  less  than  80°. 

10.2,  3.3  Optional  Cards  Following  CARD  3 

Optional  input  cards  after  CARD  3  are  selected  by  parameters  IEMSCT  on 
CARD  l  and  IPK  on  CARD  3A1. 

For  a  more  complete  description  of  the  single  scattering  input  parameters 
see  Section  4  and  Appendix  C. 

CARD  3A1  1PARM,  IPH,  IDAY,  ISOURC  (IEMSCT  =  2) 

Input  card  for  solar/lunar  scattered  radiation  when  IEMSCT  =  2. 

IPAIIM  =  0,  1,2  controls  the  method  of  specifying  the  solar/lunar  geometry 
on  CARD  3A2 . 

IPH  determines  the  type  of  phase  function  used  in  the  calculation. 

IPH  =  0  Henyey-Greenstein  aerosol  phase  function 

=  1  user-supplied  aerosol  phase  function  (see  CARD  3B) 

=  2  MIE -genera ted  database  of  aerosoi  phase  functions  for  the 
LOW  TRAN  models 

IDAY  =  day  of  the  year,  that,  is,  from  1  to  365  used  to  specify  earth  to  sun 
distance  and  (IPARM  =  l)  to  specify  the  sun's  location  in  the  sky. 

ISOURC  =  0  extraterrestrial  source  is  the  sun 

=  1  extraterrestrial  source  is  the  moon 

CARD  3A2  PARMl,  PARM2,  PARM3,  PARM4,  TIME,  PSlPO,  ANGLEM,  G 
(IEMSCT  =  2) 

Input  card  for  solar/lunar  scattered  radiation  when  IEMSCT  =  2.  Definitions 
of  PARMl,  PARM2,  PARM3,  PARM4  determined  by  value  of  IPARM  on  CARD 
3Al. 


For  IPARM  =  Q 

PARMl  =  observer  latitude  (-90°  to  +90°) 

Note  -  if  ARS(PARMl)  is  greater  than  89.  5°  the  observer  is 
assumed  to  be  at  either  the  north  or  the  south  pole.  In  this 
ease  the  path  azimuth  is  undefined.  The  direction  of  line -of - 
sight  must  be  specified  as  the  longitude  along  which  the  path 
lies.  This  quantity  rather  than  the  usual  azimuth  is  read  in. 

PAR  M2  =  observer  longitude  (0°  to  3<jU°) 
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PA I-5M 3  =  source  (sun  or  moon)  latitude 
PARM4  =  source  (sun  or  moon)  longitude 

For  1PAKM  =  1 

(IDAY  and  TIME  must  be  specified,  cannot  be  used  with  ISO  UK  (J  =  l) 

PARMi  -  observer  latitude  (-90°  to  +90°) 

PA  It  M2  =  observer  longitude  (0°  to  3(30°) 

PAR  M3,  PARM4  are  not  required 

For  IPAltM  =  2 

PARMi  =  azimuthal  angle  between  the  observer's  line -of -sight  and  the 

observer-to-sun  path,  measured  from  the  line  of  sight,  positive 
east  of  north,  between  -180°  and  180° 

PARM2  =  the  sun's  zenith  angle 

PARM3,  PARM4  are  not  required 

REMAINING  CONTROL  PARAMETERS 

TIME  =  Greenwich  time  in  decimal  hours,  that  is,  8:49  am  is  8.7  5, 

5:20  pm  is  17.  33  etc.  (used  with  IPAltM  =  l) 

PSIPO  =  path  azimuth  (degrees  east  of  north,  that  is,  due  north  is  0.  0°, 
due  east  is  90.0°  etc.  (used  with  IPAltM  =  0  or  l) 

ANGLEM  =  phase  angle  of  the  moon,  that  is,  the  angle  formed  by  the  sun, 
moon,  and  earth  (required  if  J.SOURC  =  l) 

G  =  asymmetry  factor  for  use  with  llcnycy-Greenstein  phase  function 

(used  with  1PH  =  o) 

CARD  3B1  NANGLS  (1PH  =  l) 

Input  card  for  user -defined  phase  functions  when  IPH  =  1. 

NANGLS  =  number  of  angles  for  the  user-defined  phase  functions  (maximum 
of  50) 

CARD  3B2  (1  to  NANGLS)  (IPH  =  l) 

(ANGF(I),  F(l,I),  F (2 ,1),  F(3,l),  F(4,I),  I  =  1,  NANGLS) 

Input,  card  for  user-defined  phase  functions  when  IPH  =  1. 

ANGF(I)  =  phase  angle  in  decimal  degrees  (0.0°  to  180.  0°) 

F ( l ,  I )  =  user-defined  phase  function  at  ANGF(I),  boundary  layer 

(0  to  2  km) 

F(2,I)  =  user-defined  phase  function  at  ANGF(I),  troposphere  (2  to  10  km) 

F(3,I)  =  user -defined  phase  function  at  ANGF(l),  stratosphere 

(10  to  30  lun) 

F(4,I)  =  user -defined  phase  function  at  ANGB’(r),  mesosphere 

(30  to  100  Km) 


10.2. 4  CARD  4  Vl,  V2,  DV 

The  spectral  range  and  increment  of  the  calculation. 

Vl  -  initial  frequency  in  wavenumber  (cm  *) 

V2  =  final  frequency  in  wavenumber  (cm-1),  where  V2  >  Vl 

DV  =  frequency  increment  (or  step  size)  (cm-1) 

(Note  that  v  =  104/X,  where  r  is  the  frequency  in  cm  1  and  X  is 
the  wavelength  in  pm,  and  that  DV  can  only  take  values  that  are 
a  multiple  of  5  cm-1  and  that  V  l  and  V 2  are  reset  to  the  next 
lowest  integer  multiple  of  5  cm-1.  ) 

10.2.5  CARD  5  IRPT 

The  control  parameter  IRPT  causes  the  program  to  recycle,  so  that  a  series 
of  problems  can  be  run  with  one  submission  of  LOWTRAN. 

IRPT  =  0  to  end  program 

=  1  to  read  all  data  CARDS  (l,  2,  3,  4,  5) 

=  2  not  used 

=  3  read  CARD  3  (the  geometry  card)  and  CARD  5 
=  4  read  CARD  4  (frequency)  and  CARD  5 
>  4  or  IRPT  =  2  will  cause  program  to  STOP 

Thus,  if  for  the  same  model  atmosphere  and  type  of  atmospheric  path  the 
reader  wishes  to  make  further  transmittance  calculations  in  different  spectral 
intex'vals  Vl1  to  V21  etc.  ,  and  for  a  different  step  size  (DV  1  etc.  ),  then  IRPT  is 
set  equal  to  4.  In  this  case,  the  card  sequence  is  as  follows  and  can  be  repeated 
as  many  times  as  required. 

CARD  5  IRPT  =  4 

CARD  6  Vl’  V2 ’  DV' 

CARD  7  IRPT  =  4 

CARD  8  Vl"  V2"  DV" 

CARD  9  IRPT  =  0 

The  final  IRPT  card  should  always  be  a  blank  or  zero.  When  using  the 
IRPT  =  4  option,  the  wavelength  dependence  of  the  refractive  index  is  not 
changed  (use  IRPT  =  l  option  if  this  is  required). 

Table  17  summarizes  the  user-control  parameters  on  CARD  4  am:  CARD  5. 

10.3  Non-Standard  Conditions 

Two  options  are  available  if  atmospheric  transm ittanec  calculations  are  re¬ 
quired  for  non-standard  conditions.  Here  non -standai’d  refem  to  conditions  other 
than  those  specified  by  the  six  model  atmospheres  provided  by'  LOWTRAN,  which 
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Table  17.  LOWTRAN  CARD  4  and  CARD  5  Input  Parameters:  VI,  V2,  DV, 
IRPT 


CARD  4 

Vl,  V2,  DV 

Format  (3F 10.  3) 

Vl  (cm  *)  V2  (cm  ^ )  DV  (cm  *)  Multiple  of  5  cm  * 

CARD  5 

IRPT 

Format  (15) 

COL 

5 

IRPT 

0 

End  of  program. 

1 

Read  new  CARDS  1,  2,  3,  4,  and  5. 

2 

Not  used  (stops  program). 

3 

Read  new  CARDS  3  and  5. 

4 

Read  new  CARDS  4  and  5. 

are  selected  by  the  parameter  MODEL  on  CARD  1.  These  options  enable  the 
user  to  insert: 

(1)  An  additional  atmospheric  model  (MODEL  7),  which  can  be  in  the  form 
of  radiosonde  data.  The  data  need  not  be  specified  at  the  same  altitudes  as  in 
the  standard  models. 

(2)  Meteorological  conditions  for  a  given  horizontal  path  calculation 
(MODEL  =  0  case). 

10.  3.  1  ADDITIONAL  ATMOSPHERIC  MODEL  (MODEL  =  7) 

A  new  model  atmosphere  can  be  inserted  as  CARD  2C  and  the  required 
multiples  ox  CARD  2 C,  provided  the  parameters  MODEL  and  IM  are  set  equal 
to  7  and  1  respectively  on  CARD  l,  The  number  of  atmospheric  levels  to  be 
inserted  (ML)  must  also  be  specified  on  CARD  2C.  Altitude -dependent  aerosol 
control  options  on  the  MODEJ ,  =  7  cards  provide  flexibility  to  the  user  in  model¬ 
ing  aerosol  extinction. 

The  appropriate  meteorological  parameters  and  format  for  the  atmospheric 
data  are  given  below: 

Z,  P,  T,  UP,  Rli,  WH,  WO,  AHAZE,  VISl,  IIIAl,  ISEAl,  IVULl 
FORMAT  (3E10.  3,  2E5.  1,  3E10.3,  F7.3,  311) 


Z  =  altitude  (km) 

P  =  pressure  (mb) 

T  =  ambient  temperature  (°C) 

DP  =  dewpoint  temperature  (°C) 

RH  =  relative  humidity  (%) 

WH  =  water  vapor  density  (gin  m 

WO  =  ozone  density  (gtn  m 

AMAZE  =  aerosol  number  density  (normalized  by  the  "«=-r  +n  the  veniiired 

meteorological  range  using  the  LOWTliAN  extinction  coefficients) 

VISl  =  meteorological  range  (km)  for  the  altitude,  Z 

IIlAl  =  aerosol  extinction  and  meteorological  range  control  for  the 
altitude,  Z 

ISEAl  =  aerosol  season  control  for1  the  altitude,  Z 

IVULl  =  aerosol  profile  and  extinction  control  for  the  altitude,  Z 

Note  that  it  is  necessary  to  specify  those  quantities  underlined  with  a  full  line 
and  only  one  of  the  three  quantities  underlined  with  the  dashed  line.  If  DP,  RH, 
and  WH  are  left  blank  then  a  dewpoint  of  0°C  is  assumed  unless  the  M2  option  on 
CARD  l  is  greater  than  0. 

If  the  ozone  density  (WO)  is  not  known  then  a  value  can  be  obtained  from  one 
of  the  standard  atmospheric  models  (for  the  appropriate  latitude  and  season)  by 
using  the  parameter  M3  on  CARD  l. 

Also  note  that  for  Ml  >  0  on  CARD  l,  both  pressure  and  temperature  are 
interpolated  from  the  model  atmosphere  (MODEL,  =  Ml)  for  the  altitude  Z. 

For  the  modeling  of  the  aerosol  profiles  and  extinction  coefficients,  if 
AHAZE,  VISl,  ISEAl,  and  IVULl  are  left  blank  on  the  MODEL  7  input  card, 
then  the  aerosol  control  parameters,  IHAZE,  ISEASN,  IVULCN,  and  VIS  on 
CARD  2  will  control  the  modeling  of  the  altitude-dependent  aerosol  parameters 
as  described  in  Section  10.2.  LOWTRAN  will  use  the  aerosol  models  contained 
in  the  program  and  interpolate  the  profiles  to  the  same  altitudes  as  the  radiosonde 
(or  new  model  atmosphere)  data. 

The  additional  aerosol  options  on  the  MODEL  7  card  provide  user  flexibility 
in  modeling  altitude -dependent  aerosols  such  as  low  ground  fogs,  where  finer 
altitude  resolution  is  required  to  specify  the  aerosol  profile.  These  options 
are  categorized  as  follows: 

(a)  AIIAZPJ  >  0,  VISl  =  HI  A  l  =  ISEAl  =  IVULl  =  0. 

For  this  case,  the  program  will  use  the  value  of  A1IAZF  at  the  altitude,  Z, 
to  define  the  aerosol  profile.  The  parameters  on  CARD  2  will  be  used  only  to 
select  the  type  of  aerosol  extinction  coefficients  to  be  used  in  the  0  to  2  km, 

2  to  10  km,  10  to  30  kin,  and  30  to  100  kin  altitude  regions  as  in  the 


MODEL  =  1  to  G  cases.  VIS  on  CARD  2  is  not  used.  The  user  must  scale  the 
A  HAZE  values  to  the  proper  sea  level  meteorological  range. 

(b)  A  HAZE  >  0,  cither  IIIA  ]  >  0  or  lVULl  >  0,  ISEAl  =  0, 
where  IHAl  =  1  to  9  with  the  same  extinction  coefficient  options  as  IHAZE  in 
Section  10.2,  and  IVULl  =  1  to  5  with  die  same  extinction  coefficient  options  as 
IVULCN  in  Section  10.2.  When  IIIA  1  is  defined,  it  will  select  the  type  of  extinc¬ 
tion  coefficient  to  be  used  with  AliAZE  at  the  altitude,  Z,  and  correspondingly 
when  IVULl  is  defined.  Only  four  different  altitude  regions  are  allowed  for  the 
aerosols  in  the  program.  The  boundary  altitudes  are  determined  from  the 
altitude,  Z,  on  the  MODEL  7  card  when  either  IHAl  or  IVULl  changes  value. 
These  boundaries  do  not  necessarily  have  to  correspond  to  the  default  values  in 
the  standard  models. 

(e)  AHAZE  =  0,  cither-  one  or  all  of  the  parameters  VISl,  IILAi,  ISEAl, 
and  IVULl  defined, 

where  ISEAl  =  l  or  2  with  the  same  seasonal  profile  options  as  ISEASN  in 
Section  10.2.  The  aerosol  profiles  and  extinction  coefficients  will  be  determined 
by  the  values  of  these  parameters  at  each  altitude  Z.  Again,  as  in  (b)  only  four 
altitude  regions  for  the  aerosols  are  allowed  in  the  program,  with  the  boundaries 
of  the  regions  determined  by  the  altitude  Z  when  the  control  parameters  change.. 
Note  also  that  IHAl  takes  precedence  over  IVULl  in  the  selection  of  the  type  of 
extinction  coefficients. 

Using  the  fog  models  provided  in  the  LOWTRAN  code,  these  additional 
aei-osol  options  allow  the  user  to  model  clouds. 

Note  that  IHAZE  must  be  defined  to  some  initial  value  greater  than  zero  to 
calculate  aerosol  extinction  and  that  at  least  two  altitudes  are  needed  to  define 
an  aerosol  altitude  region. 

When  using  the  Navy  maritime  aerosol  model  (IHAZE  =  3)  together  with  the 
Vertical  Structure  Algorithm  (VSA)  model,  some  difficulties  can  arise  if  run 
with  a  user-defined  atmosphere  (MODEL  =  7).  This  will  occur  when  the  user 
lets  the  maritime  aerosol  model  determine  the  meteorological  range  from  wind 
speed  and  relative  humidity.  The  problem  arises  because  when  MODEL  =  7,  the 
i.  .aritime  model  is  called  after  the  VSA  model,  which  requires  the  meteorological 
range  at  the  surface  to  be  specified.  To  get  around  this  the  user  must  run 
LOWTRAN  twice:  the  first  time  to  let  the  maritime  model  determine  the  mete¬ 
orological  range  (set  VIS  =  (j.  0  and  IVSA  =  o);  and  the  second  time,  the  actual 
run,  setting  VIS  equal,  to  the  meteorological  range  determined  by  the  maritime 
model  in  the  first  run. 
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10.3.2  HORIZONTAL  PATHS  (MODEL  =  0) 


If  meteorological  data  are  to  be  used  for  horizontal  path  atmospheric  trans¬ 
mittance  calculations,  then  set  MODEL  =  0  on  CARD  l.  The  following  param¬ 
eters  can  then  be  specified  on  CARD  3. 

CARD  3  HI,  _P,  T,  DP,  RH,  WH,  WO,  RANGE  (FORMAT  3Fl0.  3,  2F5.  1, 
2E10.3,  FlO.  3),  where  the  above  parameters  refer  to  altitude  (km),  pressure 
(mb)  ambient  temperature  (°C),  dewpoint  temperature  (°C),  relative  humidity 
(%),  water  vapor  density  (gm  m  °),  ozone  density  (gm  m  °),  and  path  length  (km) 
respectively. 

The  format  for  the  above  card  is  similar  to  that  used  to  input  radiosonde 
data  (MODEL  =  7).  Again,  it  is  only  necessary  to  specify  the  quantities  under¬ 
lined  with  the  solid  line  and  one  of  the  quantities  underlined  with  the  dashed  line. 
The  ozone  density  WO  can  be  specified  using  the  parameter  M3  on  CARD  1  if 
measurements  are  not  available.  In  the  latter  case,  a  value  will  be  calculated 
at  altitude  Hi  based  on  the  appropriate  model  atmosphere  selected  by  M3. 

The  aerosol  control  parameters  for  the  MODEL  =  0  cases  are  on  CARD  2  as 
described  in  Section  10.2. 

11.  EXAMPLES  OF  PROGRAM  OUTPUT 

Six  cases  representative  of  the  different  new  options  available  in  LOWTRAN 
6,  are  presented  in  this  section.  The  input  cards  for  the  six  cases  are  listed  in 
Table  18  and  the  listing  of  the  file  OUTPUT  (TAPE  6)  is  in  Tables  19  through  24. 
The  TAPE7  output  for  all  six  cases  is  in  Table  2  5. 

In  the  next  section,  the  output  for  case  1  will  be  discussed  in  detail.  For 
the  subsequent  cases,  only  the  details  relevant  to  the  new  features  will  be 
discussed. 

11.1  Cage  1:  Solar  Scattering 

This  case  demonstrates  the  new  air  mass  calculation  and  the  calculation  of 
the  singly-scattered  solar  radiance.  The  parameters  selected  for  this  case  are 
as  follows:  the  atmospheric  profile  is  the  1962  U.S.  Standard  and  the  boundary- 
layer  aerosol  model  is  RURAL  with  2 3 -km  meteorological  range.  The  path  is  a 
slant  path  from  500  km  to  the  ground  with  a  zenith  angle  at  500  km  of  160°.  The 
sun  is  used  as  the  source  and  the  sun  position  is  described  in  terms  of  the 
relative  path  azimuth  of  45°  and  the  solar  zenith  angle  of  60°.  The  earth-to-sun 
distance  is  that  for  day  1  (January  l).  The  aerosol  phase  functions  are  from  the 
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Table  20.  Program  Output  for  Case  2  (Contd) 
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Table  20.  Program  Output  for  Case  2  (Contd) 
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Table  21.  Program  Output  for  Case  3  (Contd) 


\  "S. 


oooo  ooooooouooo 


oeooooov?oooo  &  o  o 


0  3  0  0  0  0  v?  O  O  C’J  o  o  »:>  < 

l/'iD  lO  cD  ; 

o  o  o  o < 

tilt 
t.J  uJ  tu  U  l 

o  o  o  o  » 
en  —  o  t 
to  <7>  o  »-  < 

oooor.soooO‘»-i'-vrc4’ 

%r  ^  v  t)-  i  j  tf  i  ip 

o  o  o  o  o  t>  o  o  o 

i  <  i  I  I  I  l  I  l 

i.U  ij  10  uJ  uj  iO  u,*  uJ  uj 

OC-OOOOOOQ 

o  "•  cn  ^  o  o  m  >- 

CN  Ui  CO  O  v'M  O  O'  <*“)  CO 


03  <J  O  lO  y  tC  O  C  a)»-\TNOCVCUl 

<1  C  ^  n  PUN  Cr  O  O  01  o  ®  N  li11  lO 

n*  r-  cn  r-  cn  r-  cn  <n  y  u>  *-  to  •-  o 

w~  t-  cn  rv  n  n  t:  v  in  m  id  tc  n  b 


CO 

<1-  ^  CM  03 

in 

V- 

r- 

■c 

C) 

CD  ** 

TT 

f-  -- 

<T  r- 

Of  fN 

CN 

a  — 

u> 

C'  1/3 

in  ro  *- 

03 

cr  c- 

UT 

1-  y 

0 

CO  03  CD  r- 

CX  ♦- 

i/l  (U 

jn  r>,  in 

•*- 

▼- 

y 

y 

y 

T- 

»-  »- 

■»' 

£  5? 

UJ 

-y 

I/O 

0  < 

•-  -r-  •— 

T— 

T- 

■*-’ 

*- 

y  y 

y 

y  r- 

y  y 

1-  O 

I.t3 

r- 

<c 

»  *  * 

A 

• 

• 

4 

• 

•  • 

*  • 

2 

O 

a. 

UJ 

t_J  O 

»•• 

UJ  o, 

O 

cn 

3r. 

Iff 

0  0  O 

C' 

O 

O 

O  O 

O 

O  O 

O 

O  O 

0  0 

2 

C.1  Cl>  0 

CJ 

O 

O 

O 

-  j 

CJ 

O  O 

O 

O  O 

0  0 

m  1 

Cj  0 

O 

O 

O 

O 

O 

0 

O  O 

O 

O  O 

0  0 

n 

O  £f 

0  0  0 

O 

O 

O 

O 

0 

O  O 

O 

O  O 

0  0 

U.  CM 

Z  l- 

’  * 

in  in  to  <r  n  r* 


■(*3000000 


o  o  O  O  O  O  CV  o  o  c-  O  O  O  O  Ci 


o  O  O  o  c»  o  o  o  o  o  o  O  C3  o  o 

to  tr*  to  ir>  in  in  m  inters 
o ooooooooon 
t  lit  <  l  >  «  i  t  l 

uJ  uJ  UJ  U-J  UJ  Uj  lil  Ul  l;J  UJ  UJ 

i'-  n  ^  cv  cn  cr  cn  *-  «r  p- 

tr  pi  r-  m/ion  v-  rj  q  ifl  n 

(NO  PJ(N  -  IN  .N  CO  <T  CW 

tN0gW’MCN(N(Nr-P)r.(T)OOOO 

r^r>v-^tototOtotoiO'T*J-r3oto 

***  ^*  *'«-«-  T*  r-  y  «-  T-  T-  t-  y  o 

+  +-»■  +  +  +  -♦  +  +  +  •*■++  +  •♦■ 

III  UJ  UJ  LU  uj  LU  ul  UJ  UJ  UJ  UJ  UJ  Ul  UJ  111 

fl>  r»  ■  Lfi  rw  IO  O  in  to  *-  CN  <T  vr  to  in 

03  t-  o  r>3  0>  *-  v.  T  m  «T  •  <  —  O) 

f7-tOCJ’~  ^lfiOinaitDO«-CJU0P) 

y  y  T-v-tncocor-^CNin  —  f-roT- 


CO  r-  ?-  I  -  1C  O  iN  vf  •—  O  iu  iN  u)  O 
03  CM  to  o  <T  N  (jl 

Oi  o  o  "  t-  cn  n  ‘nn  «:■  r  r-  ^ 

t-  CN  (N  (N  (N  CN  (N  CN  tN  OKN  l-  CN  CN 


OOOOOOOOOOOC-^COO 
ooooooooooinoMDino 
0-10100030  O  t-  (N  O  O  lO  It)  o  o 

a)  totoa?  omoincNiDP)*" 
tp  to  nt  nj-  ro  ro  cm 


o  oooooooooooooo 

OOOO  OCOOOOOOOOO 


tc  0)  o  —  cn  ro  •r  iooiOoloooO 

*-  y  CNtNtNCNCNOJCnCOvj^TlOr^O 


(J)O«-NP)11ll(0Na)0)O.'WP) 

—  CMCN<NCN(NCNCNCNCN(NCOCJCOC*3 


■*■  V-  c0  ex'  U)  O)  nj  0(  O  r  -r  W  (T)  TJID 

Vr  /n  «-  n»  r-.  u)  p)  *.  o  a)  o  i  N  o  to  to 

to  *n  in  n  vr  tj  tj  v  k>  D  n  n  n  co  cm  cm 

c-  i'  t-  !•-  t-  f-  I"  r*  r-  c-  r-  r»  t-  r» 


o  o  o  o  o  o  o  o  o  o  o  o  o  ©  o  © 
O  r'  O  O  O  O  O  O  O  O  O  O  C  T*  ©  O 
0.  o  o  c*  o  o  o  o  o  o  o  o  o  o  o  o 
0000000000000000 


•—  o  *-•  tv  in  to  uj  (N  c£)  o  m  cr>  tn  o  c-  co 

in  m  vj-  ir  "3  r>»  cn  o  oi  p)  v  <d  of/  a»  •- 

—  y  —  y  —  —  -  !'j  (N  CN  (N  (V  (N  M  P) 

0000000000000000 


oooo 00000000 oooo 

o  o  o  o  o  o  o  o  o  o  o  o  t>  o  o  o 
c*  o  00  00000  00  00000 

OOOOOOOOOOOOOOOO 


O  C»  ©  C’J  O  O  O  O  O  <T-  <N  *-■  <r>(XJf-T- 

o  o  a  o  o  o  o  o  o  01  03  in  if:  vr  *-•  ^ 

O  O  O  r>  o  O  O  O  O  03  .73  03  (T:  <  J\  a>  it* 

oooooooooci'-tntncyitriatcn 


OiN'JMDnfniDMrtNOOtNOO 
O  ID  O  Ul  O'  y  *-  <J  (MD  S  MO  O  O 

o  at  03  (j)  00  t- c~  10  in  t  ^  *3  ■a- 

o  O'  c  o  cr  c  c  a-  o'  cr  cn  o  O'  o>  cr*  cr> 


O  Cl  Q  *7  OOOtPO)iJPil7|'JP)'-ID 

O)  (C)  (N  G  at  01  r-~  n  T-  —  CNID  O  ro  if!  CN 
r-  cr  pi  to  m  n  r-  1-  o  o  *- 
U)COCOOD03COa5COflDCD(DCOa?(CCCa) 


oitroipni'-inncD’'tD'-ouD(N'- 
U3i--a)ctioiO'~rNnrrin'jtTinif)N 
O  O  O  O  o  *- 

OOOO  O  OOO  00  o  00000 


’'Orficno-nina'fNin-'f'rtoiio 

»-  u>  a)  n  r->  .  1’  ff)  n  cd  in  n  -  tooin 
-  u  cji  o'  cd  ai  n  iptPininty'j'nncN 

■^■^-00000000000000 


o  in  o  m  o  in  o  in  o  if »  o  in  o  in  o  in 
oorr(\rcnnn  in  in  H)  (D  f-  P- 

O3CT><J30>O)ty>O>a>(J3cncnO3  0>O3O3j> 
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Table  22.  Program  Output  for  Case  4  (Contd) 
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Table  22.  Program  Output  for  Case  4  (Contd.) 
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Program  Output  for  Case  5  (Contd) 
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Table  24.  Program  Output  for  Case  6  (Contd) 
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Output  (Contd) 


Mie -generated  database.  The  surface  albedo  is  0.  05  and  the  surface  tempera¬ 
ture  is  300°K,  The  spectral  range  of  the  calculation  is  from  4000  to  4500  cm”* 
in  steps  of  10  cm 

The  output  for  this  case  shown  in  Table  19  will  now  be  described.  On  page 
one.  Table  19,  the  program  echoes  the  input  cards  exactly  as  they  are  read  in. 
CARD  l  selects  the  U.S.  Standard  Atmosphere  (MODEL  =  6),  type  of  slant  path 
(ITYPE  =  2),  single  scattering  option  (IEMSCT  =  2),  surface  temperature 
(TBOUND  =  300.  0),  and  surface  albedo  (SALE  =  0.  05).  CARD  2  selects  the 
RURAL  aerosol  profile  with  the  default  meteorological  range  of  2  3  km  (IHAZE  = 
l).  CARD  3  describes  the  slant  path  in  terms  of  the  observer  altitude  (Hi  - 
500  km),  endpoint  at  the  ground  (H2  =  0)  and  zenith  angle  (ANGLE  =  160.  0°).  On 
CARD  3Ai  1PARM  =  2  determines  that  the  sun  position  will  be  described  on 
CARD  3A2  in  terms  of  the  relative  azimuth  at  Hi  of  45°  (PARMl)  and  the  solar 
zenith  angle  at  Hi  of  60°  (PARM2).  IPH  =  2  specifies  that  the  Mie-generated 
phase  function  will  be  used,  while  IDAY  =  1  specifies  that  the  earth-to-sun 
distance  used  to  calculate  the  incident  solar  intensity  will  be  that  for  January  1. 
ISOURC  =  0  specifies  the  sun  as  the  source.  On  CARD  3A2,  only  PARMl  and 
PARM2,  as  described  earlier,  are  used  for  this  case.  The  remaining  param¬ 
eters  are  left  blank.  Card  4  selects  the  spectral  range  of  4000  to  4500  cm'A  in 
steps  of  10  cm  Following  CARD  4  the  program  interprets  the  input  cards  in 
a  more  readable  fashion. 

The  next  two  pages,  labeled  Atmospheric  Profiles,  list  the  profiles  of 
pressure,  temperature,  and  absorber  densities  at  the  indicated  altitudes  Z  for 
the  selected  atmospheric  model.  The  densities  of  the  absorbers  are  as  indicated 
or  are  in  relative  units  (-).  The  densities  of  water  vapor  (H^O),  uniformly- 
mixed  gases  (C02+),  ozone  (O^),  and  nitrogen  continuum  (N?)  are  in  scaled 
LOWTRAN  units.  'CNTMSLF1  and  'CNTMFRN'  refer  to  the  self-  and  foreign- 
broadened  water  vapor  continuum  density  factors.  The  column  labeled  'MOL 
SCAT1  is  the  total  air  density  relative  to  1013.25  mb  and  273.  1SK,  which  is  used 
to  calculate  molecular  scattering,  'N-l'  is  the  refractivity  or  the  index  of 
refraction  -1.0.  The  ultraviolet  ozone  [Oo(UV)]  density  is  the  actual  density  as 
is  the  nitric  acid  (HNO„)  density.  AEROSOL  l  to  4  refer  to  the  aerosol  scaling 
factors  for  the  four  aerosol  altitude  regimes.  'AERLRII'  is  AEROSOL  l  times 
the  relative  humidity  (RH)  and  is  used  to  calculate  the  path -averaged  RH  in  the 
boundary  layer.  ‘CIRRUS'  is  the  cirrus  cloud  density  profile  in  km  *  (see 
Section  7).  Since  the  cirrus  cloud  option  was  not  selected  for  this  case,  the 
cirrus  density  is  zero  everywhere. 

After  atmospheric  profiles,  messages  are  printed  concerning  the  interpreta¬ 
tion  of  the  slant  path  parameters.  In  this  case  the  observer  altitude  Hi  is  reset 
to  the  altitude  of  the  highest  atmospheric  profile  boundary  (100  km).  All  possible 
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combinations  of  slant  path  parameters  are  reduced  to  the  standard  set  of  Hi,  H2, 
ANGLE,  PHI,  IIMIN,  and  LEN,  where  PHI  is  the  arrival  angle  at  H2,  HMIN  is 
the  altitude  of  the  lowest  point  along  the  path,  and  LEN  has  a  value  of  unity  if 
the  path  goes  through  a  tangent  height  and  a  value  of  zero  otherwise. 

A  layer-by-layer  description  of  the  ray  trace  through  the  atmosphere 
appears  next.  The  ray  trace  always  begins  at  the  lowest  point  along  the  path 
HMIN  and  proceeds  upwards.  The  various  quantities  presented  are:  layer 
boundaries;  zenith  angle  at  the  bottom  of  the  layer  THETA;  curved  path  length 
through  the  layer  DRANGE;  cumulative  path  length  RANGE  from  HMIN;  earth- 
centered  angle  subtended  by  the  layer  DBETA  and  by  the  cumulative  path  BETA; 
arrival  angle  at  the  top  of  the  layer  PHI;  path  bending  within  the  layer  DBEND 
and  cumulative  BENDING;  density-weighted  path  averaged  pressure  and  tempera¬ 
ture  PBAR  and  TBAR  for  the  layer;  and  average  total  air  density  RHOBAR. 

The  next  table  lists  the  cumulative  absorber  amounts  for  the  species  de¬ 
scribed  in  the  atmospheric  profile  table  from  the  observer  to  the  different  layer 
boundaries.  'CNTMSLFl'  and  'CNTMSLF2'  refer  to  the  self -broadened  water 
vapor  amount  and  the  temperature  correction  factor  (Section  3).  A  summary  of 
the  geometry  calculation  appears  next,  followed  by  "Equivalent  Sea  Level  Total 
Absorber  Amounts"  that  are  the  total  amounts  for  the  whole  path  for  each  of  the 
forementioned  variables.  "Mean  RH"  refers  to  the  mean  relative  humidity  in 
the  boundary  layer  only. 

Next  is  the  output  of  the  single  solar  scattering  geometry  routines:  this 
includes  a  summary  of  the  path  from  each  scattering  point  to  the  sun,  'SUB¬ 
TENDED  ANGLE'  is  the  earth- centered  angle  between  the  observer  at  Hi  and 
the  scattering  point.  The  'SOLAR  ZENITH'  angle  is  the  astronomical  or  unre¬ 
fracted  zenith  angle  to  the  sun.  However,  if  the  refractive  bending  for  the  path 
to  the  sun  is  calculated  to  be  greater  than  0.  1°,  then  the  solar  zenith  angle  is 
corrected  by  the  bending  and  that  path  is  calculated  again  (up  to  four  times)  until 
the  correction  for  refraction  is  less  than  0.1°. 

The  'RELATIVE  AZIMUTH'  is  negative  because  internally  the  azimuth 
angles  are  measured  positive  counterclockwise,  while  for  the  input  cards, 
azimuth  angles  are  specified  positive  clockwise,  'SCTTR  ANGLE'  is  the  scat¬ 
tering  angle  at  that  point.  A  message  will  be  printed  if  any  scattering  point  is 
in  the  shade. 

Finally,  the  results  of  the  calculation  are  printed.  In  this  case,  it  is  the 
radiance  and  the  total  transmittance  for  the  path.  The  first  two  columns  are 
the  frequency  and  wavelength  in  wavenumbers  (cm  /  and  micrometers  (nm) 
respectively.  The  radiance  is  broken  down  by  source :  atmospheric  radiance 
(including  boundary  emission.',  path  scattered  radiance,  ground  reflected 
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radiance,  and  the  total  radiance,  each  listed  per  an  1  and  pm.  The  next  column 
"integral"  is  the  integrated  radiance  from  the  initial  frequency.  The  last,  column 
lists  the  total  transmittance  from  Ill  to  112.  A  short  summary  is  printed  at  the 
end  of  the  table,  including  the  actual  boundary  temperature  used  (if  any)  and  the 
minimum  and  maximum  values  of  the  total  radiance,  which  are  useful  in  deter¬ 
mining  the  limits  for  plotting. 

The  last  item  printed  is  the  value  of  IRPT  on  CARD  5.  In  this  case,  IRPT  = 

1  indicates  that  another  full  set  of  input  cards  follows. 

11.2  Case  2:  Cirrus  Cloud  Model 

Case  2  illustrates  the  cirrus  cloud  model.  The  cirrus  model  is  selected  by 
IC'IR  -  1  in  CARD  2  and  the  cirrus  parameters  are  read  in  on  CARD  2A.  Since 
the  parameters  CTHIK  and  CALT  are  input  as  zero,  the  default  values  appropriate 
to  the  Tropical  model  atmosphere  are  used.  The  message  regarding  the  prob¬ 
ability  of  clouds  occurring  is  informative  only,  it  has  no  effect  on  the  results. 

ITYPE  =  2  on  CARD  1  selects  a  slant  path  from  Hi  to  H2.  On  CARD  3  the 
parameters  specified  are  Hi  =  7  km,  H2  =  12  km,  and  P.ANGE  =  500  lun.  This 
value  of  RANGE  is  the  straight  line  distance  between  Ill  and  112  and  is  used  to 
calculate  the  value  of  ANGLE,  which  is  the  zenith  angle  at  Hi.  The  path  is  then 
traced  using  Hi,  H2,  and  ANGLE.  Because  of  refraction  the  calculated  slant 
path  of  545.094  km  is  greater  than  the  input  RANGE, 

Under  Atmospheric  Profiles,  cirrus  density  is  given  at  11.  0  km;  at  all  other 
altitudes  it  is  zero.  When  calculating  path  amounts,  the  program  will  linearly 
interpolate  the  cirrus  density  between  10  and  11  Ion,  and  between  11  and  12  km. 
The  transmittance  due  to  cirrus  is  independent  of  frequency  and  is  printed  only 
once,  at  the  end  of  the  table. 

11.3  Case  3:  Navy  Maritime  Aerosol  Model 

In  Case  3,  the  Navy  maritime  aerosol  model  is  selected  by  IIIAZE  =  3  on 
CARD  2,  Since  the  current  wind  speed  WSS  and  the  24-h  average  wind  speed 
WHH  are  left  as  zero,  these  parameters  default  to  the  values  appropriate  for  the 
tropical  model  atmosphere  (MODEL  =  1  on  CARD  1).  The  air  mass  character 
parameter  ICSTL  defaults  to  3  for  all  atmospheric  models.  The  meteorological 
range  VIS  is  determined  by  WSS,  WHH,  ICSTL,  and  the  relative  humidity.  The 
relative  humidity  in  the  boundary  layer  is  determined  by  the  water  vapor  density 
profile  for  the  selected  model  atmosphere.  If  a  value  fer  VIS  is  read  in  on 
CARD  2,  that  value  would  override  the  value  calculated  by  the  Navy  maritime 
model.  The  scaling  factors  for  the  Navy  maritime  model  are  listed  under 
Atmospheric  Profiles  as  AEROSOL  1. 
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11.4  Case  4;  Rain  Model 


Case  4  illustrates  the  use  of  the  rain  model,  whieh  is  selected  by  a  non-zero 
value  of  the  rain  rate  parameter  RA1NRT  on  CARD  2.  Selecting  MODEL  =  0  on 
CARD  l  alters  the  form  of  CARD  3,  which  now  requires  the  altitude  pressure, 
temperature,  water  vapor  and  ozone  amounts,  and  path  length  for  a  constant 
pressure  path.  The  transmittance  due  to  rain  is  independent  of  frequency  and 
is  printed  only  once,  at  the  end  of  the  transmittance  table. 

11.5  Case  5:  Vertical  Structure  Algorithm  (VS A) 

For  Case  5,  the  vertical  structure  algorithm  is  selected  by  specifying  IV SA 
on  CARD  2.  On  CARD  2B,  the  cloud  ceiling  ZCVSA,  cloud  thickness  ZTVSA, 
and  the  inversion  height  are  left  as  zero.  The  combination  of  VIS  =  5  km 
(determined  by  IHAZE  =  5:  urban  aerosol  model,  5-km  VIS)  and  ZCVSA  >  0 
selects  a  VSA  Case  2  boundary  layer  aerosol  profile.  The  cloud  ceiling  and 
cloud  depth  default  to  values  for  Case  2,  that  is,  1.8  and  0.2  km,  respectively. 

The  VSA  model  generates  a  boundary -layer  aerosol  profile  at  nine  altitudes 
between  the  ground  and  the  cloud  top  shown  in  Table  23.  This  profile  uses 
the  aerosol  type  shown  under  IHAZE  and  the  scaling  factor  shown  under 
EXTINCTION.  It  also  includes  the  relative  humidity,  which  is  set  to  100  per¬ 
cent  inside  a  cloud  or  fog. 

The  VSA  boundary-layer  aerosol  profile  is  incorporated  into  the  overall 
atmospheric  profile  by  creating  an  internal  MODEL  =  7.  The  VSA  profile  is 
written  to  an  array.  Then  the  overall  atmospheric  profile  is  added  to  that  array, 
starting  at  0.  01  km  above  the  top  of  the  VSA  profile.  Some  of  the  profile  bound¬ 
aries  are  eliminated  in  order  to  keep  the  total  number  of  profile  boundaries 
under  34.  This  array  is  then  read  in  by  the  subroutine  NSMDL  (non-standard 
model)  as  if  it  were  a  user -defined  profile  (MODEL  =  7).  Since  the  VSA  model 
defines  two  aerosol  types  in  the  boundary  layer  and  the  program  allows  at  most 
four  aerosol  regimes,  the  stratospheric  aerosol  regime  is  continued  up  to  the 
top  of  the  atmosphere  and  the  upper  atmospheric  regime  is  eliminated. 

11.6  Case  6:  Directly  Transmitted  Solar  Irradiance 

Case  6  calculates  the  directly  transmitted  solar  irradiance  from  6000  to 
7  500  cm  1  for  an  observer  at  the  ground  and  a  solar  zenith  angle  of  This 

case  is  selected  by  IEMSCT  -  3  on  CARD  1.  An  ITYPK  =  3  slant  path  is  re¬ 
quired:  CARD  3  specifies  the  observer  altitude  III  =  0.0,  the  apparent  solar 
zenith  angle  at  Hi,  ANGLE  =  00°,  and  the  day  of  the  year  1DAY  =  74  (the  Ides 
of  March). 
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The  spectral  output  consists  of  the  transmitted  solar  irradiance  and  the 
incident  solar  irradiance,  corrected  for  variations  in  the  earth-to-sun  distance 
according  to  the  day  of  the  year.  Also  given  are  the  cumulative  integrals  of 
these  two  quantities  and  the  total  transmittance  for  the  path  from  Hi  to  the  sun. 

11.7  Tape  7  Output 

The  spectral  results  from  each  successful  case  run  are  written  to  a  file  on 
UNIT  =  7,  along  with  a  10 -line  header  identifying  the  parameters  for  the  case. 
This  file  can  then  be  used  to  further  process  the  results.  For  example,  the  plot 
and  filter  programs  described  in  Appendices  A  and  B  work  from  this  file.  A 
listing  of  this  file  for  the  six  sample  cases  is  given  in  Table  2  5. 

The  data  on  the  file  consists  of  a  10-line  header  followed  by  the  spectral 
data,  The  header  consists  for  the  most  part  of  images  of  the  input  control  cards, 
with  the  default  values  supplied  for  parameters  that  have  defaults,  for  example 
VIS  on  CARD  2.  The  first  four  lines  of  the  header  are  CARD  1,  CARD  2, 

CARD  2A,  and  CARD  2B.  When  the  cirrus,  VSA,  and  single  scattering  options 
are  not  selected,  the  images  of  their  respective  control  cards  are  supplied  with 
the  parameters  set  equal  to  -99.  The  fifth  line  of  the  header  is  the  image  of 
CARD  2 Cl  and  gives  the  number  of  atmospheric  layers  and  title  of  the  atmos¬ 
pheric  profile.  This  data  is  given  for  all  values  of  MODEL  even  though  it  is 
read  in  only  foi  MODEL  -  7.  CARDS  2C2  and  2D  are  not  included  in  the  header. 
Line  6  contains  the  values  of  the  path  parameters  on  CARD  3  after  they  have 
been  evaluated.  CARDS  3Bl  and  3B2  are  not  included.  Lines  7  and  8  are  the 
images  of  CARDS  3Al  and  3A2,  with  all  the  parameters  set.  equal  to  -99  if  the 
single  scattering  option  is  not  selected.  Line  9  is  the  image  of  CARD  4,  the 
frequency  card,  and  line  10  gives  the  value  of  IRPT  from  CARD  5. 

Line  11  gives  the  values  of  the  transmittance  due  to  rain  and  cirrus  clouds, 
in  2F9.4  format.  Since  these  transmittances  are  independent  of  frequency,  they 
are  given  only  once  per  case. 

The  following  lines  contain  the  spectral  output  and  differ  depending  on  the 
type  of  calculation:  transmittance,  radiance,  radiance  with  single  scattering, 
or  directly  transmitted  solar  irradiance.  For  a  transmittance  calculation 
(IE MS C T  =  0)  each  line  of  spectral  data  has  the  following  format: 

1.  frequency  (cm  *)  1' 7 . 0 

2.  wavelength  (micrometer)  P8.3 

3.  transmittance,  total 

4.  transmittance,  water  vapor 

5.  transmittance,  uniformly  mixed  gases 

G,  transmittance,  ozone 
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10F8.  4 


7.  transmittance,  nitrogen  continuum 

8.  transmittance,  water  vapor  continuum 

9.  transmittance,  molecular  scattering 

10.  transmittance,  aerosol 

11.  transmittance,  nitric  acid 

12.  aerosol  absorption 

13.  integrated  absorption  {cm- *)  Fi2.4 


For  either  radiance  or  radiance  with  single  scattering,  the  format  is: 


1.  frequency  (cm 

2.  wavelength  (Mm) 

”2  -I  ~1  - 1,  - 

3.  atmospheric  radiance  (W  cm  ster  (cm  )  ) 

4.  atmospheric  radiance  (W  cm  ster  Mm  ) 

5.  path  scattered  radiance  (W  cm-2ster-  *  (cm  "  *)-  *) 

6.  path  scattered  radiance  (W  cm  ster  Mm  ) 

7.  ground  reflected  radiance  (W  cm-2ster-*(cm-*)-*) 

8.  ground  reflected  radiance  (W  cm-2ster-*  Mm'^) 

9.  total  radiance  (W  cm-2ster'  '^(cm  “*)“*) 

-2  - 1  -I 

10.  total  radiance  (W  cm  ster  Mm  ) 

11.  integrated  radiance  (W  cm-2ster-*)  _ 

12.  total  transmittance.  Hi  to  H2 

13.  TEBl 

14. -  TEB2SV 


F7 . 0 
F8.  3 


9E9.  2 


3F8.4 


If  the  single  scattering  option  is  not  selected  (that  is,  IEMSCT  =  l)  then  the 
path-scattered,  ground -reflected  AND  TOTAL  radiances  are  set  to  zero.  TEBl 
and  TEB2SV  are  also  defined  only  with  the  single  scattering  option:  TEBl  is  the 
total  transmittance  for  the  L-shaped  path  from  Hi  to  112  to  the  sun.  TEB2SV  is 
the  total  transmittance  for  the  path  from  Hi  to  the  last  boundary  before  H2  to 
the  sun. 

For  directly  transmitted  solar  irradiance,  the  format  is: 

1.  frequency  (cm- *) 
wavelength  (Mm) 

-2 

transmitted  solar  irradiance  (W  cm  (cm 
transmitted  solar  irradiance  (W  cm-2  Mm-*) 
incident  solar  irradiance  (W  cm -2(cm-  *)- *) 
incident  solar  irradiance  (W  cm-2  Mm-*) 
integrated  transmitted  irradiance  (W  cm-2) 


^  ■ 

3. 

4. 

5. 

fi. 

7. 

8. 

9. 


V1) 


integrated  incident  irradiance  (W  cm 
total  transmittance.  III  to  sun 


F7. 0 
F8.  3 


8E  9 .  2 


27X,  F9.  4 


The  end  of  the  spectral  data  is  marked  by  a  frequency  of  -9999. 
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Appendix  A 

LOWTRAN  6  Plot  Program 

The  LOWTRAN  6  plot  program,  LOWPLT,  is  an  independent  program  de¬ 
signed  to  plot  the  transmittance  and/or  radiance  output  of  a  LOWTRAN  6  run. 
The  program  will  plot  one  file  of  data,  several  files  of  data,  or  the  same  file  of 
data  as  many  times  as  desired.  Another  feature  of  the  plot  program  is  the  ability 
to  choose  any  of  the  many  variables  produced  by  a  transmittance  run,  radiance 
run,  or  a  radiance  with  scattering  run.  Finally,  the  user  is  given  the  option  of 
plotting  these  variables  on  individual  plots  or  on  the  same  plot.  When  plotting 
different  variables  on  the  same  plot,  it  is  important  to  define  the  proper  scale 
for  the  Y-axis.  The  LOWTRAN  6  plot  program  uses  Tape  7  as  generated  by  a 
LOWTRAN  6  run. 

This  plot  program  utilizes  plot  routines  unique  to  the  CDC-6G00  CALCOMP 
plotter.  Therefore  some  minor  changes  may  be  necessary  to  accommodate  a 
different  computer  plotting  system,  (See  Section  A3  for  more  details  on 
CALCOMP  system  differences. ) 

Al.  INSTRUCTIONS  FOR  USING  THE  LOWTRAN  6  PLOT  PROGRAM 

The  plot  program  extracts  the  data  to  be  plotted  from  Tape  7.  To  initiate 
the  program  five  input  cards  are  required. 
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Al.l  Input  Data  and  Formate 


The  data  necessary  to  specify  a  given  plot  are  given  by  the  following  five 
.  cards: 

CARD  1  PROGID,  SCALE 

(FORMAT  (3A10,  FlO.4)) 

CARD  2  XSIZE,  PFRBEG,  PFREND,  DELTAX,  ITYP,  IXAXIS,  NTJMFIL 
(FORMAT  (4F10.4,  315)) 

CARD  3  YSIZE,  YRM1N,  YRMAX,  DELTA Y,  ICRV,  IYAXIS,  NMYDEC 
I  (FORMAT  (FlO.4,  3E10.2,  315)) 

CARD  4  IRAD,  ITRP 

(FORMAT  (515)) 

CARD  5  ISAMFL,  ISAMPT 
(FORMAT  (215)) 

I 

A1.2  Utilizing  Input  Cards 


To  produce  multiple  plots;  input  cards  2  through  5  must  be  repeated.  To  end 
plotting;  the  first  value  on  the  following  input  card  2  should  be  negative. 
Definitions  of  the  input  card  variables  will  be  discussed  in  Section  A2. 

A2.  BASIC  INSTRUCTIONS 

The  various  quantities  to  be  specified  on  each  of  the  five  control  cards  are 
defined  in  this  section. 

A2.1  CARD  1  PROGID,  SCALE 

The  variable  PROGID  is  a  30  Hollerith  character  identification  header  that 

is  printed  as  a  banner  at  the  start  of  the  plot. 

SCALE  =  Multiplicative  factor  to  increase  or  decrease  the  plot  size 
(usually  =  1.0) 
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A2.2  CARD  2  XSIZE,  PFRBEG,  PFREND,  DELTAX,  ITYP,  IXAXIS,  NUMFIL 
XSIZE  =  Length  of  X-axis  in  inches 

PFRBEG  =  Beginning  wavenumber  on  plot  in  cm  *  or  wavelength  in  pm 
PFREND  =  Ending  wavenumber  on  plot  in  cm  ^  or  wavelength  in  pm 
DELTAX  =  The  interval  between  the  numbered  values  on  the  X-axis 

ITYP  =  0  Radiance  per  pm  vs  pm 

=  1  Radiance  per  cm'l  vs  cm-* 

=  2  Transmittance  vs  pm 
=  3  Transmittance  vs  cm-* 
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IXAXIS  =  0  X-axis  will  be  linear 

=  1  X-axis  will  be  logarithmic 

NUMFIL  =  0  uses  next  available  file  of  data 

>  0  uses  file  of  data  specified  by  NUMFIL 

A2.3  CARD  3  YSIZE,  YRMIN,  YRMAX,  DELTAY,  ICRV,  IYAXIS,  NMYDEC 
YSIZE  =  Length  of  Y-axis  in  inches 

YRMIN  =  Minimum  transmittance  or  radiance  value  to  be  plotted 
YRMAX  =  Maximum  transmittance  or  radiance  value  to  be  plotted 

In  a  log  plot  YRMIN  and  YRMAX  are  input  as  the  values  of  the  respective  expo¬ 
nents.  In  a  linear  plot  they  are  entered  as  the  actual  minimum  and  maximum 
value, 

DELTAY  =  The  interval  between  the  numbered  values  on  the  Y-axis 

ICRV  =  0  Normal  plot  of  line 

>  0  Calls  special  plotting  routine  to  plot  dashed  and  dotted  lines 
=  1  Solid  line  without  symbols 

=  2  Dashed  line  without  symbols 

=  3  Dotted  line  without  symbols 

=  4  Alternating  dashes  and  dots  without  symbols 

=  5  Alternating  dashes  and  2  dots  without  symbols 

=  6  to  10  Same  as  ICRV  =  1  to  5  with  symbols  at  every  point 

>  10  Alternating  dashes  of  different  lengths 

<  0  Data  points  only  using  symbol  number  =  |ICRV  | ,  where  sym 
bol  number  refers  to  the  computer -system  plotting  table 

IYAXIS  =  0  Y-axis  will  be  linear 

=  1  Y-axis  will  be  logarithmic 

NMYDEC  =  Number  of  digits  to  the  right  of  the  decimal  point  on  the  Y-axis 

A2.4  CARD  4  IRAD,  ITRP 
If  IEMSCT  =  0, 

IRAD  =  ITRP  =  0  plot  of  total  transmittance 

ITRP  =  1  Plot  of  water  vapoi’  band  transmittance 

=  2  Plot  of  uniformly  mixed  gases  transmittance 
=  3  Plot  of  ozone  transmittance 
=  4  Plot  of  nitrogen  continuum  transmittance 
=  5  Plot  of  water  vapor  continuum  transmittance 
=  6  Piot  of  molecular  scattering  transmittance 
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=  7  Plot  of  aerosol  transmittance 
=  8  Plot  of  nitric  acid  transmittance 
=  9  Plot  of  aerosol  absorption 

If  IEMSCT  =  1, 

IRAD  =  0  Plot  of  total  transmittance 

=  1  Plot  of  atmospheric  radiance  per  cm  1  vs  cm  1 

=  2  Plot  of  atmospheric  radiance  per  Mra  vs  Mm 

If  IEMSCT  =  2, 

IRAD  =  0  Plot  of  total  transmittance 

=  1  Plot  of  atmospheric  radiance  per  cm  *  vs  cm  * 

=  2  Plot  of  atmospheric  radiance  per  Mm  vs  Min 

■-  3  Plot  of  scattered  radiance  per  cm  ~  *  vs  cm ”  * 

=  4  Fiot  of  scattered  radiance  per  Mm  vs  Min 

=  5  Plot  of  reflected  radiance  per  cm'"*  vs  cm"* 

=  G  Plot  of  reflected  radiance  per  Mir.  vs  Mm 

-  7  Plot  of  total  radiance  per  cm”*  vs  cm”* 

-  8  Plot  of  total  radiance  per  Mm  vs  Mn> 

If  IEMSCT  =  3, 

IRAD  =  0  Plot  of  total  transmittance 

_  i 

=  1  Plot  of  transmitted  solar  irradiance  per  cm  1  vs  cm 
=  2  Plot  ox  transmitted  solar  irradiance  per  Mm  vs  Mm 
=  3  Plot  of  incident  solar  irradiance  per  cm”*  vs  cm"* 

=  4  Plot  of  incident  solar  irradiance  per  Mm  vs  Mm 


A2.5  CARD  5  ISAMFL,  ISAM  FT 

ISAMFL  and  ISAMFT  are  set  to  n-  lie  the  next  set  of  data  to  be  plotted. 
ISAMFL  is  used  when  the  user  wishes  to  plot  a  second,  third  or  ,  variable 
from  the  same  file  of  LOW TRAN  G  Tape  7  data.  The  user  should  be  sure  that 
the  range  of  the  variable  falls  within  the  range  as  specified  on  the  plot  axis  by 
CARD  2  and  CARD  3. 

ISAMFL  ~  o  Normal  advance  to  next  file  of  data 
=  1  Rewind  and  stay  on  same  file 

ISAMPT  -  0  Normal  advance  to  new  plot 

-  1  Plot  data  on  same  p/iysical  plot 


A3.  THE  AFGL  COMPUTER  SYSTEM  CALCOMP  PLOT  ROUTINES 


The  program  LOWPLT  uses  several  AFGL  computer-system  CALCOMP 
plot  calls  that  may  differ  from  the  system  plotting  routines  available  to  the  user 
It  is  anticipated  that  with  suitable  adaptation  of  the  CALCOMP  plot  calls  includ¬ 
ing  the  calls  to  PLTID3,  PLOT,  ENDPLT,  NUMBER,  SYMBOL,  and  LINE  plot¬ 
ting  can  be  accomplished  with  minimal  difficulty. 

The  general  functions  of  the  AFC'L  system  CALCOMP  plot  routines  are  as 
follows: 

1.  CALL  PLTID3  (PROG1D,  XMAX,  YMAX,  FACTOR,  IDC) 

This  subroutine  must  be  the  first  routine  called  as  it  initializes  the  plot. 

PROCTD  =  A  30-character  holleri.th  array  used  as  an  identifier  on  the 
plot  header 

XMAX  -  Max  dimension  in  X  inches  (X  limit  for  entire  plot)  (real) 

YMAX  =  Max  dimension  in  Y  inches  (Y  limit  for  entire  plot)  (real) 

FACTOR  =  A  multiplicative  factor  to  change  size  of  plotting  (usually  1.0) 
in  both  the  X  and  Y  directions 


IDC  Directs  plot  output  to  different  plot  queues 

=  1  -  all  plot  files  going  to  tapes  (user  specified) 

2  -  all  scratch -tes t- runs  (no  file  created) 

3  -  for  pens  -  standard  11 -in.  black  ballpoint  pen  plots 

4  -  for  pen  -  11 -in.  red  ink  pen 

5  .tor  pen  -  32 -in.  red  ink  pen 


2.  CALL  ENDPLT  -  PHYSICAL  END  OF  PROGRAM 

Subroutine  ENDPLT  must  be  the  last  plotting  subroutine  called  in  all  levels. 


ENDPLT  will  empty  plotting  buffers  and  write  identifying  information  at  the  end 
of  the  page.  A  call  to  ENDPLT  causes  termination  of  the  job. 

3.  CALI,  LINE  (X,  Y,  N,  K,  J,  L,  XMIIv,  Da,  YMIN,  DY,  SYMSZE) 

Subroutine  LINE  produces  a  single  line  by  connecting  the  points  defined  in 


the  dimensioned  variables  X  and  Y. 

X  =  Array  of  X  values  (real) 

Y  =  Array  of  Y  values  (real) 

N  -  Number  of  points  to  ue  plotted  (integer) 

K  =  Repeat  cycle  (usually  K  =  l)  (integer) 

When  K  =  2  the  first,  third,  fifth,  etc.  ,  points  will  be  plotted. 

When  K  -  3  the  first,  fourth,  seventh,  etc.,  points  will  be  plotted,  etc  ■ 

Example  CALI,  LINE  (X(  l),  X(2 ),  NV. . )  This  is  usually  used 

when  X  and  Y  arrays  are  one  mixed  array. 


J  =  Control  for  using  symbols  (integer) 

J  =  0  will  produce  a  line  plot  without  symbols 

J  =  1  will  produce  a  line  plot  with  a  symbol  at  every  point 

J  =  2  will  produce  a  line  plot  with  a  symbol  at  every  second 

point  etc.  ,  a  negative  J  will  suppress  the  lines  between 
the  points. 

L  =  A  number  describing  the  symbol  to  be  used.  Only  symbols  whose 

integer  equivalent  is  0  to  13  are  centered  around  the  point 
(X,Y).  For  all  others  the  point  (X,  Y)  is  at  the  lower  left 
corner  of  the  symbol. 

XMIN  =  Starting  value  of  X-axis,  in  units  of  X  (real) 

DX  =  The  difference  between  the  final  value  on  X-axis  and  the  starting 

value  of  X-axis  divided  by  the  length  of  X-axis  in  inches  (real). 

\'MIN  =  Starting  value  of  Y-axis,  in  units  of  Y  (real). 

DY  =  The  difference  between  the  final  value  on  Y-axis  and  the  starting 

value  of  Y-axis  divided  by  the  length  of  Y-axis  in  inches. 

SYMSZE  =  A  number  defining  the  size  of  the  symbol  to  be  used  at  the  point 
(X,  Y)  (real).  (Default  size  is  0.08  in.  ) 

4.  CALL  PLOT  (X,Y,IC) 

Subroutine  PLOT  is  used  to  move  the  pen  and  to  redefine  a  new  origin. 

X  =  X  coordinate,  in  inches  (real) 

Y  =  Y  coordinate,  in  inches  (real) 

1C  =  If  IC  =  2,  pen  down  as  pen  moves  to  (X,  Y)  (integer) 

If  IC  =  3,  pen  up  as  pen  moves  to  X,Y. 

If  IC  -  -2  or  -3  a  new  origin  is  defined  at  X,Y. 

Note:  The  pen  will  move  to  location  X,  Y,  on  the  page  in  all  eases.  X  and  Y 
are  defined  with  respect  to  the  previously  defined  origin.  Page  frame  limits 
should  be  considered  in  all  cases.  If  1C  ^  ±  2  or  3,  an  error  message  is  printed. 

5.  CALL  NUMBER  (X,  Y,  HGHT,  FPN,  THETA,  N) 

Subroutine  NUMBER  will  interpret  and  plot  a  real  (floating  point)  number. 

X  =  X  coordinate  of  lower  left-hana  corner  of  first  digit,  in  inches, 

relative  to  the  current  origin  (real) 

Y  =  Y  coordinate  of  lower  left-hand  corner  of  first  digit,  in  inches, 

relative  to  the  current  origin  (real) 

HGHT  *  Height  of  numbers  to  be  plotted,  in  inches  (real) 

FPN  =  Number  to  be  plotted  (real) 

THETA  =  Orientation  of  the  number  with  respect  to  the  X  axis,  counter¬ 
clockwise  in  degrees  (real) 

N  =  Number  of  digits  after  the  decimal  point  (integer),  N  =  -1  will 

suppress  the  decimal  point 

6.  CALL  SYMBOL  (X,  Y,  HGHT,  BCD,  THETA,  N) 

Subroutine  SYMBOL  will  draw  a  series  of  symbols  as  defined  in  the  symbol 
table  in  the  CALCOMP  instruction  manual. 
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X 


=  Coordinate  of  the  lower -left  corner  of  the  first  character,  in 
inches,  relative  to  the  current  defined  origin  (real) 

Y  =  Coordinate  of  the  lower -left  corner  of  the  first  character,  in 

inches,  relative  to  the  current  origin  (real) 

HGHT  =  The  height  of  the  characters,  in  inches  (real).  For  pen  plots  the 
width  of  each  character  will  be  equal  to  the  height. 

BCD  =  This  parameter  and  the  last  parameter  in  CALL  to  SYMBOL  (called 
N)  determine  the  type  of  annotation  the  routine  produces.  If 
BCD  is  the  text  to  be  used  as  annotation,  usually  BCD  or  A  type 
format,  the  characters  must  be  left-justified  and  contiguous  in 
a  single  variable,  in  an  array,  or  in  a  Hollerith  literal.  Param¬ 
eter  N  must  contain  the  number  of  characters  to  be  plotted.  If 
BCD  is  a  single  character  of  text,  the  text  must  be  right-justified 
and  parameter  N  =  0. 

THETA  =  The  angular  orientation  with  respect  to  the  X  axis  counterclockwise, 
degrees  (real) 

N  =  This  parameter  plus  parameter  BCD  determines  type  of  lettering/ 

symbols  produced  by  routine  SYMBOL. 

N  >  0  -  defines  character  count  in  array  BCD,  left-justified. 

N  =  0  -  defines  single  characters  to  be  plotted,  right-justified. 

N  <  0  (negative)  -  determines  the  condition  of  the  pen  in  the  move  from  its 
present  position  to  the  place  where  the  symbol  is  to  be  produced 

If  N  =  -1,  the  pen  is  up  during  the  move,  after  which  a  symbol  is  produced. 

If  N  =  -2  or  less,  the  pen  is  down  during  the  move,  after  which  a  symbol 
is  produced. 

Three  of  the  CALCOMP  plot  calls  from  the  AFGL  plotting  library  used  in 
LOWPLT,  are  not  consistent  with  the  standard  CALCOMP  versions.  These 
calls  and  their  possible  modifications  are  as  follows: 

1.  CALL  PLT1D3  -  should  be  replaced  with  CALL  PLOTS  (0,  0,  LDEV ) 
where  LDEV  is  the  plot  device  designation; 

2.  CALL  ENDPLT  -  should  be  replaced  with  a  plot  completion  call  appro¬ 
priate  to  the  user's  system;  and 

3.  CALL  LINE  -  a  program  modification  must  be  made  so  that  the  starting 
value  and  the  scaling  factor  immediately  follow  the  data  values  in  the  X -array 
and  the  Y-array  respectively. 


A4.  SAMPLE  PLOTS 


Five  plots  were  run  from  the  output  of  the  sample  cases  shown  in  Section  11. 
The  input  file  used  to  generate  these  plots  is  listed  in  Table  Al. 

The  plot  depicted  in  Figure  Al  is  a  plot  of  three  variables  from  CASE  1. 

The  continuous  line  plot  represents  the  atmospheric  radiance,  dashed  line  rep¬ 
resents  the  path  scattered  radiance,  and  dotted  line  represents  the  total  radiance. 
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Table  Ai.  Sample  Input  File  for  Plot  Program 


LW  ABREU 

0W6TEST 

PLOTS 

1  .00 

8.0 

4000.0 

450C.0 

50.0 

1 

0 

0 

S  .  0 

-11. 

-8.0 

1  .0 

1 

1 

1 

1 

0 

1 

1 

8.0 

4000. 0 

4500.0 

50 . 0 

1 

0 

1 

6.0 

-It. 

-8.0 

1  .0 

2 

1 

1 

3 

0 

1 

1 

8 . 0 

400C  .  0 

4500.0 

50 . 0 

1 

0 

1 

6.0 

-11. 

-8.0 

1 .0 

3 

1 

1 

7 

0 

0 

0 

8.0 

900.0 

1150.0 

50.0 

3 

0 

2 

6 . 0 

0.0 

1  .0 

0 . 1 

3 

0 

1 

0 

\ 

0 

0 

8.0 

900.0 

1150.0 

50.0 

3 

0 

3 

6.0 

0.0 

1 . 0 

0 . 1 

1 

0 

1 

0 

1 

1 

1 

0.0 

900. 0 

1150.0 

50.0 

3 

0 

3 

6.0 

0.0 

1 .0 

0 . 1 

2 

0 

1 

0 

7 

0 

0 

8 . 0 

900. 0 

1150.0 

50 . 0 

3 

c 

5 

6 . 0 

0.0 

1  .  0 

0 . 1 

5 

0 

1 

0 

5 

0 

0 

8.0 

6000.0 

7500 . 0 

150.0 

1 

r> 

V 

6 

e.o 

-11.0 

-5.0 

1  .0 

2 

1 

■» 

( 

1 

0 

1 

1 

8.0 

6000 . 0 

7500.0 

150.0 

1 

0 

6 

6.0 

-11.0 

-5.0 

1  .0 

4 

1 

i 

3 

0 

0 

0 

-1.0 

The  plot  in  Figure  A2  is  a  plot  of  the  water  vapor  transmittance  for  CASE  2, 
The  plots  in  Figure  A3  are  taken  from  CASE  3,  The  continuous  line  represents 
the  transmittance  due  to  water  vapor  and  the  dashed  line  represents  the  aerosol 
transmittance.  In  Figure  A4  the  water  vapor  continuum  transmittance  from 
CASE  5  is  plotted.  The  plot  in  Figure  A5  taken  from  CASE  6  shows  the  incident 
and  transmitted  solar  irradiance. 
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Figure  Al.  Three  Variables  are  Plotted  From  Sample  Case  1.  The 
continuous  line  represents  the  atmospheric  radiance,  the  dashed 
line  represents  the  path  scattered  radiance,  and  the  dotted  line  is 
the  total  radiance 
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Figure  A2.  The  Water  Vapor  Transmittance  From  Sample 
Case  2  is  Plotted 
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Figure  A3.  Two  Variables  From  Sample  Case  3  are  Shown, 
The  continuous  line  represents  the  transmittance  due  to 
water  vapor  and  the  dashed  line  represents  the  aerosol 
transmittance 
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Figure  A4.  The  Water  Vapor  Continuum  Transmittance  is 
Depicted  by  the  Dashed/2  Dots  Lane 
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Appendix  B 


I.OWTRAW  3  Filter  Function  Program,  LCVVPIL 

Bl.  INTRODUCTIOI' 

The  LOWTRAN  6  Filter  Function  Program,  LOWFIL,  is  primarily  designed 
to  calculate  the  effective  atmospheric  transmittance  that  would  be  measured  by  a 
filtered  transmissometer.  It  can  also  be  used  to  calculate  the  effective  radiance 
seen  by  a  radiometer  measuring  the  radiance  emitted  or  scattered  by  the  atmos¬ 
phere. 

This  program  is  written  as  an  independent  program  package  for  use  with  the 
LOWTRAN  6  code.  It  assumes  that  the  basic  mass  storage  file,  TAPE7,  written 
by  LOWTRAN,  is  available.  It  was  written  to  be  fairly  flexible  in  its  use, 
although  it  is  recognized  that  users  may  wish  to  modify  the  code  to  tailor  it 
more  specifically  for  their  own  applications. 

The  effective  average  atmospheric  transmittance  measured  by  u  transmis¬ 
someter  looking  at  a  blaokbody  source  is  given  by: 

/t(u)  F(u)  B(u,TBR)du 

r  =  — ~ - —  ,  (Bl) 

J  F(u)  BKTBB)  dv 


where 
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t(u)  is  the  spectral  transmittance  of  the  atmosphere, 

F(t')  is  the  total  combined  response  function  of  the  sensor,  including  the 
transmittance  of  any  filter,  the  effect  of  the  optical  components  and  the  detector 
response, 

B(k,  Tgg)  is  the  spectral  emissivity  of  the  blackbody  source, 

Tbb  is  the  temperature  of  the  blackbody,  and 

v  is  the  frequency  and  the  range  of  integration  includes  all  the  non-zero 
values  of  F(u). 

The  atmospheric  or  scattered  solar  radiance  measured  by  a  radiometer  is 
proportional  to: 


-f\  <* 


)  F  (u)  du 


where  I  ( v )  is  the  spectral  radiance.  In  general  when  using  Eq.  (B2)  to  determine 
the  radiance  measured  by  a  radiometer  it  would  be  necessary  to  include  in  F  (v), 
the  radiometer  angular  field -of-view,  calibration  constants,  and  all  other  fre¬ 
quency  independent  quantities  that  can  be  neglected  for  Eq.  (Bl),  since  they 
cancel  out.  The  effective  normalized  radiance  is  given  by: 


/’< » 


)  F  (u)  d. 


*Eff  ~  F 


where  F  =  the  maximum  value  of  F(i>).  When  used  with  radiance  output, 
m  aX 

LOWFIL  will  calculate  both  R  and  The  units  of  R  will  depend  on  the  units 

of  F(v),  the  sensor  response  function;  Ig^  will  be  in  W/cm^-ster. 

The  effective  filter  weighted  average  transmittances  for  the  different  atmos¬ 
pheric  components  used  in  LOWTRAN  are  also  computed,  (for  example,  the 
transmittance  due  to  ozone,  the  water  vapor  continuum,  aerosol  etc.).  It  should 
be  recognized,  that  while: 


r(y)  -  rH20(u)  r CQ2  . 


THN03  (u) 


in  general: 


7  t  rH20  tC02 


HN03 


where  the  average  transmittance  for  each  of  the  atmospheric  components  is 
found  by  evaluating  an  expression  analogous  to  Eq.  (Bl). 


In  addition  to  the  eight  average  component  transmittanoes,  the  wate;  con¬ 
tinuum  and  band-type  absorption  are  combined  into  a  total  water  transmittance. 


TWater^  TH20  Band^  rli20  cont.  ^ 


(B6) 


Also  the  other  gaseous  transm ittarcex  are  combined: 


TGas(t,)  =  tC02+(m)  •  Tozone(u)  *  TN2(u)  *  Tmol.  scat(u)  *  THN03(u)  (B7) 

These  combined  transmittance  components  have  proved  useful  in  developing 
simple  analytic  expressions  for  the  effective  transmittance  measured  by  a  system 
(for  example,  Shettle  and  Wise^) 


T  ~  T 


Water  '  TGas  ‘  TAerosol 


<B8) 


where  analytic  expressions  are  derived  separately  for  TWiter*  rGas*  and 

/T* 

Aerosol' 


B2.  COMPUTATIONAL  DETAILS  OF  THE  PROGRAM 

The  program  starts  by  reading  in  a  set  of  filter  functions,  or  systems 
response  functions.  There  can  be  up  to  15  different  filters  each  specified  at  up 
to  80  different  wavelengths  (these  limits  can  be  modified  by  adjusting  the  dimen¬ 
sion  statements)  in  each  set. 

The  program  is  written  so  that  the  systems  response  functions  or  filter 
function  F(u),  can  be  input  either  as  a  function  of  wavelength  (pm)  or  as  a  func¬ 
tion  of  wavenumber  (cm  ^).  If  they  are  input  as  a  function  of  wavelength,  the 
wavelengths  are  converted  to  wavenumbers  and  the  order  of  the  points  is  reversed 
to  be  in  ascending  order  in  wavenumber. 

Next,  a  specified  number  of  LQWTRAN  cases  are  read  in  from  TAPE7  one 
at  a  time.  For  each  LOW  TRAN  case,  the  different  filters  are  cycled  over  one- 
by-one.  The  subset  of  the  LOWTEAN  wavenumbers,  which  bracket  the  values 
for  which  the  filter  is  defined,  are  found  for  each  filter  in  turn.  Then  the  values 
of  the  filter  function  are  interpolated  to  this  subset  of  the  LOWTEAN 
wavenumbers. 


■  1.  Shettle,  E.  P.  ,  and  Wise,  J.O.  (1983)  Simple  Analytic  Expressions  for 
Atmospheric  Transmittance  (to  be  published). 
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Finally  the  integrals  in  Eqs.  (Bl)  through  (B3)  are  evaluated  using  the 
trapezoidal  rule.  Alter  a  set  of  filters  is  used  with  the  specified  number  of 
LOWTRAN  cases,  the  next  set  of  filter  functions  is  read  In,  and  repeated  until 
the  input  is  terminated. 

B3  USE  OF  THE  FILTER  FUNCTION  PROGRAM 

B3.1  General  Remarks 

The  filter  function  program  LOWFIL  utilizes  TAPE7,  the  standard  mass 
storage  output  file  as  written  by  the  LOW  TRAN  6  computer  code.  The  control 
cards  for  the  filter  are  read  from  TAPES,  which  in  the  listing  that  follows  is 
equated  to  the  INPUT  file. 

The  program  is  written  so  that  a  set  of  up  to  15  filter  functions  (system 
response  functions)  can  be  read  in  at  once,  and  used  with  one  or  several  consec¬ 
utive  LOWTRAN  calculations.  The  filter  functions  can  be  input  either  as  a  func¬ 
tion  of  wavelength  (pm)  or  as  a  function  of  wavenumber  (cm  *).  The  wavelengths 
(wavenumbers)  should  be  in  ascending  order,  but  can  be  at  uneven  intervals.  An 
option  is  included  to  rewind  the  LOW  TRAN  TAPE?  output  file,  so  more  than  15 
filter  functions  (system  response  functions)  can  be  used  or  the  same  filters  can 
be  re-used  with  a  different  blackbody  source  temperature,  It  should  be  empha¬ 
sized  that  the  spectral  range  of  the  LOWTRAN  calculations  should  include  the 
full  sensitive  range  of  the  systems  response  functions  to  be  used. 

B3.2  Control  Cards 

There  are  three  (3)  basic  control  cards: 

CARD  1  NF,  NEW,  IFT,  TEMP,  IPRINT,  NLOW 
(FORMAT  (315,  FlO.2,  215)) 

Repeat  CARD  2  and  CARD  3  "NF"  times,  where  NF  is  specified  on  CARD  1. 

CARD  2  IDFIL,  KODE,  IFWV,  NW 
(FORMAT  (2A10,  315)) 

CARD  3  (WAVE  (I),  FF  (I),  I  =  1,  NW) 

CARD  3  is  a  free  format  with  as  many  cards  as  needed  for  all  the  values  of 
WAVE  and  FF. 

B3.2.1  CARD  1  NF,  NEW,  IFT,  TEMP,  IPRINT,  NLOW 

NF  indicates  number  of  different  filters 

NF  >  0  Read  in  NF  filters  (NF  £  15) 

=  0  Use  preceding  filter  for  next  LOWTRAN  output. 

<  0  Stop  filter  program 
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NEW  is  an  option  to  reuse  the  previous  LOWTRAN  data  set  for  the  next  set 
of  filters. 

NEW  =  0,  No,  Read  next  NLOW  sets  of  LOWTRAN  data 
=  1,  Yes,  Rewind  the  LOWTRAN  TAPE  7  output  file 

IFT  is  an  option  to  enter  blackbody  temperature  of  source 
IFT  =  0  No  blackbody 

=  1  Fold  in  blackbody  emmissivity 

TEMP  =  blackbody  source  temperature  in  degrees  Kelvin 

IPRINT  controls  the  information  printed 

IPRINT  =  10,  Print  LOWTRAN  transmittances  and  results  below 

=  5,  Print  filter  function  with  blackbody  function  folded  in 

=  0,  Only  print  the  filter  weighted  transmittances 

NLOW  indicates  number  of  LOWTRAN  data  sets  to  use  with  this  set  of  filters 
B3.2.2  CARD  2  IDFIL,  KODE,  IFWV,  NW 

If  NF  is  set  to  zero  on  CARD  1,  a  new  input  set  of  CARDS  2  and  3  are  not 
read.  If  NF  is  >  zero  then  "NF"  input  sets  of  CARD  2  and  3  are  read. 

IDFIL  =  20  hollerith  character  identification  for  the  given  filter 

KODE  =  Filter  identification  number  {5  digits) 

IFWV  is  an  option  to  input  the  filter  function  either  vs  wavelength  or  vs 
wavenumber. 

IFWV  =  0  Wavelength 
=  1  Wavenumber 

NW  =  Number  of  wavelengths  or  wavenumbers  for  the  filter  (NW  £  80) 
B3.2.3  CARD  3  (WAVE  (I),  FF  (I),  I  =  1,  NW) 

The  variable  NW  on  input  card  2  gives  the  number  of  wavelength  or  (wave- 
number)  filter  function  pairs  necessary  on  CARD  3  or  multiples  of  CARD  3. 

These  should  be  in  ascending  order  by  wavelength  (or  wavenumber). 

WAVE  =  Wavelength  or  wavenumber,  depending  on  IFWV 
FF  =  Corresponding  filter  function 

B4.  SAMPLE  PROGRAM  OUTPUT 

This  section  presents  a  sample  run  of  the  Filter  Function  Program  to 
illustrate  its  use.  The  contents  of  TAPE7  are  assumed  to  be  those  produced 
by  running  LOWTRAN  with  the  test  cases  described  in  Section  11.  The  input 
control  cards  are  shown  in  Table  Bl. 
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Table  Bl.  Sample  Input  Data  for  Filter  Function  Program 


1 

1 

1  5700.00 

25 

1 

2.29- 

2.41 

MICRONS  1 

0 

13 

2.29  0. 

2  . 

30  0.2  2.31 

0.8  2 

32  1. 

2  2.33 

3.0 

2. 34 

17.5 

2.35  46 

.0  2. 

36  49.8  2.37 

50.7  2 

38  18.0  2 , 39 

8.1 

2.40 

6.7 

2.41  0. 

2 

0 

1  1133. 

03 

4 

9.12- 

10.57 

MICRONS  2 

0 

30 

9.123 

0 . 0 

9.173  0.0 

9  •  223 

.3 

9-273 

1  .0 

9-323 

1  .  B 

9.373 

3.3 

9.423 

7.8 

9.472  14.3 

9.523 

28 . 1 

9.573 

49 . 3 

9.623 

64.6 

9.673 

69.6 

9.723 

69. 1 

9. 773  66.0 

9.023 

65.2 

9.873 

67 . 1 

9.923 

72 . 9 

9.973 

77.7 

10.023 

73.9 

1C. 073  77.0 

10.123 

69.4 

10.173 

50 . 0 

10.223 

24.1 

10.273 

9. 1 

10- 323 

3.7 

10.373  2.6 

10-423 

1  .4 

10-473 

.5 

10.523 

0,0 

10 . 573 

0.0 

9.  *a- 

11.08 

MICRONS  3 

0 

33 

9.461 

0.0 

9.531  0.0 

9.581 

.0 

9.631 

.  1 

9.681 

1  .0 

9.731 

2.1 

9.731 

2 . 3 

9.831  5.2 

9.881 

9.4 

9-931 

17.1 

9-981 

32 . 1 

10.031 

51  .6 

10.081 

65.2 

1C. 131  70.2 

10-181 

71 .0 

10.231 

68.2 

10- 281 

65.8 

10.331 

67.8 

10.381 

72.0 

10.431  75.5 

10 .481 

77.8 

10.531 

77.3 

10.581 

71.2 

10.631 

55.0 

10.581 

2G  .  5 

10.731  13.0 

10.781 

6.  1 

10.831 

3,6 

10.881 

2.0 

10.931 

.7 

10.981 

.0 

11.031  0.0 

11.081 

0.0 

0.000 

-1 

It  should  be  noted  that  one  system  response  function  is  used  with  the  first 
LOWTRAN  test  case,  and  the  same  two  filter  functions  are  used  with  each  of  the 
last  four  LOWTRAN  cases.  The  first  case  uses  a  blackbody  source  temperature 
of  57  00°K  corresponding  to  the  sun,  since  this  is  a  solar  scattering  case.  For  the 
other  cases  the  blackbody  source  temperature  is  1133°K,  the  typical  operating 
source  temperature  of  the  AFGL/OPA  transmissometer  system.  The  resulting 
output  is  shown  in  Table  B2. 
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B5.  STRUCTURE  OF  THE  FILTER  FUNCTION  PROGRAM 

A  description  of  all  the  subroutines  making  up  the  Filter  Function  Program 
is  given  in  Table  B3,  in  the  order  in  which  they  are  called. 
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Table  B2.  Output  for  Filter  Function  Program  (Contd) 
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Table  B3.  Description  of  LOWTRAN  Filter  Program  Subroutines 


LOWFIL 

WAVEN 

BRACK 

INTLOG 

BIGFIL 

INTRAD 

BLKBDY 

INTGRT 

COMBT 


Main  driver  program.  Reads  control  cards  and  TAPE7  output 
from  LOWTRAN 

Changes  a  system  response  function  vs  wavelength  (Mm)  to  a 
response  function  vs  wavenumber  (cm"1) 

Finds  the  LOWTRAN  wavenumbers,  which  bracket  the  system 
response  function 

Interpolates  a  pair  ox  vectors,  F(I)  =  f(xj_),  and  X(I)  =  Xj, 

i  =  1,  2,  ...  N,  to  a  new  set  of  coordinates;  FNEW  (J)  =  i(x-) 
and  XNEW  (J)  =  x-,  j  *  1,  2,  . .  .  M  J 

Routine  to  find  the  maximum  value  of  the  system  r  esponse  function 

Integrates  the  emitted  or  scattered  radiance  from  LOWTRAN 
times  the  systems  response  function 

Weights  the  system  response  function  by  a  blackbody  radiance 

Finds  the  average  value  of  the  transmittance  from  LOWTRAN 
weighted  by  the  systems  response  function 

Finds  the  weighted  average  transmittance  due  to  water  vapor 
(both  band-type  and  continuum  combined)  and  all  the  other 
gases  combined  together 


Appendix  C 

Single  Scattering  Geometry 


The  single  scattering  algorithm  requires  the  scattering  angle  and  the  equiv¬ 
alent  absorber  amounts  for  the  primary  solar  paths  from  each  scattering  point 
along  the  line -of -sight.  This  calculation  involves  tracking  two  optical  paths  for 
each  scattering  point.  One  path  leg  extends  from  the  observer  to  the  scattering 
point,  and  the  other  from  the  scattering  point  toward  the  sun  (or  moon),  ending 
at  the  top  of  the  atmosphere  (nominally  100  km).  Each  path  will,  in  general,  be 
bent  by  refraction  but  each  remains  in  a  fixed  vertical  plane  containing  the  path 
endpoints  and  the  center  of  the  earth.  The  intersection  of  the  two  planes  at  each 
scattering'  point  defines  the  scattering  angle. 

The  optical  path  within  each  plane  can  be  described  in  a  spherical  coordinate 
system  in  terms  of  the  initial  altitude  and  the  zenith  angle  at  that  altitude.  The 
new  air  mass  module  (Section  2)  is  used  repeatedly:  first  to  trace  the  line-of- 
sight  path  from  the  observer  to  the  endpoint,  tracking  the  altitude,  zenith  angle, 
and  absorber  amounts  at  each  scattering  point.  Next,  it  is  called  at  each  scat¬ 
tering  point  to  trace  the  path  from  the  scattering  point  to  the  sun. 

Aside  from  the  in-plane  tracking  of  each  path  leg,  it  is  necessary  to  deter¬ 
mine  the  angle  of  intersection  of  the  two  vertical  planes,  which  is  also  the  azi¬ 
muthal  angle  separating  the  projections  of  the  two  path  legs  onto  the  earth's 
surface.  The  projection  of  a  path  segment  onto  the  earth  forms  part  of  a  great 
circle  on  the  earth,  since  a  vertical  plane  always  includes  the  earth  center.  The 
projected  path  segments  for  a  typical  geometrical  configuration  are  shown  in 
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Figure  Cl  from  a  perspective  above  the  earth's  surface  looking  down.  The  rela¬ 
tive  azimuth  angle  tf.-  between  the  projected  path  legs  at  the  projected  scattering 
r,  p 

point  is  found  by  first  computing  the  absolute  azimuth  angle  for  the  observer -to- 
subsolar-point  projection  OS.  The  absolute  line -of -sight  azimuth  </&  is  re¬ 
quired  as  an  input  to  the  code.  Within  the  program,  absolute  azimuth  angles  are 
measured  from  the  local  east  direction,  with  positive  angles  indicating  north-of- 
east.  in  this  discussion.  In  the  LOWTRAN  input  directives,  however,  absolute 
azimuths  are  specified  in  the  clockwise  sense  from  north,  that  is,  positive  east- 
of-north.  A  conversion  from  one  convention  to  the  other  occurs  within  the  code. 

The  relative  azimuth  angle  at  the  observer  i //  is  calculated  next.  A  trans- 

T#  O 

formation  of  the  coordinate  system  of  the  optical  path  is  then  made  to  facilitate 

the  calculation  of  the  relative  azimuth  angle  i b  at  each  scattering  point.  This 

r  >  P 

angle  allows  precise  calculation  of  the  scattering  angle  at  each  point. 


Figure  Cl.  Looking  Down  on  the  Scat¬ 
tering  Path.  All  points  are  projected 
on  the  earth  and  all  line  segments  are 
parts  of  great  circles.  The  relative 
azimuth  can  be  seen  to  vary  as  one 
moves  the  scattering  point  along  OP 


The  subsolar  point  and  the  projected  observer  position  together  define  the 
great  circle  which  satisfies  the  condition 

tan  t>  =  A  cos  <j>  +  E  sin  ,  (Cl) 

where  0  and  <!>  are  the  standard  latitude  and  longitude  angles  shown  in  Figure  C2 

and  A  and  B  are  constants  to  be  found.  The  case  where  the  subsolar  point  and 

observer  lie  on  the  same  longitude,  tan  $  -  tan  </>  ,  is  treated  separately  since  <f> 

s  o 
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Figure  C2.  Latitude  and  Lcngitude  Angles, 
Where  <t>  =  0  Passes  Through  the  Greenwich 
Time  Base 


is  fixed  in  this  case  and  0  \4‘)  is  indeterminate.  Except  for  this  special  case, 

A  and  B  are  determined  for  the  particular  great  circle  passing  through  these  two 
points  by  using  the  angular  surface  locations  of  the  subsolar  (#gi  4Q )  and  observer 
(6o.  <Pq)  points  to  obtain 


A 


tan#  sind  -  tan#  sin<£ 
s _ o _ o _ s. 

cos <t>  ^  sin4  -  co sind 

SO  OS 


(C2) 


ana 


B  = 


tan#  cos^  -  tan#  cos4 
_ Q _ s  s _ o_. 

cos0o  sin</>  -  ccs0  sin<£ 
so  OS 


(C3) 


The  absolute  azimuth  angle  i/.  ,  for  the  plane  containing  the  great  circle  seg- 

cl  t  o  s 

ment  OS  can  be  written  as 

ds 


tan  1 / 


a,  os  ds- 


N  _  - 1  d# 

coso  do 


(C4) 
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where  d<9/d<£  along  this  path  is  given  by 


dfl  _  B  cos4>  -  A  sin# 
d<t>  [1  +  tan2  6} 


<C5) 


The  relations  given  in  Eq.  (C4)  can  be  derived  from  considerations  based  on 
Figure  02.  Using  the  observer  point  location,  the  equations  above  are  combined 
to  obtain 


tan  ^ 


a,  os 


A  sin</>Q  ~  B  cos <j>Q 

f  1  +  tan2  6  ]  cosO 
o  o 


tan0Q  cosl^g  -  4>Q]  -  tandg 

[1  +  tan2  0  j  cos®  sin[0  -  ] 

o  o  so 


(CG) 


The  relative  azimuth  between  the  line -of-sight  OP  plane  and  the  OS  plane  is  then 
given  by 


0 


r,o 


,  os 


a,  op 


<C7) 


It  is  then  necessary  to  calculate  the  angle  Aq  subtended  at  the  earth  center 
by  radial  lines  from  the  sun  and  the  observer.  Using  the  dot  product  between 
position  vectors  on  the  unit  sphere,  it  can  be  shown  that 

cosA  =  cosO  cost  cos  0  cost  + 


cos9  sin$  cos  9  sin  t>  + 

s  s  00 


sind  sin  0 
s  o 


=  cosO  cosO  cos[0  -  4>  ]  +  sinO  sinO 
SO  so  so 


(C8) 


The  two  angles  0  and  Aq,  together  with  the  multi -path  altitude  and  zenith 
iniormation,  specify  the  three-dimensional  single  scattering  geometry  completely 
for  any  scattering  source  point.  However,  the  angle  1 /■  between  the  PS  and  OP 

planes  is  yet  to  be  obtained. 
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It  proves  convenient  to  make  a  coordinate  transformation  so  that  the  line -of  - 

sight  optical  path  is  more  easily  parameterized.  In  the  original  configuration, 

the  longitude  and  latitude  of  each  scattering  point  are  not  known  and  cumbersome 

to  calculate.  The  transformation  used  moves  the  subsolar  point  and  observer 

location  on  the  sphere  so  that  tne  line-of-sight  runs  east  from  an  observer  on  the 

equator  at  (O'  =  0,  =  0).  The  angles  ^  (short  notation  for  i b  )  and  A  are 

O  O  O  I  i  o  o 

left  invariant  by  the  transformation  so  that  the  relative  position  of  the  sun  in  the 
sky  is  the  same  for  the  observer.  No  physical  changes  or  approximations  are 
introduced  through  this  transformation.  As  one  can  see  from  Figure  C3,  by 
using  Eq.  (C5)  and  Eq.  (C8)  for  the  new  coordinates,  this  transformation  must 
satisfy  the  following  equations  for  d^  and  A^: 

tan  d1  =  tan  i/*  sin  <t>'  (C9) 

s  o  s 


and 


cos  A  =  cos  O'  cos  <t> '  .  (Cio) 

o  s  s 

These  conditions,  when  inverted,  define  the  latitude  6'  and  longitude  of  the 

s  s 

new  subsolar  point  in  terms  of  the  fixed  angles  tf/  and  A^,  that  is, 

tan  .</>',  =  -  tan  A  cos  <//  (Cn) 

a  O  O 


and 


sin  O'  =  sin  A  sin  i /  .  (C12) 

s  o  o 

A  line-of-sight  scattering  point  P  now  has  tne  new  latitude  and  longitude  coordi¬ 
nates  =  0,  where  j3  is  the  angle  subtended  at  earth  center  by  '"3  dial 

lines  to  P  and  the  observer  O. 

It  is  now  possible  to  compute  the  relative  azimuth  angle  between  two  path 
legs,  that  is,  between  the  vertical  plane  containing  PS  and  the  vertical  plane 
containing  OP.  r|  le  azimuth  angle  at  (0,  -p)  for  the  primary  solar  path  connect¬ 
ing  (0,  -p)  and  (O'  ,  <jd  )  is  given  by 

sinAQ  sim /. 

far  [  V-p  (p)]  cosP  sinA  cosi/  -  sin|3  cos  A 


(Cl  3) 


S'  ^.o  /l  ^>.p 
op  £ 

(0,-0) 


east  along  equator 


Figure  C3.  Looking  Down  on  the 
Transformed  Scattering  Geome  ¬ 
try.  The  line-of-sight  lies  on  the 
equator.  The  subsolar  point  at 
(0g  >  4>'s)  has  the  same  relative 
position  as  in  Figure  Cl 


It  is  also  necessary  to  calculate  the  angle  A(j3)  subtended  at  the  earth  center  by 
radial  lines  to  the  sun  and  the  scattering  point  (0,  -(3).  This  Is  given  by 


cos  [A ((3)]  =  cos  Aq  cos|3  -  sin(3  sinA^  eosi iQ 


(Cl  4) 


where  in  the  absence  of  refractive  bending  A(jj)  would  be  the  zenith  angle  ol  the 
sun  as  viewed  from  the  scattering  point.  In  the  presence  of  refraction  A (/l)  can 
be  used  as  an  initial  guess  for  the  actual  zenith  angle,  and  the  primary  solar 
path  can  be  tracked  iteratively  to  correct  for  refraction  and  obtain  the  correct 
zenith  angle.  The  scattering  point  altitude  Hl'(j3)  and  approximate  zenith  angle 
A(/j)  are  used  to  specify  the  starting  point  of  the  primary  solar  path  leg  from  P 
to  space. 

The  set  of  angles  that  are  discussed  above  is  sufficient  to  track  both  legs  of 
each  L -shaped  single  scattering'  path.  One  additional  angle,  the  included  angle 
for  single  scattering  y</3),  is  needed  for  the  phase  function  in  the  scattering 
source  term.  That  angle  is  given  by 


cos  y  =  n  •  n 

ps  op 


(CIS) 


where  each  path  direction  vector  refers  to  the  local  direction  at  the  scattering 
point.  In  terms  of  the  two  local  zenith  angles  and  at  the  scattering  point 
and  the  relative  azimuth  between  vertical  planes  </  (;T/,  the  included  angle  becomes 


cos  ItOD]  =  sinu  sinu  cos  [il-  (/>)]  +  nose  eoso 
'  '  ‘  ps  op  p  '  J  ps  op 


(Clfi) 
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The  scattering  angle  7  and  the  cumulative  absorber  and  scattering  amounts  are 
computed  lor  each  scattering  point,  which  contributes  to  the  intensity  in  the  sum 
given  by  Eq.  (38).  LOW  TRAN  specifically  uses  as  scattering  source  points  the 
intersections  01  the  line -of -sight  with  the  model  atmospheric  layer  boundaries. 


Appendix  D 

Development  of  the  Standard  LOWTRAN  Phase  Functions 


In  developing  the  set  of  standard  aerosol  phase  functions  stored  in  LOWTRAN. 
the  need  for  reasonable  accuracy  was  balanced  against  the  desire  to  keep  the 
memory  requirements  compatible  with  the  other  parts  of  the  LOWTRAN  code. 

The  phase  functions  for  each  of  the  LOWTRAN  aerosol  models  were  calculated 
at  27  different  wavelengths  between  0.2  and  40  hm.  These  frequencies  were 
chosen  so  that  all  the  significant  features  in  the  refractive  indices  of  the  aerosol 
particles  were  included  and  the  phase  functions  for  other  wavelengths  could  be 
found  by  interpolation. 

The  phase  functions  were  calculated  using  Mie  scattering  theory  for  spherical 
particles.  While  the  assumption  of  sphericity  is  valid  for  the  liquid  aerosols,  it 
is  not  for  the  dry  and  dustlike  particles.  However,  there  is  no  practical  method 
for  exact  calculation  of  scattering  by  nonspherical  particles  and  replacing  irreg¬ 
ularly  shaped  particles  by  equivalent  spheres,  generally  gives  reasonable 
(although  not  completely  correct)  results.  The  errors  introduced  by  this  approx¬ 
imation  are  less  than  the  other  uncertainties  in  the  properties  of  the  aerosol 
models  and  the  other  approximations  made  in  the  subsequent  use  of  the  phase 
functions. 

While  the  phase  functions  were  determined  for  all  the  models  at  27  wave¬ 
lengths  it  was  recognized  that  to  store  such  a  large  data  set  would  add  signifi¬ 
cantly  to  the  length  of  LOWTRAN.  This  data  set  consisted  of  7  02  phase  functions 
or  23,868  data  points  (20  models  X  'ln  wavelengths  X  34  angles),  which  exceeds 
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the  combined  length  of  all  data  subroutines  currently  in  LOWTRAN  6.  This  full 
set  of  phase  functions  is  provided  as  a  separate  data  file  on  the  LOWTRAN  tape 
obtained  from  the  National  Climatic  Center,  and  they  are  tabulated  by  Shettlc 
et  al.  ^  A  reduced  set  of  VO  phase  functions  was  chosen,  each  of  which  could  be 
used  to  approximate  a  number  of  phase  functions  from  the  full  set.  A  look-up 
table  is  used  with  a  cross-reference  from  each  of  the  702  pairs  of  26  aerosol 
models  and  27  wavelengths  to  one  of  the  7  0  standard  phase  functions. 

This  reduced  set  of  phase  functions  was  chosen  by  the  following  several 
step  process.  The  rms  difference  (&..)  and  the  correlation  (r-  )  between  all 
pairs  of  phase  functions  were  calculated: 


34 


■  £  (m  uyyi  -  m  [r,<V])2  .  i  ■  *. 


7  02 


(Dl) 


k=  .1 


rij 


£»*  '£wuFi(Sk)pj<<’k>  -£*kpi(ek'£»kpj<\: 

k _ k _ k _  k _ 

ai  * 


(D2) 


where 


4 

W,  =  - g- 

k  [1R(\)  +  P^)} 


ai =  <!Ewkpi2(V  •  E 

k  k 


wk  " 


'  YV  P.(0.  ) 

/Ls  k  i  k 

.  k 


1/2 


a.  - 
3 


Ewkpj2(V  •  Ewk  -  [  Ewr 

k  k  L  k 


■yv 


,  i  1/2 


The  logarithm  in  L'q.  (Dl)  and  the  weighting  factor,  w^,  in  Pep  (D2)  were 

both  used  to  reduce  the  effect  of  a  few  extreme  values,  since  in  some  cases  the 

values  of  the  phase  function  versus  angle  range  over  5  or  G  orders  of  magnitude 

The  value  of  6.^  can  be  shown  to  be  approximately  equal  to  the  average  relative 

difference  between  P.(0)  and  P.(0). 

^  3 

1.  Shettle,  P.P.,  Abrcu,  L.W.,  and  Moose,  1L  (1083)  Angular  Scattering  Prop¬ 
erties  of  the  Atmospheric  Aerosols,  AFGL-TK-83-  [folje  published). 
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These  two  measures  of  the  difference  between  the  phase  functions  were  used 
to  group  them  into  70  subsets  of  similar  scattering  functions.  Where  possible,  a 


criteria  of  6 ^  0.20  and  r . .  2:  0.  980  was  used  for  all  i  and  j  in  a  subset.  In  some 


cases  it  was  necessary  to  relax  this  criteria  to  6..S  0.40  and  r..  >  0.900.  For 

13  ij 

a  few  of  the  scattering  functions  it  was  not  possible  to  find  any  other  phase  func¬ 


tions  that  were  similar  even  with  the  weaker  criteria.  These  were  left  as  indi¬ 
vidual  phase  functions. 

The  mean  value  of  all  the  scattering  functions  within  each  of  the  70  subsets 
were  calculated  at  each  of  the  scattering  angles.  The  original  7  02  phase  func¬ 
tions  were  then  compared  with  these  70  resulting  "mean"  phase  functions  to 
determine  the  closest  match  using  Eqs.  <Di)  and  (D2).  This  resulted  in  some 
changes  in  the  composition  of  some  of  the  subsets.  The  mean  scattering  function 
was  then  calculated  for  each  of  these  7  0  modified  subsets  of  scattering  functions. 
These  final  phase  functions  were  used  as  the  set  of  standard  LOWTRAN  phase 
functions.  These  standard  LOWTRAN  scattering  functions  are  shown  in  Figures 
Dl  through  DlO. 

With  this  final  set,  50  percent  of  the  original  7  02  phase  functions  matched 


one  of  the  standard  phase  functions  within  a  criteria  of  6..  s  0.  10  and  r..  2:  0.995, 

k)  ij 

and  95  percent  met  a  criteria  of  6^  ^  0.20  and  m.  >  0.  975.  The  poorest  agree¬ 
ment  between  the  original  phase  functions  and  the  corresponding  member  of  the 


Figure  Dl.  Standard  LOWTRAN 
Phase  Functions  Numbers  1  to  7 


Figure  D2.  Standard  LOWTRAN 
Phase  Functions  Numbers  8  to  14 


SCATTERING  ANGLE 

Figure  D3.  Standard  L.OWTRAN 
Phase  Functions  Numbers  15  to  2  l 


Figure  D5,  Standard  LOWTRAN 
Phase  Functions  Numbers  2  9  to  35 


SCATTERING  ANGLE 


Figure  D4.  Standard  LOWTRAN 
Phase  Functions  Numbers  22  to  28 


Figure  D6.  Standard  LOWTRAN 
Phase  Functions  Numbers  36  to  42 
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SCATTERING  ANGLE 


Figure  D7.  Standard  LOWTRAN 
Phase  Functions  Numbers  43  to  49 


SCATTERING  ANGLE 

Figure  D9.  Standard  LOWTRAN 
Phase  Functions  Numbers  57  to  63 


SCATTERING  ANGLE 

Figure  D8.  Standard  LOWTRAN 
Phase  Functions  Numbers  50  to  56 


SCATTERING  ANGLE 

Figure  DlO.  Standard  LOWTRAN 
Phase  Functions  Numbers  64  to  70 
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Table  Dl.  Cross  Reference  Table.  Standard  LOWTRAN  phase  functions  for  each 
aerosol  model  and  wavelength 

(a)  Wavelength  range  0.  2  to  6.  0  Mm 


>0EL 

JO. 

MODEL  0, 

NAME 

.20  0 

,30  0 

.55  0, 

,69  1 

.06  1 

.54 

2.00  2, 

.50 

2.70 

3.00  3 

.20  3, 

,39  5, 

.00  6, 

.00 

1 

rural  rhoo 

3 

5 

4 

4 

4 

6 

8 

22 

21 

22 

26 

26 

26 

27 

2 

RURAL  RH70 

3 

5 

4 

4 

4 

6 

8 

22 

21 

21 

22 

26 

26 

27 

3 

RURAL  RH80 

3 

5 

5 

5 

19 

6 

8 

21 

7 

21 

22 

8 

26 

1 

4 

RURAL  RH99 

58 

58 

62 

62 

62 

63 

63 

63 

60 

64 

64 

64 

21 

70 

5 

MAR  IT  RHOO 

59 

11 

1 1 

1 1 

20 

20 

20 

20 

28 

16 

16 

16 

16 

37 

6 

MARIT  RH70 

9 

59 

11 

13 

13 

26 

26 

26 

27 

29 

46 

28 

28 

29 

7 

MAP  IT  RH80 

9 

9 

14 

14 

14 

15 

15 

15 

66 

65 

27 

26 

29 

27 

8 

MARIT  RH99 

57 

57 

69 

69 

69 

63 

68 

68 

61 

70 

60 

21 

7 

66 

9 

URBAN  RHOO 

2 

18 

18 

19 

19 

6 

22 

22 

22 

22 

22 

22 

29 

29 

10 

URBAN  RH70 

2 

3 

18 

16 

19 

6 

22 

22 

21 

22 

22 

22 

27 

27 

1  1 

URBAN  RHBO 

2 

3 

18 

18 

19 

6 

22 

21 

7 

22 

22 

22 

27 

27 

12 

URBAN  RH99 

58 

58 

62 

62 

62 

62 

63 

63 

60 

64 

64 

64 

21 

60 

13 

OCEAN  RHOO 

10 

59 

11 

1 1 

13 

20 

20 

28 

28 

16 

16 

16 

16 

37 

14 

OCEAN  RH70 

10 

10 

14 

14 

13 

13 

26 

26 

27 

29 

29 

2B 

20 

29 

15 

OCEAN  RHBO 

10 

69 

14 

14 

14 

15 

15 

7 

66 

65 

27 

26 

29 

1 

16 

OCEAN  RH99 

47 

57 

69 

69 

69 

68 

6B 

68 

61 

70 

60 

21 

7 

66 

17 

TROPO  RHOO 

29 

16 

32 

32 

36 

36 

23 

23 

23 

38 

38 

38 

25 

25 

lb 

TROP0  RH70 

29 

IS 

32 

32 

32 

36 

23 

23 

23 

38 

38 

38 

38 

25 

19 

TROPO  RH30 

29 

28 

37 

37 

32 

36 

36 

23 

23 

23 

23 

23 

38 

25 

20 

TROPO  RH99 

15 

26 

28 

28 

37 

37 

32 

32 

32 

23 

23 

23 

23 

38 

21 

STRT  H2504 

20 

20 

37 

37 

24 

23 

38 

25 

25 

35 

35 

35 

40 

40 

22 

AGED  VOLC 

7 

20 

16 

37 

32 

32 

24 

23 

23 

23 

38 

33 

25 

25 

23 

FRESH  VOLC 

17 

51 

13 

13 

20 

20 

20 

28 

37 

37 

37 

37 

32 

32 

24 

RAD  FOG 

47 

30 

55 

55 

55 

13 

15 

12 

42 

12 

1 

26 

44 

52 

25 

ADVEC  FOG 

48 

53 

31 

31 

31 

41 

41 

41 

49 

17 

17 

17 

56 

50 

26 

MET  OUST 

59 

59 

11 

11 

30 

20 

28 

28 

46 

40 

46 

46 

33 

33 
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Table  Di.  Cross  Reference  Table.  Standard  LOWTRAN  phase  functions  for  each 
aerosol  model  and  wavelength  (Contd) 


(b)  Wavelength  range  7.2  to  40.  0  Mm 


MODEL 

no. 

MODEL 

NAME 

7.2 

7  , 

8.7 

9.2 

10.0 

10-6 

12.5 

15.0 

17.2 

18.5 

21.3 

30.0 

o 

o 

V 

1 

RURAL  RH00 

29 

1 

34 

34 

34 

34 

33 

33 

36 

36 

36 

36 

23 

s 

RURAL  RH70 

29 

1 

34 

33 

33 

33 

46 

46 

34 

34 

36 

36 

23 

3 

RURAL  RH80 

27 

1 

46 

46 

46 

29 

29 

46 

34 

34 

34 

36 

23 

4 

RURAL  RH99 

65 

65 

65 

65 

66 

66 

65 

27 

29 

29 

46 

33 

34 

5 

MARIT  RH00 

37 

37 

32 

36 

32 

32 

32 

36 

23 

23 

23 

38 

25 

6 

MAriT  RH70 

28 

28 

37 

37 

37 

37 

34 

36 

36 

36 

23 

23 

38 

7 

MARIT  RH80 

29 

29 

29 

29 

29 

29 

46 

34 

36 

36 

36 

23 

23 

8 

MARIT  RH99 

1 

1 

1 

1 

66 

66 

1 

29 

46 

46 

33 

34 

36 

9 

URBAN  RHOO 

46 

29 

34 

34 

34 

34 

34 

34 

36 

36 

36 

23 

23 

10 

URBAN  RH7 0 

29 

27 

33 

33 

46 

46 

46 

33 

34 

34 

36 

36 

23 

1 1 

URBAN  nH80 

27 

27 

46 

29 

29 

29 

29 

46 

33 

34 

34 

36 

36 

12 

URBAN  RH99 

7 

65 

65 

65 

66 

67 

65 

7 

27 

27 

29 

46 

33 

13 

OCEAN  RHOO 

37 

37 

32 

36 

32 

32 

32 

36 

23 

23 

23 

38 

25 

14 

OCEAN  RH70 

28 

26 

37 

37 

37 

37 

33 

32 

36 

36 

23 

23 

38 

15 

OCEAN  RH80 

29 

29 

29 

29 

29 

29 

46 

34 

36 

36 

36 

23 

23 

16 

OCEAN  RH99 

1 

1 

1 

1 

66 

66 

1 

29 

46 

46 

33 

34 

36 

17 

TROPO  RHOO 

25 

25 

25 

25 

25 

25 

35 

35 

35 

35 

35 

40 

39 

18 

TROPO  RH70 

25 

25 

25 

25 

25 

25 

35 

35 

35 

35 

35 

40 

39 

19 

TROPO  RHBO 

25 

25 

25 

25 

25 

25 

35 

35 

35 

35 

35 

40 

40 

20 

TROPO  RH99 

36 

38 

38 

25 

25 

25 

25 

35 

35 

35 

35 

35 

40 

21 

STRT  H2S04 

39 

39 

39 

39 

39 

39 

39 

39 

39 

39 

39 

39 

39 

22 

AGED  VOLC 

35 

35 

35 

25 

35 

35 

40 

40 

40 

40 

40 

39 

39 

23 

fresh  volc 

24 

24 

23 

23 

36 

38 

38 

39 

38 

38 

38 

25 

35 

24 

rad  FOG 

44 

44 

45 

45 

45 

45 

24 

23 

38 

38 

38 

25 

35 

25 

ADVEC  FOG 

43 

43 

67 

61 

42 

54 

42 

12 

12 

12 

12 

1 

44 

26 

MET  DUST 

34 

34 

34 

34 

36 

36 

36 

36 

23 

23 

23 

38 

38 

i  no 


standard  scattering  functions  was  6..  =0.29  and  r..  =  0  940.  A  comparison  of 

^ J  ^  J 

the  original  phase  function  and  the  standard  function  used  for  this  case  is  shown 
in  Figure  Dll. 
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Figure  Dll.  Comparison  of  One  of  the 
Original  Aerosol  Model  Phase  Func¬ 
tions  With  the  Corresponding  Standard 
LOWTRAN  Phase  Function  (see  Table 
Dl),  Showing  the  Poorest  Agreement 

A  look-up  table  is  stored  in  LOWTRAN,  identifying  for  each  of  the  original 
phase  functions  (26  models  x  27  wavelengths),  which  standard  phase  functions 
should  be  used  (see  Table  Dl).  For  wavelengths  and/or  relative  humidities 
other  than  used  in  the  original  set  of  Mie  calculations,  LOWTRAN  will  inter¬ 
polate  between  the  appropriate  standard  phase  functions  as  indicated  by  the  table. 


